la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUT 


REPORT  1  Jfpiorsd  ior 

ntr  AO-A248  7B8  — 


2b.  DECLASSIFICATION /DOWNGRA 


Ttasao 


4.  PERFORMING  ORGANIZATION  RE 


m 


l^lpUMBEV 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S} 


8a.  NAME  OF  PERFORMING  0RGA|iUZAJ4ON«<^ 
University  of  Notre  Dame 

Department  of  Physics 


6c.  ADDRESS  {City,  State,  and  ZIP  Code) 
Notre  Dame,  Indiana  46556 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


.trOFFlCTSYMBOL  I  7a.  NAME  OF  MONITORING  ORGANIZATION 


(If  applicable) 


Office  of  Naval  Research 
Resident  Representative 


7b.  ADDRESS  (C/ty,  State,  and  ZIP  Code) 

536  S.  Clark  Street  Room  286 

Chicago,  Illinois  60605-1588 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

1114  SC/ONR  N  00014-89-.J-1136 


to.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT 
ELEMENT  NO.  ]  NO 


8c.  ADDRESS  (Cfty,  State,  and  ZIP  Code) 

800  N.  Quincy  Street 
Arlington,  Virginia  22217-5000 


1 1 .  TITLE  (Include  Security  Classification) 

Theory  of  Semiconducting  Superlattices  and  Microstructures 


12.  PERSONAL  /  'JTHOR(S) 
Dow,  John  D 


13a.  TYPE  OF  REPORT 
Final 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  WVEREO  |14.  DATE  OF  REPORT  (Year,  Month,  Day)  |l5  PAGE  COUNT 

FROM  88-10-1.5Td>l-10-14i  92-03-01  |  75 


17 

COSATI  CODtS  1 

FIELD 

GROUP 

SUB-GROUP  1 

18  SUBJECT  TERMS  (Continue  on  revene  if  necessary  and  identify  by  block  i  umber) 

Superlattices  Hjalmarson-Ffe.nkel  Theory 

Deep  level  theory  Mobius  transforms  Aroorphus  silicon 

ZnSe,  Gap _  Theory  of  dopants  Ti  Tc  Superconductors 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  final  report  of  work  summarizes  thj  theoretical  development  of  Semiconducting 
Superlattices  and  Microstructures.  These  include; 

The  first  theory  of  Hjalraarson-Frenkel  core  excitons  in  superlattices,  and 
applied  it  to  strained-layer  systems.  This  theory  will  be  useful  in  characterizinj 
phenomena  in  such  superlattices,  using  synchrotron  radiation. 

A  unified  theory  of  dopants  in  Il-VI  materials,  with  emphasis  on  ZnSe  and 
candidates  for  blue-green  lasers.  This  explains  why  Ga  is  no  longer  an  effective 
dopant  in  Znj.j^MnxSe  for  x^D.l,  with  a  theory  that  also  provides  a  natural 
explanation  of  why  ZnTe  is  almost  unique  among  the  II-VI  semiconductors  in  that 
it  is  naturally  p-type.  This  work  provides,  a  theoretical  framework  for  viewing 
and  understanding  the  recent  successes  at  3M  by  Haas  et  al.  in  obtaining  blue 
emission  from  ZnSe.  It  also  provides  strong  evidence  against  the  thirty-year-old 


20  DISTRIBUTION /AVAIUBILITY  OF  ABSTRACT 
(3  UNCLASSIFIED/UNLIMITED _ □  SAME  AS  RPT  □  OTIC  USERS 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

George  B.  Wright 


DD  FORM  1473,  84  MAR  83  APR  edition  may  be  used  until  exhausted 

All  other  editions  are  obsolete 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
□  DTic  USERS  Unclassified 


122b  TELEPHONE  (IrKlude  Area  Code) 

(703)696-4202 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


eUA  aiwrwit  KhtUnt  OMtat:  1Ma-a074M7 


#19  Abstract  cont4nue<i, 

.  e^i^at'fon  that  self-compensation  is  responsible. for  the  inability  to  make 
a  p-n  junctiolT  in  II-VI  semiconductors.  v 

-  A  general  theory  of  substitutional  impafi'tiea  in  superlat^ices  that  predicts 

anomalous  changes  of  dopant  character  layer  thicknessel',  and  establishes 

an  intellectual  framework  for  understan'ding  the  properties -bf  impurities  in 
superlattices.  • 

i  '■ 

-  The  mathematical  theory  of  Mobius  transforms  its  applications  to  physical 
inverse  problems.  This  is  an  important  new  area  of  theoretical  physics, 
especially  for  the  Navy:  with  such  transforms,  one  should  be  able  to 
detect  the  thermal  radiation  of  a  distant  object  and,  from  the  radiation, 
reconstruct  the  shape  and  temperature  profile  of  the  object. 

-  A  new  molecular  dynamics  of  III-V  semiconductors,  based  on  Sankey's  ideas 
for  Si.  This  method  solves  the  quantum  local  density  equations  very 
rapidly  for  the  forces  on  atoms,  and  then  moves  the  atomic  nuclei  according 
to  Newton's  equations.  The  result  is  a  program  that  a  computer  can  digest 
that  can  compute  realistic  trajectories  of  atoms  in  condensed  matter.  This 
will  allow  us  to  simulate  many  interesting  microscopic  processes,  such  as 
semiconductor  growth,  impurity  diffusion,  chemical  reactions,  etc. 

-  A  new,  relaxed-lattice  model  of  isoelectronic  traps  in  GaP,  explaining  some 
old  mysteries  of  the  data. 

-  A  theory  of  the  localized  states  in  the  HOMO-LUMO  gaps  of  macromolecules. 

The  theory  has  the  property  that  it  can  often  be  evaluated  approximately 
for  very  large  molecule,  whose  electronic  structure  is  poorly  known.  The 
localized  states  associated  with  frontier  orbitals  often  are  responsible 
for  the  critical  chemical  reactions  in  which  macromolecules  partake . 

-  A  new  picture  of  doping  amorphous  Si,  together  with  an  explanation  of  its 
band-tail  states.  This  picture  appears  to  agree  with  recent  data. 

-  An  addition  to  the  field  of  high-critical-temperature  superconductors  by 
computing  the  electronic  structures  of  La2CuO^  and  other  high-T^  materials. 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNfflCANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


OFFICE  OF  NAVAL  RESEARCH 


FINAL  REPORT 


for 


Contract  N00014-89-J-1136 


Task  No.  NR  616-027 


Theory  of  semiconducting  superlattices  and  microstructures 


John  D.  Dow 

Freimann  Professor  of  Physics 
Department  of  Physics 
University  of  Notre  Dame 
Notre  Dame,  Indiana  46556 

March  1,  1992 


Reproduction  in  whole,  or  in  part,  is  permitted  for  any  purpose  of  the  United 
States  Government. 


*  This  do>;v.ment  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited. 


92-*09820 


4  16  041 


page  2 


THEORY  OF  SEMICONDUCTING  SUPERLATTICES  AND  MICROSTRUCTURES 

John  D.  Dow 
Department  of  Physics 
Univsrsicy  of  Notre  Dame 
Notre  Dame,  Indiana  46556 

A.  Core  excitons  ir.  superlattices 

We  have  developed  the  first  theory  of  Hjalmarson-Frenke’  core 
excitons  in  superlattices,  and  applied  it  to  strained- layer  systems.  We 
expect  that  this  theory  will  be  useful  in  characterizing  phenomena  in 
such  superlattices,  using  synchrotron  radiation. 

B.  Theory  of  dopants  in  ZnSe  and  related  materials 

We  have  provided  a  unified  theory  of  dopants  in  II -VI  materials, 
with  emphasis  on  ZnSe  and  candidates  for  blue-green  lasers.  We  have 
explained  why  Ga  is  no  longer  an  effective  dopant  in  Znj^_j^Mn^Se  for 
x>0.1,  with  a  theory  that  also  provides  a  natural  explanation  of  why  ZnTe 
is  almost  unique  among  the  II-VI  semiconductors  in  that  it  is  naturally 
p-type.  This  work  provides,  we  believe,  a  theoretical  framework  for 
viewing  and  understanding  the  recent  successes  at  3M  by  Haase  et  al.  in 
obtaining  blue  emission  from  ZnSe.  It  also  provides  strong  evidence 
against  the  thirty-year-old  explanation  that  self -compensation  is 
responsible  for  the  inability  to  make  a  p-n  junction  in  II-VI 
semiconductors . 

£•  Theory  of  deep  levels  in  superlattices 

We  have  developed  a  general  theory  of  substitutional  impurities  in 
superlattices  that  predicts  anomalous  changes  of  dopant  character  with 
layer  thicknesses,  and  we  have  established  an  Intellectual  fiamework  for 
understanding  the  properties  of  impurities  in  superlattices. 

D.  Mobius  transforms 

We  have  developed  the  mathematical  theory  of  Mobius  transforms  and 
applied  it  to  physical  inverse  problems.  This  is  an  important  new  area 
of  theoretical  physics,  especially  for  the  Navy:  with  such  transforms, 
one  should  be  able  to  detect  the  thermal  radiation  of  a  distant  object 
and,  from  the  radiation,  reconstruct  the  shape  and  temperature  profile  of 
the  object. 


E.  Molecular  dynamics  of  semiconductors 

We  developed  a  new  molecular  dynamics  of  III-V  semiconductors,  based 
on  Sankey's  ideas  for  Si.  This  method  solves  the  quantum  local  density 
equations  very  rapidly  for  the  forces  on  atoms,  and  then  moves  the  atomic 
nuclei  according  to  Newton's  equations.  Tlie  result  is  a  program  that  a 
computer  can  digest  that  can  compute  realistic  trajectoiies  of  atoms  in 
condensed  matter.  This  will  allow  us  to  simulate  many  interesting 
microscopic  processes,  such  as  semiconductor  growth,  impurity  diffusion, 
chemical  reactions,  etc. 
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F.  Isoe) ectronic  traps  In  GaP 

We  have  developed  a  new,  relaxed- lattice  model  of  isoelectronic 
traps  in  GaP,  explaining  some  old  mysteries  of  the  data. 

G.  Macromolecules 

We  have  de'^'eloped  a  theory  of  the  localized  states  in  the  HOMO-LUMO 
gaps  of  macromole’cules .  The  theory  has  the  property  that  it  can  often  be 
evaluated  approximately  for  very  large  molecules,  whose  electronic 
structure  is  poorly  know.  The  localized  states  associated  with  frontier 
orbitals  often  are  responsible  for  *:he  critical  chemical  reactions  in 
which  macromolecules  partake. 


Doping  of  a-Si 

We  have  provided  a  new  picture  of  doping  of  amorphous  Si,  together 
with  an  explanation  of  its  band- tail  states.  This  picture  appears  to 
agree  with  recent  data. 


High-T^^  superconduc tors 

We  have  added  to  the  field  of  high- critical -temperature 
superconductors  by  computing  the  electronic  structures  of  La2CuO^  and 
other  high-T^  materials. 
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CORE  EXCITONS  IN  STRAINED-LAYER  SUPERLATTICES 


JOHN  D.  DOW  and  JUN  SHEN 

Department  of  Physics 
University  of  Notre  Dame 
Notre  Dame,  Indiana  46SS6  USA 

SHANG  YUAN  REN 

Department  of  Physics 
University  of  Notre  Dame,  Indiana,  USA 
and 

University  of  Science  and  Technology  of  China 
Hefei,  Anhui,  People's  Republic  of  China 


The  physics  ot  core  exdtoiu  in  semiconductors  is  reviewed,  with  emphasis  on  the  fact  that 
Hjaltnarson-Frenkel  '4eep'  core  exdtons  are  observed,  and  co-exist  with  Wannier-Mott  'shallow' 
exdtons  which  are  not  ".ormally  resolved  experimentally.  The  theory  of  Hjalmarson-Frenkel  exd¬ 
tons  is  extended  to  exdt.'sns  in  superlatuces,  and  the  Gaid  core  exdton  in  CaAs\-,P,/Ga? 
strained-layer  superlattices  ts  predicted  to  change  from  a  resonance  in  the  conduction  band  (with 
apparent  negadve  binding  ener'^)  to  a  bound  state  in  the  gap  (positive  binding  energy),  as  the 
GaAsi^^Pj  layer  thickness  decreasits. 


t.  Intsviduction 

Franco  Bassani  pioneered  the  theory  of  core  exciu  ns  and  was  one  of  the  first  theorists  to 
call  attention  to  the  fact  that  the  Si  Ip  core  exdton  exhibits  a  binding  energy  significantly 
different  from  the  shallow  donor  binding  energy  in  Si  [1].  Since  he  was  also  one  of  the 
first  theorists  interested  in  artifidal  superlattices,  we  shall  honor  him  by  discussing  the 
physics  of  core  excitons  in  superlattices. 

When  a  soft  x-ray  exdtes  an  atom  in  a  semiconductor,  iv  creates  a  core  hoie  of  very 
small  radius  and  an  electron  with  the  same  wave  vector  (see  i'ii^res  1  and  2).  Because 
the  core  orbital  has  such  a  small  radius,  it  does  not  overlap  the  corresponding  orbital  on 
an  adjacent  site,  and  so  the  core  band  is  fiat,  with  infinite  mass.  ,<%s  a  result,  the  hole  is 
immobile,  and  is  a  fixed  point  charge  around  which  the  dectron  can  'evolve.  One  expects 
the  dectron's  orbit  to  be  hydrogenic,  with  an  envelope  wave  junction  obeying  the 
effecrive-mass  Schrddinger  equation 

((  — A^/2m'')V^— e‘/<r)^r)  =  £’V'{r). 

Here  E  is  the  energy  of  the  orbiting  electron  with  respect  to  the  conduction  band  edge,  e 
is  the  didectric  constant  of  the  semiconductor,  and  m '  is  the  conduction  bancl*:.  effective 
mass  [2]  (we  have  assumed  an  isotropic,  non-degenerate  effective  ntass,  for  simplidty). 
The  ground  state  of  such  an  dectron  in  the  presence  of  its  hole  should  be  the  hyf,  ogenic 
Is  state,  with  a  binding  energy  rdative  to  the  conduction  band  minimum  of 

=  (13.6elO(m'/mo<^), 

439 

S.  Girlanda  tt  at.  (tdt,),  Prograi  on  EUaron  ProptnUs  of  Solids,  439-409. 
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SI  CoraExcItons 


Wav*  vector 

Figure  1.  Energy  versus  wave  vector:  illustration  of  a  Si2p  core  excitation  in  Si.  The  core  hole 
band  is  flat  (lower  heavy  line),  and  the  excited  electron  is  either  (i)  associated  with  the  conduction 
band  minima  nenr  X  for  the  shallow  Wannier>Mott  exciton  (dashed  line),  or  (ii)  in  the  antibonding 
s-lihe  A I  deep  level  resonant  with  the  conducuon  band  for  the  Hjalmarson-Frenket  core  exciton 
(upper  heavy  line). 

typically  tens  of  meV.  Such  a  correlated  electron-hole  pair  is  a  Wannier-Mott  exciton  (3), 

Analogous  impurity  levels  to  the  Wannier-Mott  core  excitons  have  been  thoroughly 
studied  and  were  thought  to  be  completely  understood  until  recently  (4);  the  shallow 
hydrogenic  donor  states  (2|.  The  electrons  of  these  states  obey  the  same  effective-mass 
Schrddinger  equation  as  the  elecuons  of  the  core  excitons,  and  exhibit  the  predicted 
binding  energies.  In  the  case  of  P  substituting  for  r  Si  atom  in  bulk  Si,  the  physics  of 
shallow  donor  levels  has  been  well-established  for  .’ecades.  Therefore  one  expected  that 
the  Si  2p  core  ex  iton  would  exhibit  the  same  binding  energy  as  the  P  donor.  It  did  not, 
being  in  many  experiments  two  order?:  of  •<tAgni  ude  larger:  0.1  to  0.9 eK  [5-9). 

The  failure  of  the  Si  Ip  core  exciton  biao'  .  neigy  '-•)  match  the  P  donor  binding 
energy  was  particularly  perplexing  because,  to  ..n  ad.jqu-w  ';;'proximaiion  (the  Z  + 1  rule 
(1,10)  or  the  optical  alchemy  approximation  !'l)),  trx;  ceve  hole  h.'.i  r.ero  radios  and  the 
same  charge  distribution  as  a  proton,  rn  the  case  nf  ;ore-excitea  Si,  the  excited  atom 
with  its  ‘proton’  and  electron  is  'P"  (similarly  core-excited  Ga  and  In  are  ‘Ge*  and  'Sn, 
respectively). 

Hence  in  the  case  of  the  Si  2p  core  exciton,  one  expected  an  absorption  spectrum  that 
reflected  the  spectrum  of  the  'P  donor’  impurity  on  Si,  with  its  small  hydrogenic 
effective-mass  binding  energy:  13.6eK  (w’/mo<^).  The  observation  of  much  larger  SiTp 
core  exciton  binding  energies  ranging  from  0.1  to  0.9  eK  [S]  became  known  as  tire  'St  core 
exciton  problem'  -  a  problem  whose  explanation  was  complicated  by  the  fact  that  the 
experimental  results  apparently  did  noi  agree.  Indeed,  the  experiments  calling  for  a  bind¬ 
ing  energy  almost  the  size  of  the  Si  band  gap  were  particularly  difficult  to  accept  in  the 
light  of  theories  existing  at  that  time.  Bassani  [I]  was  one  of  'he  first  theorists  to  recog¬ 
nize  that  these  observations  called  for  a  new  theory  of  core  excitons. 

The  differences  between  P  in  Si  and  a  core-excited  Si  atom.  Si’,  are  small:  (i)  because 
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GaAs  Core  Excitons 


Hjalmarson-Frenkel 

Wannier-Mott 


Wave  vector 

Figure  2.  Illustration  of  the  Ga  3<f  core  excitons  of  CaAs.  The  shallow  Wannier-Mott  state  lies  only 
ImtV  below  the  conduction  band  edge.  The  deep  Hjalmanon-Frenkel  slate  lies  OAeV  above  the 
edge  and  is  a  T;  resonance  (13). 


of  the  Franck-Condon  principle  1 12],  there  is  no  lattice  relaxation  surrounding  the  core 
hole,  whereas  the  lattice  certainly  relaxes  around  a  P  impurity  •  an  effect  that  has  always 
been  thought  to  be  negligible,  and  (ii)  the  differences  in  the  core  hole  and  proton  charge 
distributions  have  always  been  believed  to  be  negligible.  Hence  a  successful  theory  of  the 
Si  core  exciton  problem  would  have  either  relied  on  dynamical  effects  such  as  time- 
dependent  screening  or  explicitly  embraced  the  sameness  of  Si'  and  P  -  and  ascribed  the 
apparent  differences  in  binding  energies  to  the  fact  that  experiments  probe  different 
aspects  of  the  same  defect  (I).  The  latter  approach  now  provides  a  simple  and  natural 
solution  to  the  Si  core  exciton  problem,  but  had  been  rejected  by  many  theorists  in  the 
1970’s  because  it  implied  that  the  theii-current  understanding  of  P  in  Si  had  been  incom¬ 
plete. 


Z  Co-Elxistencc  of  ‘Deep’  and  ‘Shallow’  States 

The  many  apparently  contradictory  facts  about  the  Si  2p  exciton  can  be  understood  once 
one  recognizes  that  impurities  such  as  F  in  Si  have  both  'deep*  and  ‘shallow’  states  (4)  • 
and  that  there  are  corresponding  deep  or  ‘Hjalmarson-Frenker  [13]  and  shallow  or 
‘Wannier-Mott’  excitons  (3).  Infrared  and  transport  experiments  are  sensitive  to  the  .shal¬ 
low  impurity  levels  while  core  exciton  experiments  resolve  the  ‘deep’  Hjalmarson-Frenkel 
excitons  rather  than  the  shallow  Wannier-Mott  states.  For  P  in  bulk  Si  the  deep  level 
lies  slightly  above  the  conduction  band  edge,  so  that  the  corresponding  exciton  has  an 
apparent  negative  binding  energy  (4,14)  (see  Figure  1).  However,  if  the  Hjalmarson- 
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Frenkel  Si  7p  exciton  is  near  a  surface,  its  energy  is  rather  dramatically  penurbcd,  so  that 
its  apparent  binding  energy  can  be  0.8eF  JIS)  (sec  Figure  3). 

Si  (100)  -  (2  X 1)  exciton 


1.00 


5"  U.75 

3, 

>> 

o> 

u  0.50 


0,25 


0.00 

Layer  1  Layer  2  Bulk 

Figure  3.  Si2p  Hjalmtrson-Frenkel  core  exciton  energies  at  the  (2X 1)  asymmetric-dimer  (100)  sur¬ 
face  of  Si,  after  Ref.  (15|.  'Up'  {‘down’)  denotes  the  upper  (lower)  surface  site  (assuming  the  Si  oc¬ 
cupies  the  lower  half  plane),  and  results  are  given  for  the  surface  laser  (Layer  1)  and  the  upper¬ 
most  sub-surface  layer  (Layer  2).  These  levels  coalesce  intu  a  resonance  in  the  conducuon  band  for 
excitons  in  the  bulk  (layer  oo). 

Thus  the  explanation  of  the  wide  range  of  observed  exciton  energies  is  dtat  the  different 
measurements  had  different  surface  sensitivities,  and  probed  core  exciton  binding  energies 
in  Si,  ranging  from  slightly  negative  values  characteristic  of  the  bulk  to  large  positive 
values  (nearly  the  size  of  the  band  gap)  at  the  surface. 

A  central  element  of  this  picture  for  the  Si  Ip  core  exciton  was  the  prediction  that  the 
Hjalmarson-Frenkel  exciton  lies  slightly  above  the  conduction  band  edge  in  Si,  with  an 
apparent  negative  binding  energy  •  a  prediction  that  could  be  tested  by  measuring  its 
energy  in  Stj,(7e|„,  alloys  as  a  function  of  alloy  compo-sition  x.  The  theory  [14]  predicted 
that  the  Hjalmarson-Frenkel  state  would  descend  into  the  gap  and  assume  a  positive 
binding  energy  for  x«<0.27,  and  then  re-enter  the  conduction  band  for  larger  x.  The  com¬ 
position  x~tl.n  [16]  corresponds  to  a  aoss-over  in  the  conduction  band  minima  of 
Si„Ce\-x  being  near  the  (100)  A'-point  of  the  Bnllouin  zone  of  St  for  x  =  1  to  the 
(111)  f^point  of  Ge.  Bunker  et  al.  (16]  performed  an  experiment  to  test  these  ideas,  and 
concluded  that  the  picture  is  correct  for  the  Si'ip  core  exciton,  but  that  the  predicted 
energies  of  the  Hjalmarson-Frenkel  exciton  were  slightly  lower  than  observed  (see  Figure 
4).  The  resulting  Si  Ip  core  exciton  binding  en-rrgy  is  about  -0.06  eV. 

The  Hjalmarson-Frenkel  core  exciton  level  slightly  above  the  band  gap  in  bulk  Si 
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Alloy  CompotaioJi  x 

Figure  4.  Energies  (in  eV)  versus  ailoy  composition  x  in  Si^Gei-j,  alter  Ref.  [16|.  conducuon 
band  edges  Oigbt  solid  line),  (2)  Warmier-Mott  exciion  lei'els  (dashed  line),  and  experimental  con¬ 
clusions  (16)  (heavy  solid  line),  obtairted  by  analyzing  dati  with  an  adjusted  version  of  the  theory 
of  Ref.  (14). 

corresponds  to  a  ‘deep’  level  of  F  in  Si  that  is  resonant  with  the  conduction  band.  This 
deep  level  had  not  been  a  part  of  the  theoiy  of  the  P  impurity,  and  yet  was  demanded  by 
a  fact  that  was  well-known,  but  not  fully  appreciated:  The  central-cell  defect  potential  of 
a  P  substitudonal  impuriiy  in  Si  must  be  greater  than  AeV  deep,  because  the  difference 
between  P  and  Si  s-state  atomic  energies  is  4eK  In  a  solid,  it  is  impossible  to  have  a 
4eV  perturbation  potential  that  does  not  alter  the  electronic  structure  on  a  'ief'  energy 
scale  •  and  the  old  shallow  donor  tlteoiy  indicated  alterations  of  the  electronic  structure 
only  on  the  lens  of  meV  scale  of  the  donor  binding  energy.  Thus,  in  retrospect  (4),  P  in 
Si  must  have,  in  addition  to  the  shallow  donor  levels,  four  deep  levels  that  lie  above  the 
conduction  band  edge  and  are  rconances:  one  r-like  level  that  corresponds  to  the 
Hjalmarson-Frenkel  core  exciton  and  a  triply  degenerate  /t-like  level  at  higher  energy. 
These  four  levels  are  localized  in  space  and  arise  from  the  four  perturbed  bonds  fomied 
by  the  P  impurity  when  it  replaces  Si. 

The  same  theoretical  picture  of  Hjalmarson-Frenkel  core  excitons  in  Si  described 
Ga  2d  core  excitons  in  the  bulk  113,17,18)  (Figure  5)  and  Ga  2d  and  In  dd  excitons  at  sur¬ 
faces  (Figures  6  and  7)  (17,19,20)  in  GaAs,  GaSb,  JnAs,  InSb,  and  InP  -  a  dramatic  suc¬ 
cess  because  the  surface  exciton  energies  differ  markedly  from  those  in  the  bulk,  yet  gen¬ 
erally  agree  with  the  theory. 

The  major  difference  between  C-.c  Si  2p  core  exciton  on  the  one  hanri  and  the  Ga  2d 
and  In  Ad  excitons  on  the  other  is  symmetry.  Dipole  selection  rules  imply  that  tlic  core- 
exdtcd  Ga  and  In  generate  electrons  in  odd-parity  states  and  therefore  have  the  energies 
of  the  p-like  Tj  deep  levels  of  Ge  and  Sn,  respectively.  In  Si,  the  cote-cxdted  electron 
occupies  an  j-like  A  i  level. 
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Figure  5.  Comparison  of  measured  (18)  and  predicted  energies  (in  eV)  of  Hjalmarson-Frenkel 
Ga  3d  core  excitons  in  GoP,  GaSb,  anc  oc.is,  after  Ref.  (13).  The  zero  of  energy  is  the  conduction 
band  edge. 


3.  Core  Excitons  in  Superiattices 

Superiatticcs  arc  particularly  interesting  materials  for  studying  core  excitons  becatise  the 
superlatticc  band  edges  arc  very  sensitive  to  the  layer  thicknesses,  while  deep  impurity 
levels  are  not.  Therefore  a  C'  V  core  excitation  in  a  GaAs  quantum  well,  for  example, 
would  have  its  core  exciuti,.n  energy  (i.e.,  the  Ge  impurity  Tj  level)  remain  almost 
independent  of  the  width  of  the  GaAs  well  while  the  conduct:  in  band  edge  would 
dramatically  increase  in  energy  as  the  well  became  thinner  (due  to  quantum 
confinement).  As  a  result,  the  apparent  core  exciton  binding  energy  (i.e.,  the  energy  of  the 
Ge  impurity  level  with  respect  to  the  conduction  band  edge)  increases  markedly  as  the 
well-width  decreases.  The  physics  is  even  more  complex  when  the  superlatticc  is  not 
lattice-matched,  and  strain  is  important. 

In  this  paper  we  apply  tl.ese  ideas  to  Ga3cl  core  excitons  in  GaAsQiPot/CaP 
strained-layer  superiattices.  We  believe  that  these  arc  the  first  theoretical  results  for  this 
type  of  Hjalmarson-Frenkel  core  exciton  in  III-V  superiattices.  Our  goal  is  to  show  that  a 
Hjalmarson-Frenkel  Ga3d  core  exciton  at  a  site  near  the  center  of  a  CoA jo.6^o.4  layer 
can  emerge  from  the  conduction  band  into  the  fundamental  band  gap  of  the  superlattice 
as  the  layer  thickness  decreases.  That  is,  its  apparent  binding  energy  with  respect  to  the 
conduction  band  edge  changes  from  negative  to  positive  with  decreasing  layer  thickness 
(see  Figure  8). 

The  theory  is  identical  to  a  theory  for  a  Ge  impurity  at  a  central  Ca  site  of  a 
GaAsQiPoA^GaP  superlaitice.  The  theoretical  formalism  has  been  published  for  such  a 
defect  in  strain-free  GaAs/Alj,Cai.j,As  superiatticcs  (21-24),  and  can  be  modified  in  a 
straight-forward  manner  to  treat  strained-layer  superiattices.  The  strain  energy  is 

E,  =  (a’/64)(2y(6cii  +  12ci2l((<f.y-<fo)/‘^o)^'*'2;y*(ci)  -Ci2|(6tfy*)^) 

where  the  sums  are  over  the  bonds  (of  length  d,j)  connecting  atoms  t  and  j  and  over  the 
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Figure  6.  Experimental  and  theoretical  Go  2d  Hjalmarson-Frenlcel  core  exciton  energies  (propell¬ 
ers)  at  the  (110)  surfaces  of  GaAs,  GaSb,  and  GuP,  after  Ref.  (19].  The  valence  and  conduction 
band  edges  are  denoted  E,  and  £,.  The  horizontal  lines  denote  the  bottom  of  the  surface  state 
band. 

distontons  Qiji,  of  bond  angles  defined  by  atoms  ij,  and  k.  Here  cn,  c\2,  <fo>  ^  i^re 
the  elastic  constants,  perfect-crystal  bond  lengths,  and  perfect-crysi^  lattice  constant  of 
the  relevant  material.  There  are  three  types  of  contributions  to  the  strain  energy:  (i)  from 
GaP,  (ii)  from  GaAs\^j,P,,  and  (iii)  from  the  interfaces.  Three  parameters  define  the 
geometry  of  the  strained  superlattice:  (i)  dperp,  the  lattice  consuni  in  the  direction  per¬ 
pendicular  to  the  growth  or  z  direction  (this  is  the  same  in  both  CaAs\.,Px  and  CaP 
layers),  (ii)  c(GaP),  the  strained  lattice  consunt  in  the  z  or  growth  direction  in  a  GaP 
layer,  and  (lii)  c(GaAs\.,Px),  the  corresponding  lattice  constant  in  GaAs)-xPx.  The 
strain  energy  is  minimized  with  respect  to  these  three  parameters,  and  Operp.  c(GaP),  and 
c{GaAs].xPx)  are  determined.  Tncn  the  off-diagonal  (only)  matrix  elements  of  the  Ham¬ 
iltonian  are  scaled  according  to  the  rules  of  Harrison  (i.e.,  ar  d~^  (25))  and  Slater  and 
Koster  (i.e.,  the  appropriate  angle  dependences  (25,26)).  Finally  the  superlattice  elcctromc 
structure  and  Green's  functions  are  calculated,  in  the  usual  fashion  [21-24]. 

The  predictions  of  the  theory  arc  contained  in  the  p-liJce  Ti  'Ce  impurity’  states, 
because  the  3d  core  hole  is  associated  with  Tj  electron  sutes  as  a  consequence  of  the 
dipole  selection  rules  for  optical  absorption.  Ihe  predicted  band  gaps  for  [001] 
GaAsQiPQi/ GaP  superlatticcs,  are  expected  to  lie  slightly  below  (27,28)  the  expenmental 
gap,  but  within  about  ^O.leF'  of  it.  The  band  edge  of  the  bulk  materials,  both 
unstrained  and  strained  as  in  the  superlattice,  are  compared  with  the  supcrlattice  band 
edges  in  Figure  9  (we  have  assumed  a  valence  band  offset  for  a  strained 
GaAsoiPoi/ CaP  heterojunction  of  60%  of  0.33 eF  (29)).  The  corresponding  ‘binding 
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Figure  7.  Experiment*!  and  iheoreticai  In  4d  Hjalmaison-Frcnkel  core  exnlon  energies  (propellers) 
at  the  (IIO)  surfaces  of  InAs,  InSb,  and  JnP,  after  Ref.  (19). 


Figure  8.  Predicted  apparent  binding  energy  of  a  Ca  3d  core  exdton  p.-like  lexel  at  the  center  of  a 
Cadso*Po«  layer  in  an  A^XIO  GaAjotPo*/ CaP  superlattice  as  a  function  of  N.  Note  that  for 
smaU  (lu'ge)  N  the  binding  energy  is  positive  (negative). 

mergies'  of  the  Hjalmarson-Frcnkel  exdions  are  displayed  in  Figure  8  (28).  As  the 
GoAsoiPot  layer  becomes  thinner,  the  conduction  band'xigc  moves  up  in  energy  with 
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Band  Edga  Ralationa  In  (GaAt^,P^J/(GaP),  Supartattica 


Figure  9.  Band  edges  in  unstrained  and  strained  (as  in  the  superlattice)  GoAsotPo  t  and  GiiP,  and 
in  the  SXS  GaAst^t^Pat/ GaP  superlatdce. 


respect  to  the  core  exdton,  until  the  corc-exciton  level  descends  into  the  gap.  For  simpli¬ 
city  of  presentation  in  Figure  8,  we  have  plotted  only  the  />j-like  core  exciton  level  for 

an  exciton  far  from  an  interface.  In  the  superlattice.  and  especially  near  interfaces,  the  72 
levels  (which  are  degenerate  in  the  bulk)  split,  with  the  splittings  being  largest  near  the 
interfaces. 

Sit«>D«p*nd*nc«  of  Cort  Exciton  L*y*t» 


SUt 


Figure  10.  Predicted  site  dependence  of  the  Ca'id  core  exciton  levels  (in  eV)  in  a  5X5 
GoAscuPa*/ GaP  superlattice.  The  band  gap  of  the  superlatuce  is  denoted  £^p(5f.)  The  Ca 
atoms  ..re  at  odd-numbered  sites.  When  only  two  lines  are  plotted  at  a  site,  two  of  the  exdton 
levels  are  nearly  degenerate. 


Figure  10  illustrates  how  the  core  exciton  levels  should  vary  in  energy  for  a  5X5 
GaAsoiPo^/  CaP  superlattice,  as  the  position  of  the  core-excited  Ca  atom  is  changed.  Of 
course,  expieriments  measuring  such  excitons  should  detect  a  broadened  line  whose  shape 
reflects  the  distribution  rf  sites  at  which  the  excitons  arc  created  and  the  splittings  at 
each  site. 
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4.  Summary 

We  hope  that  this  work  will  stimulate  efforts  to  measure  core  excitons  in  superlattices 
and  to  show  that  suitable  manipulation  of  band  edges  in  superiattices  by  controlling  layer 
thicknesses  will  cause  Hjalmarson-Frcnkel  excitons  to  move  into  and  out  of  the  band 

gap- 
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A  theory  of  deep  impuniy  texel.'i  in  superlaiiices  is  outlined,  and  applied  to  ZnSe/Zn„5Mni,,i.Se  uith  the  substitutional  dopant 
Ga7„.  Ga  is  predicted  to  undergo  a  .shallow -deep  transition  at  a  function  of  in  bulk  Zn|_.Mn,Se.  and  so  Ga  dopes  ZnSe  n-type 
hut  prevents  Zn,),«Mno<.Se  from  being  doped  n-iype  even  by  mivdulaiion  doping.  In  ZnSe/Zno«Mn,|..Se  superlattices,  the  band 
edge.s  are  quite  sensitive  to  changes  in  the  layer  thicknesses,  but  the  deep  levels  are  not.  As  a  result,  shallow-deep  transitions  as 
functions  of  layer  (hickne.ss  are  predicted  to  occur.  The  physics  of  shallow-deep  transitions  in  superlattices  is  elucidated,  and  its 
relevance  to  the  Il-Vl  doping  problem  is  discussed. 


1.  Introduction 

Every  s-  and  p-bonded  substitutional  ‘impurity 
in  a  semiconductor  produces  four  “deep"  levels 
tnat  lie  near  dr  in  the  fundamental  band  gap  of 
the  host.  These  levels  are  due  to  the  central-cell 
defect  potential,  and  may  all  lie  resonant  with  the 
host  bands,  in  which  case  the  impurity  is  termed 
"shallow."  Or  at  least  one  of  these  levels  may  lie 
within  the  gap.  in  which  ca.se  the  impurity  is 
"deep"  (1).  Normally  one  such  deep  ie\el  is  .A,- 
symmetric  or  s-like  and  three  arc  o-like  (and  pos¬ 
sibly  degenerate,  depending  on  the  site  symmetr%'). 

2,  Sliailow-deep  n-ansitions 

In  a  \or)'  crude  (but  instructive!  appro.ximation 
(2).  the  deep  levels  are  insensitive  to  changes  of  the 
host  cor.’ position,  atomic  ordering  (e.g..  super- 
hmice  versus  random  alloy),  or  pressure,  and  re- 
min  their  absolute  energies.  In  contrast,  the  con¬ 
duction  and  valence  band  edges  are  sensitive  to 
such  changes,  and  so  it  is  rather  common  that  a 
band  edge  pa.sses  through  a  deep  level,  changing 
the  character  of  the  impurity  from  shallow  lo 
deep.  This  is  believed  to  be  the  case  [3]  for  S 
(Si  on  a  cation  site)  in  .Al.Ga,_,.As;  for  .v<6.2 
the  .A, -.symmetric  deep  level  lies  in  the  conduction 
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band,  making  Si  a  shallow  donor;  but  for  .v  >  0.3 
the  deep  level  is  in  the  fundamental  band  gap. 
allowing  the  Si  atom  to  trap  an  extra  electron 
rather  than  donate  one  to  the  conduction  band, 
rendering  the  material  semi-insulating  rather  than 
n-type  [3).  This  shallow-deep  transition  is  particu¬ 
larly  interesting  in  superlattices.  where  the  band 
edges  are  sensistive  to  the  choice  of  layer  thick¬ 
ness.  but  the  deep  levels  are  not  (2). 

3.  Shallow-deep  transitions  in  su)x;rlanices 

One  example  of  such  a  iran.sition  is  the  Ga^n 
impurity  near  the  cen'er  of  a  ZnSe  layer  in  a 
ZnSe/Znc jMn,><Se  [001)  .superlattice.  This  im¬ 
purity  is  a  shallow  donor,  with  its  .A, -symmetric 
deep  level  in  the  conduction  band,  for  iltick  ZnSe 
lovers.  Fin.  1  illustrates  how  the  conduction  band 
edge  of  3  .Vx  10  superlatiice  passes  through  tlio 
Ga  deep  level  as  .V  decrca.'cs  from  -V«10  to 
“  1  (a  single  layer  of  /’nSe).  !n  the  thin  .super- 
lattices  (for  A'  <  3.  according  to  the  theory  (2.4)). 
the  superlatticc's  conduction  band  edge  lies  above 
the  deep  level;  and  Ga  b'comes  a  deep  impuri’v: 
the  extra  electron  (relative  to  Zn)  of  neutral  Ga  is 
trapped  in  the  deep  level,  which  can  also  trap  an 
additional  election  of  oppo.viie  spin.  For  thick 
ZnSe  layers  (.V  >  3).  the  extra  v'.eciron  of  neutral 


so 
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Fig,  1.  Will-ccmer  Ca^n  A.i  deep  levels  in  (a)  thin  (I  x  10)  and  in  (b)  thick  (10  x  10)  quantum  t^'cll  ZnSe/Zn,ijMno  jSe  |(X)I1 
superlatiices.  The  thick  solid  lines  are  the  band  edges  of  ZnSe  and  Zn^  jhlnj,  jSe  bulk  semiconductors  respectively.  The  dashed  lines 
are  the  superlaiiice  h.ind  edges.  The  short  solid  lines  are  the  deep  levels.  The  G.a7,„  deep  level  is  above  the  superlaiiice 

conduction  band  edge  in  a  10  X  10  superlaiiice  and  is  in  the  gap  of  a  1  X  10  superlaiiice.  The  extra  valence  electron  in  the  Ga^n 
resonant  deep  level  will  fall  to  the  conduction  band  edge  in  the  10  x  10  superlaiiice.  while  in  ihc  1  x  10  superlaiiice  the  extra  electron 
will  occupy  the  d<-p  level  uhich  can  also  trap  another  electron  of  opposite  spin. 


Ga  spills  out  of  the  deep  level  (which  lies  above 
the  superlattice's  conduction  band  edge)  and  the 
Ga  is  autoionized.  creating  a  long-ranged  Coulomb 
potential  vvhich  binds  the  electron  at  zero  temper¬ 
ature  in  a  shallow  donor  level. 

Fig.  2  illustrates  how  the  Ga  deep  level,  the 
conduction  band  minimum  (CBM),  the  valence 
band  maximum  (VBM),  and  the  shallow  level  are 
predicted  to  var>'  with  ZnSc  layer  thickness  N  in 
an  A'  X  10  ZnSe/Zno.jMno.jSe  superlattice.  TTie 
predictions  use  an  empirical  tight-binding  Ham¬ 
iltonian  [5.6]  together  with  the  Green's  function 
method  (1). 

This  behavior  of  the  Ga  deep  level  as  a  func¬ 
tion  of  layer  thickness  N  is  similar  to  that  found 
as  a  function  of  alloy  composition  .v  in  Zn,.,- 
Mn.Se:  for  .v  >  0.1  the  Ga  deep  level  lies  in  the 
band  gap.  not  in  the  conduction  band  and  traps 
electrons  rather  than  donating  them.  This  means 
that  doping  of  Zn,.,Mn  .Se  for  .v  >  0.1  with  Ga 
should  produce  semi-insulating  rather  than  n-type 
material,  which  appears  to  be  the  case  c.xpcrimcn- 
tally  [7],  Even  modulation  doping  of  Zn,  .  ,Mn^Se 
with  Ga  will  not  produce  n-iype  material  for 
.V  >  0.1.  because  Ga  is  a  deep  trap  in  both  layers 
of  a  Zn|.,Mn,Se/Zn,.,Mn,,Sc  supcrl.atticc  for 
,v>.v>0.1. 


•  Fig.  3  illustrates  the  predicted  dependence  on 
alloy  composition  .v  of  the  levels  of  a  Ga2»„  impur¬ 
ity  in  the  ZnSe  layer  of  a  1  X  10  ZnSc/Zn,_*- 
Mn,Sc  superlattice.  For  .x  ■■  0,  the  superlattice 
reduces  to  bulk  ZnSe.  and  Ga  has  a  shallow 
hydrogenic  ground  state  donor  level  slightly  below 
the  conduction  band  minimum,  which  provides 
n-type  doping.  The  Ga  deep  level  lies  above  the 
conduction  band  minimum.  As  the  alloy  composi¬ 
tion  X  of  the  ZnSc/Zni_*.Mnj..Sc  superlattice  in¬ 
creases.  the  band  gap  opens  up  and  the  conduc¬ 
tion  band  edge  (measured  with  re.spcct  to  the 
valence  band  m,T.\:mum)  moves  to  higher  energy 
until,  near  .v  -  0,4,  the  b.and  edge  passes  through 
the  deep  level.  For  .v  >  0.4,  the  siabJe  ground  si  lic 
of  the  neutral  Ga  impuriiV  in  the  1  x  10 
ZnSe/Zn,  .  .Mn  ,Se  superlattice  has  the  deep  level 
occupied  by  one  electron.  This  deep  level  can  trap 
a  second  eicctrvjn  of  opposite  spin,  and  so  it 
removes  electrons  from  the  conduction  band, 
making  the  material  semi-insulaiing  rather  than 
n-iypc. 

Shallow-deep  transitions  can  occur  when  the 
valence  band  edge  passes  through  a  deep  level, 
much  the  same  as  when  a  conduction  band  edge 
docs.  The  vuicncc-band  shallow-deep  tran.siiions 
normally  have  a  much  more  dramatic  effect  on  the 
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C*Zn  'i’  {ZnS«)„/ {Zn,,Mn,,S#)„  Suptrliillc* 
Shallow-Deep  Translllon 


ZnS«  Liyir  Thlckntit  N 

Fig.  2.  Dependence  of  deep  levels  and  superlaitice  band  edges 
on  ZnSe  layer  ihicjiness  .V  (number  of  ZnSe  molecules  thick) 
for  A'xlO  ZnSe/ZnojMno jSe  superlaitices.  The  thick  lines 
are  the  superlaitice  conduction  band  (CBM)  and  valence  band 
(\'BM)  edges  respectively,  The  top  of  the  ZnSe  vafence  band 
(without  strain)  is  taken  to  be  the  zero  of  energy.  The  "thin  solid 
line' is  the  predicted  weli-cenicr  Oaj„  A)  deep  level,  The 
shallow  level,  which  follows  the  conduction  band  edge,  is 
depicted  by  a  dashed  line,  A  shallow-deep  transition  is  predic¬ 
ted  to  occur  around  ,V  ■  3, 


doping  character  of  a  material,  however,  because 
they  invariably  involve  p-like  deep  levels  capable 
of  containing  six  electrons  (whereas  the  conduc¬ 
tion-band  shallow-deep  transitions  normally  in¬ 
volve  .A  I -symmetric  levels  capable  of  trapping  only- 
two  electrons). 

We  believe  that  such  shallow-deep  transitions 
arc  responsible  for  the  different  doping  characters 
[8.9]  of  ZnSe  (which  can  be  easily  doped  n-type 
but  not  p-type)  and  ZnTc  (which  can  be  doped 
p-type):  deep  levels  that  lie  in  the  gap  of  ZnSe  and 
trap  holes  insie.i-J  lie  below  the  valence  band 
ma.ximum  in  ZnTe  and  donate  holes.  Clearly  one 
way  to  enhance  the  p-type  dopability  of  a  11- VI 
semiconductor  is  to  manipulate  the  semiconduc¬ 
tor's  valence  band  maximum,  moving  it  up  in 
energy  until  it  covers  the  deep  hole  traps.  For 
example,  the  p-dopability  of  CdTe  can  be  im¬ 
proved  [10]  by  fabricating  a  Cdle/ZnTc  strained- 
layer  siiperlattice.  In  this  case  the  strain  splits  the 


Dependence  on  x 
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X 

Fig.  3,  Predicted  dependence  on  Mn  foncenireimn  x  of  the 
Gaj„  deep  level,  and  the  5.hallo»  Jonor  level  in  1x10 
ZnSe/Zn,.,.Mn,Se  (001)  supcrlatiice>.  The  A)  deep 
level  it  resonari  with  the  conduction  nand  when  x  <  O,*!  (mak¬ 
ing  Ga  a  shallow  donor  impurity),  and  is  u  deep  trap  occupied 
by  the  extra  electron  for  x  >  0.4, 
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(a)  fb) 

Fig.  4.  Schematic  energy  band  siruciiire  (energy  in  e\'  versus 
w.-»ve  vector)  of  CdTe,  illusiruiing  liw  .strain  ((ualtoilvely 
changes  the  valence  hand  level  structure  with  respect  to  the 
deep  level  energy  and  covers  the  deep  level,  (a)  The  bulk 
semiconductor  vsith  a  deep  hole  trap  i :  ,:ii  alto  eontaini  at  least 
one  hole)  within  0.2  eV  of  the  v.ilence  b.md  edge.  The  Pi  .ind 
r„  bands  are  p-like  bands  tlia:  arc  spl.;  due  to  the  large 
spin-orbii  inieracn.in  in  CdTe,  'b)  A  2s4  superlaitice  has  an 
internal  strain  that  further  .splits  the  valence  band  and  cvwcrs 
up  the  deep  level,  autoioniring  the  hole. 
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valence  hand  ma.vinium  of  the  CdTc  and  covers 
up  deep  hole  traps  in  the  gap  slightly  above  the 
valence  band  maximum  (fig.  4).  A  more  complete 
di^eus.sion  of  this  p-doping  problem  will  be  pub- 
luhed  elsewhere  (1 1], 


4.  .SuinmaPi 

The  physics  of  shallow-deep  transitions  plays  a 
major  role  in  determining  the  doping  properties  of 
ll-\'l  semiconductors.  Band  edges  pass  through 
deep  levels  and  change  the  doping  character  of  the 
impurity  from  n-type  (donor)  to  semi-insulating 
(trap)  or  from  p-type  (acceptor)  to  semi-insulat¬ 
ing.  We  believe  that  by  better  understanding  and 
using  these  shallow-deep  transitions,  it  will  be 
possible  to  circumvent  many  of  the  doping  prob¬ 
lems  that  currently  plague  Il-Vl  .semiconductors. 
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The  physics  governing  deep  levels  in  superlattices  and  quantum  wells  is  elucidated, 
with  emphasis  on  the  importance  of  shallow-deep  transitions  caused  by  a  band  edge 
pissing  through  a  deep  level,  and  the  accompanying  change  in  doping  character  of  the 
impurity. 
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L  INTRODUCTION 

The  mcdeai  definition  of  a  deep  ener^  level  is  a 
level  thut  arigin-ates  from  the  central-cell  defect  po¬ 
tential  of  the  impurity  that  produces  it.'  By  this 
definition,  wliicii  supplants  the  older  energy  crite¬ 
rion  of  a  levc>l  within  the  fundamental  band  gap  by 
at  least  0.1  eV,  many  deep  levels  in  semiconductors 
lie  resonant  with  the  host  energ}*  bands  and  so  do 
not  fall  vrlthin  the  bexrd  gap.  Indeed,  for  s-^  and  p- 
bonded  substitutional  ..ip'orities  in  semiconductors, 
one  expects  four  dQ'ep  levu'?  near  or  within  the  fun- 
damenml  gap,  asso.•;:.ate.^i  with  the  four  perturbed 
impurity  bonds;  one  of  s-like  sjTnmetry  and  three 
of  p-like  symmetTN’.  A’'  impurity  is  termed  a  deep 
impurit}'  if  one  or  mors  of  mese  deep  levels  lies 
within  the  mndamentai  band  gap,  v.nd  so  can  crap 
an  electron  or  a  hole—rem;  ing  carriers  rather  tiian 
donating  them. 

In  this  paper  we  consider  aov'  dcip  levels  behave 
differently  in  superlattices  and  in  quanrom  wells 
from  in  bulk  semiconductors,  frnJ  that  the  ef¬ 
fects  of  a  superlaitice  on  the  ar  jc'utv,  energy’  of  a 
deep  level  tend  to  be  (i)  small,  s  viting  t  •  splitting 
the  level  by  an  amount  of  order  0..1  cV  or  iess,  and 
(ii)  in  most  cases  localized  to  wiu.ivi  abuui  three 
biatomic  laye-s  of  an  interiace  benv  ew  die  simer- 
lattice’s  conetituent  materials.  Therei'crc.  at  :.rst 
glance,  the  physics  of  deep  levels  in  superlamicar 
would  appear  to  be  uninteresting.  However,  the  ef¬ 
fect  of  the  superlattjcs  on  the  host  ba.id  edges  is 
quite  dramatic  for  small-period  superiattices,  pro¬ 
ducing  quantum  confinement  effeca  and  lypica.ly 
shifting  the  #uperia:*ice  band  edges  away  from  the 
large-period  superlaitice  edges  by  several  tenths  ci 
an  eV.  As  a  result  of  the  coniinement-ir.-hired 
movement  of  band  edges  in  superiattices.  these  fi  r.-s 
can  often  pass  through  a  deep  level,  convertinr  a 
shallow  impurity  into  a  dee?  impurity.  (Note  l  ie 
distinction  between  a  deep  imp.iri:y.  which  ha.',  it 
least  one  of  .ts  devp  levels  in  the  fundament.’'! 
and  a  deep  !c:  '  hich  may  or  may  not  be  w  .he 
gap.)  Such  shaJow-dee?  transitions  as  a  fur.ct..in  of 
decreasing  superlaitice  layer  thic.kness  are  ihe  luosc 
dramatic  superiactice  eiTects  on  deep  levels. 
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n.  SHALLOW-DEEP  TRANSITIONS  IN 
BULK  ALLOYS 

The  theory’  of  doping  in  semiconductors'  has  to  be 
revised  to  account  for  the  deep  levels,  which  play  a 
major  role  in  determining  the  doping  character  of 
an  impurity  and  its  stable  ground  state.  The  pre¬ 
vious  theory’  of  doping,  e:  .octive-mass  theory’,'  sim¬ 
ply  assumed  that  substitutional  impuritie.':  from 
Columns  zo  the  right  (left)  in  the  Periodic  Table  of 
the  host  atom  would  be  donors  (acceptors).  This  is 
not  always  the  case.  Fig.  1  shows  how  the  relative 
energies  of  the  dee?  levels  and  -he  band  edges  de¬ 
termine  the  character  of  a  dopant  that  by  the  old 
rules  should  be  a  single  donor,  P  in  Si.  In  most 
cases  considered,  the  s-  and  p-like  bonding  orbitals 
of  the  impurity  lie  well  below  the  valence  band 
maximum' and ’are  fully  occupied  by  electrons:  the 
antibonding  p-like  orbitals  are  irapty  and  high  in 
the  conduction  band.  The  case  o:*P  in  Si  corresponds 
to  the  conventional  shallow  donor  limit:  its  s-like 
dee?  level  lies  above  th?  conduction  band  minim^ 
(CBM),  and  the  extra  electron  of  neutral  P,  which 
would  occupy  that  level.  spiPs  out  and  falls  to  the 
conduction  band  edge,  icnliing  the  P.  The  e.xtra 
Coulomb  potential  of  P’  then  binds  the  electron  (at 
zero  temperature)  in  a  large-radius  hydrogenit 
shallow  impurity  level,  which  is  easily  ionized  ther¬ 
mally;  in  this  case  P  is  a  shallow  dono-  impurity 
(Fie.’  la).  Si  on  a  Ga  site  in  Ga.As  and  Ca  on  a  Zn 
site’  in  ZaSe  are  similarly  shallow  impurities. 

If  the  conduction  band  edge  of  Si  were  above  the 
deep  level  (Fig.  lb),  then  the  e.xtra  P  electron  would 
occupy  the  deep  level,  and  the  P  '•vouid  be  neutral 
jnst>''Jd  of  ionized.  The  deep  level  would  be  capable 
oi  n-appi.-:?  another  electron  of  opposite  spin  (con¬ 
vening  P'‘  ro  ?".•  or  another  \ich  icon*. trting  P“  to 
P*).  In  this  hypothetical  case  of  a  dee?  level  in  the 
srap.  P  would  be  a  deer  impurity.  Of  course,  P  in  Si 
Ts  no:  a  deep  impunv  •3--  3urJ-:er  et  al.'  have  shown 
that  in  Si.Ge-,.,,  there  is  a  range  of  alloy  composi¬ 
tions  .r.  such  that  the  conduction  band  edge  is  above 
the  ?  5-'iike  deep  level  ,ind  P  is  in.leeii  a  dee?  trap. 
Other  similar  deep  traps  are  ^xygen  on  a  P_^3ite  in 
GuP-  Si  on  3  Gu  site  .n  a.nd  (ja  on  a 

cation  s:te  in  Zn.-,. jSe. 

A  third  pos.sibiiity  •>.':iits.  namely  that  the  s-liKS 
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Fig.  1  —  Illustration  of  how  the  dee?  level  structure  relative  to  the  band  edges  detertaines  whether  an  inpurin’  is  (a)  sh^low,  (bi 
deep,  or  (ci  false  valence.  Electrons  (holes)  are  denoted  by  tilled  circles  (open  mangles).  It  is  assumed  that  the  iapunpf  is  from  one 
Column  to  the  right  in  the  Penodic  Table  of  the  atom  it  replaces,  such  as  P  in  Si.  Dee?  (shaJlowl  levels  are  solid  (dashed)  lines  and 
the  band  gap  is  smped.  CBM  and  \rBM  denote  the  conduction  band  minimum  and  the  valence  band  ma-ximum.  respectively. 


P  deep  level  might  lie  below  the  vaJence  band  edge 
of  Si  (Fig.  Ic).  In  this  case  the  hole  that  would  nat¬ 
urally  occupy  that  deep  level  were  it  in  the  gap  would 
instead  bubble  up  to  the  valence  bund  maximum: 
the  P  impurit}-  would  be  negativniy  charged,  and 
the  charged  impurity  would  bind  a  hole  in  the  re¬ 
sulting  shallow  acceptor  level.  In  this  hj^pothetical 
case,  we  would  say  that  the  P  impurity  has  a  false 
valence  of  minus  two  relative  to  its  normal  valence, 
in  that  it  would  act  as  though  it  came  from  Column 
in  of  the  Periodic  Table  instead  of  Column  V,  hav¬ 
ing  the  character  of  a  shallow  acceptor.  Of  course, 
the  resultant  P  deep  level  does  not  actually  lie  be¬ 
low  the  valence  band  maximum,  and  so  this  case — 
which  appears  to  be  quite  preposterous  to  persons 
who  thini:  in  terms  of  the  old  effective-mass  ideas — 
does  not  occur  for  P  in  Si.  But  it  does  occur  :*cr  In 
on  a  Te  site  in  Pbj_^n,Te:  In  is  well  known  ;a  be 
a  donor  in  PbTe  and  an  acceptor  in  SnTe — a  be- 
haxior  that  is  now  known  to  be  a  consequence  of  the 
small  band  gap  passing  through  a  p-like  deep  level 
as  a  function  of  alloy  composition''  (rather  than  la 
changing  its  site;.  As  a  rule,  false  valence  i.Tr.puri- 
ties  are  only  to  be  e.\-pected  in  smai;  band-gap  semi 
conductors  because  false  valence  is  associated  •-vit.'; 
a  level  moving  across  the  gap  and  is  most  eaoiiy 
achieved  if  the  gap  is  small. 


Shallow-deep  transitions  were  first  demonstrated 
as  a  function  of  alloy  composition  by  Wolford  and 
Streetman  for  the  N  isoeleca-onic  impurity  substi¬ 
tuting  for  P  and  As  in  GaL-\.s,.,P.  alloys  (Fig.  2).* 
Their  system,- tic  studies  on  carefuily  ion-implanted 
alloys  showed  that  the  Column  VI  dopants  S  and  Se 
produced  hydrogenic  ejaective-mass  donor  levels  that 
followed  the  band  edges  as  a  function  of  ailoy  com¬ 
position  X.  but  that  osygen  produced  a  deep  level 
that  was  unattached  to  any  band  edge  and  varied 
linearly  with  x.  Column-V  nitrogen  behaved  like 
Column- VT  oxy’gen.  as  a  functior.  of  x.  although  its 
energy  level  in  Gap  is  shallower  energetically  than 
S’s  or  Se’s,  and  for  x  <  0.2  the  .V  level  disappears 
into  the  conduction  band.  These  data  showed  that 
the  N  level  is  actually  a  deep  level,  similar  to  ox¬ 
ygen's  and  distinct  mom  the  si  illow  levels,  and  that 
the  deep  level  passes  into  the  conduction  band  rather 
easily.  (For  simplicity,  we  call  this  a  shallow-deep 
transition  even  though  the  .V.  being  isoelectronic  to 
As  and  P.  does  not  produce  a  shallow  donor  level.) 
Tnsy  also  indicated  that  S  and  Se  must  also  have 
deep  levels  similar  to  those  of  .V  and  oxygen,  lying 
above  the  conduct!  vn  ba.-.d  m.ir.iinum.  T .nally,  since 
N  is  isoelectronic  to  P  and  As.  the  aeiect  potential 
responr-ibie  for  producing  the  deep  level  is.  to  a  good 
appro.ximation.  com'ineu  to  the  central  ceil — and  so 
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Fig.  2  —  Schtsiade  summary  of  the  energy  levels  of  GaAS).^. 
with  N,  0,  S.  and  Se  aiuoa*subttitutional  impurities,  versos  al* 
loy  composition,  as  determined  by  Wolford  «t  al.*  The  sero  of  en¬ 
ergy  is  the  valence  band  mariraum  T  and  X  label  the  conduction 
C,  band  minima  ate  ■■^'and'i  ••  (2r/ct)(1.0.0)  of  Brillouin 
rone,  respectively.  Thr'ihin  Solid  lines  iabekd  N  are  the 
measured  nitrogen  and  oxygen  deep  ievels,  and  thv  oed  lines 
are  the  S  and  Se  shallow  levels. 


the  data  demonstrated  that  deep  levels  are  due  to 
the  central-cell  potential,  gi'.nng  rise  to  the  modem 
definition  of  a  deep  level.  All  of  these  ideas  were 
revolutionar.'  and  controversial  at  the  time  they  were 
first  proposed,  but  are  now  well-established.  Some 
other  impunties  that  exhibit  shallow-deep  transi¬ 
tions  as  functions  of  alloy  composition  are  Si  on  a 
Ga  site  in  and  Ga  on  a  cation  site  in 

Zni_;rMn,Se^  iSee  Fig.  3).  For  small  x,  these  impur¬ 
ities  are  donors,  causing  the  host  semiconductor  to 
be  rt-tjTie.  For  larger  z  they  become  deep  traps,  in¬ 
hibiting  conductivity,  and  producing  semi-insulat¬ 
ing  behavior. 


m.  SIL\JLLOW-DEEP  TRA.NSITIO.NS 
VERSUS  L.AYER  THICKNESS  IN 
STJPERLAITICES 

In  the  alloy  hosts,  the  band  ed^es  var;.'  with  alloy 
compost  ion  and  pass  through  the  deep  levels,  lead- 
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rig.  3  —  Schematic  energy  level  diagram  of  Gat,  in  (a)  ZnSe  and 
(b)  Zco»Ms.'  illustrating  the  s-like  and  p-like  deep  levels.  Holes 
(elecrronsi  are  denoted  by  open  triangles  (dosed  circlesi.  The  lowest 
shallow  level  is  cashed.  In  ZnSe  the  s-like  deep  level  of  Ga  :s 
above  the  conduction  band  minimum  tCBMi.  causing  its  electron 
to  be  aumiomaed  and  napped  in  the  resulting  shallow  level  I'a: 
sero  temperature I.  In  diis  case.  Ga  is  a  shallow  donor  impunty, 
maidng  21nSe  n-type.  The  s-like  levei  L  in  the  gap  for  Znt>Mn4  tibe. 
is  a  deep  trap  (for  either  an  electron  or  a  bole),  and  makes  tha 
materiai  semi-insulating. 


ing  to  shallow-deep  rransitions,  which  are  now  well- 
documented. 

Analogous  transitions  occm*  in  superlatticet  as 
functions  of  the  layer  thicknesses.’ ■*“**  To  under¬ 
stand  this,  consider  Si  in  the  middle  of  a  GaAs  uiyer 
in  a  GaAs/Afo.7Gao,3As  superb  ctice,  and  recognize 
that  a  Si  on  a  Ga  site  in  GaAs  is  a  shallow  donor 
impurit}',  but  that  Si  on  an  AC  site  in  AfAs  is  a 
deep  trap.  As  an  initial  approximation,  think  of  the 
absolute  energy  of  the  Si  deep  level  as  being  inde- 
pendenr.  of  the  GaAs  layer  th.o.kness. 

If  the  GaAs  layer  is  thick,  then  the  band  gap  3f 
the  superlattice  is  approximately  equal  to  the  band 
gap  of  GaAs:  if  the  layer  is  thin,  then  the  superlat¬ 
tice  band  gap  is  almost  the  AC  As  gap.*^  Hence,  by 
reducing  the  thickness  of  the  GaAs  layers,  it  is  pos¬ 
sible  to  move  the  superlattice  conduction  band  edge 
from  the  energj’  of  the  bulk  GaAs  conduction  edge 
to  that  of  AfAs.  This  is  a  quantum  confinement  ef¬ 
fect.  Since  the  s-!i;-;e  deep  level  of  Si  is  relatively 
consumt  in  energy,  the  superiattice  band  edge  moves 
up  through  the  Si  deep  level  as  the  GaAs  layer-wic-h 
decreases — and  Si  is  predicted  to  underg:o  a  shal¬ 
low-deep  transition.  This  is  illustrated  in  Fig.  4 
for  the  thick  GaAs  wells  in  an  13  x  18  GaAs/ 
Afo.TGa<,,j.As  superiattice  and  for  thin  GaAs  wells  in 
a  2x34  superiattice.  In  thin  Ga-As  quantum  wells. 
Si  on  a  Ga  site  is  predicted  to  be  a  deep  trap. 

Figures  5  and  6  show  how  relatively  insensitive 
the  deep  levels  in  a  superlaKice  are  to  the  super- 
lattice  ordering,  in  comparison  with  the  band  edges. 
Note  that  the  Si  s-like  level  shifts  only  slightly  near 
an  iniertace  (<0.1  eV)  and  assumes  its  bulk  vaiue 
near  the  center  of  a  lay".  The  superiattice  effects 
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Fig,  4  —  Schematic  illustration  of  the  shallow-deep  transition  as  a  function  of  layer  thickness  for  a  SL-,.  i.T.ouncy  near  the  center  of 
a  GaAs  quantum  well  m  a  (0011  GaAs/Alo  rOa^jAs  superlatiice.'  The  bulk  energy  band  edges  vs  position  in  the  (OOll  direction  axe 
denowd  by  thick  solid  lines.  The  superlanicc  band  edges  are  denoted  by  dashed  lines.  For  the  18  x  18  suoerlattice,  the  well  is  IS  hi- 
atomic  layers  thick,  and  the  Si  deep  level  lies  above  the  suDeriattice's  conduction  band  minimum  (CBM;  dashed  line; — causing  Si  to 
be  a  shallow  donor.  In  the  2  x  34  superlactice,  the  s-Iike  deep  level  is  beiow  the  CBM  of  the  superlattice,  a.-.d  Si  is  a  deep  crap.  .Vote 
that  the  CBM  of  the  superlactice  is  the  first  confinement  level  in  a  Kronig-Penney  model. 


have  virtually  disappeared  for  impurities  more  than 
three  bi-aiomic  layers  from  an  interface.  (See  Fie. 
5.) 

The  p-like  levels  have  similar  behavior,  splitting 
near  an  interface  (of  order  0.1  eV)  with  the  highest- 
energy  superlaaice  state  of  a  cation  vacancy  in  3 
X  10  GaAsMfo.:Ga,.;As  corresponding  to  the  p-like 


orbital  oriented  toward  the  Ga.As  and  the  lowest 
being  drected  toward  the  ACo-Csj^^As.  (See  Fig.  6.) 
The  variations  in  the  deep  level  energies  that  do 
occur  are  caused  by  the  greater  clectropositirfty  of 
the  A€  and  the  valence  band  offset.  For  deep  levels 
near  the  valence  band  ma-dmum  of  ihe  superlat¬ 
tice,  the  superiattice-induced  splittings  of  the  va- 
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Fig,  6  ~  PredicMd  energy  levels  of  a  cation  vacancj-  in  bulk  GaAs.  at  various  sites  /3  in  a  10  x  10  (0011  GaAs/Afi:Ga,jAs  superlattice, 
and  in  bulk_Afs-Ga.]As,  after  Ref.  7.  Note  the  significant  splitting  of  the  valence  band  ma.xt.nuro  in  the  superlattict.  and  the  corre¬ 
sponding  splitting  of  the  p-like  Tj-derived  deep  levels,  even  when  the  vacancy  is  distant  iota  the  interface  i  because  the  host  spectrsJ 
density  of  the  superlattice  is  likewise  split). 


lence  band  and  the  host  spectral  density  are  re¬ 
flected  in  the  deep  levels,  even  for  levels  distant  from 
an  interface. 

In  contrast  to  the  deep  levels,  the  band  edge.''  shift 
and  split  considerably  from  their  bulk  energies  in 
the  superlattice. 


IV.  APPUCATION  TO  THE  H-VI  DOPING 
PROBLEM 

These  ideas  about  deep  levels  can  be  used  to  cir¬ 
cumvent  the  n-VI  doring  problem.  Stated  simply, 
the  n-VI  semiconducters  are  prime  candidates  for 
optical  semiconductors  operating  in  the  green,  blue, 
jmd  ultraviolet  ponions  of  the  spectrum,  but  are 
limited  by  the  difficulty  of  doping  them  both  n-  and 
p-cy-pe — which  is  required  for  fabricating  light- 
emitting  diodes.  The  p-type  doping  is  especiaily 
problematic  in  H-VTs.  k  notable  e.xception  is  ZnTe. 
which  is  rather  easily  doped  p-cj'pe, !  ut  not  n-type — 
in  contrast  to  ZnSe  and  most  ocher  H-VI  semicon¬ 
ductors  which  are  r.-c>pe  dopable  but  not  p-typie.  Tr.us 
a  major  mystery  has  been  why  ZnTe's  doping  prop- 
erties  are  so  differen:  from  ZnSe’s  and  those  of  other 
n-VI’s. 

We  proposed  the  following  e-vplanation  of  ZnTe’s 
unusual  procihity  toward  p-type  doping:**^*  All  but 
the  most  carefully  grown  ZnSe  and  ZnTe  have  de¬ 
fects  which  produce  p-like  deep  levels  at  energies 
slight  above  the  valence  band  maximum  of  ZnSe, 
but  below  the  valence  band  edge  of  ZnTe.  As  a  re¬ 
sult.  the  deep  level  traps  holes  in  ZnSe  (making  the 
material  semi-insulating;,  but  provides  free  holes  in 
ZnTe  (making  ZnTe  p-t>-pe).  Namely,  if  random 
ZnSei.jTe,  alloys  could  ’oe  growr..  these  deep  levels 
would  cause  the  defect  to  undergo  a  deep  to  shallow 
transition  with  increasing  x. 


We  speculate  that  the  defect  responsible  for  the 
relevant  deep  level  is  antlsite  Zn,  namely  Zns«  and 
Zut,,  or.  in  the  case  of  Li  doping,  an  antisite  dopant 
Lis,  or  LiT.-*^  analysis,  however,  does  not  de¬ 
pend  on  either  of  these  ider.tincaiicns.  and  the  pic¬ 
ture  we  propose  depends  only  on  there  being  a  de¬ 
fect  with  a  dee?  level  slightly  above  the  valence  band 
maximum  capable  of  trapping  holes  in  ZnSe;  the 
corresponding  deep  level  in  ZnTe  must  lie  below  the 
valence  band  ma.viinum  and  donate  holes  to  the  va¬ 
lence  band.  (Fig.  7.) 

Our  goa*  is  to  construct  a  superlattice  with  a  band 
gap  naar  that  of  ZnSe,  but  with  a  superlattice  va¬ 
lence  band  maximum  chat  lies  above  the  ZnSe  hole 
trap  and  hence  will  convert  it  into  .1  shallow  accep¬ 
tor.  One  way  to  do  that  is  to  embed  thin  layers  of 
GaAs  in  ZnSe.  The  valence  band  edge  of  GaAs  is  at 
higher  energy  than  that  of  ZnSe.  while  the  conduc¬ 
tion  band  minima  of  the  two  .mat.^rlals  are  at  al¬ 
most  the  same  energy.'*  Therefore  by  choosing  the 
thickness  of  the  GaAs  we  can  "tune"  the  superlat- 
nce's  valence  band  edge  so  that  it  lies  above  the  ZnSe 
hole  trap,  without  greatly  altering  the  energy  of  the 
conduction  band  edge.  (See  Fig.  S.)  The  resulting 
predicted  (low  temperature)  band  gaps  are  given  in 
Fig.  9.  Band  gaps  even  further  toward  the  blue  are 
possible  with  i.GaA3)i.,{ZnSei,/Zn3e  superlat¬ 
tices.'* 


V.  SLISDLARY 

In  summary  we  have  .«hov.T.  that  in  tj'pical  su¬ 
perlattices  the  deep  level  energies  do  not  change 
much:  iNmically  of  order  O.i  eV  for  an  impurity  at 
an  ir.ten’ace.  The  j-iike  levels  shift  and  the  ,o-like 
levels  split  with  th-,j  p-orbita.s  directed  toward  one 
material  having  absolute  energies  close  to  the  bulk 
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3.0 


0.0 


-Lj 

. :..i 

ZnSe; 

2:8  eV  GaAs: 

1.3  eV  SL:2 

.3  eV 

/ 

Deep  Level 


Fig.  8  —  Schematic  .-ilustration  of  how  thin-layers  of  Ga.As  in  a 
GaAs/ZnSe  superiattice  can  heio  cover  up  a  deep  level,  wnich  is 
a  bole  crap  in  ZnSe,  causing  it  to  become  a  shallow  acceotor.  The 
boles  bubble  up  w  the  valence  ba.nd  maximum  of  the  superiat¬ 
tice.  after  Ref.  18. 
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Fig.  9  —  Predicted  banc  gacs  of  .V.,,..  x  10  (OOli  Ga-^-s/Zr.Ss 
superiattices.  snownng  how  the  band  can  can  bo  in  tne  green  :tr 
soike  suoeriattices  tsinail  .V, ’.tter  ?.ef.  1?.  .V.^.,  is  t.he  num- 
-er  ot  GaAs  ri-layers  in  a  perioc  of  :.~.e  suoeriaiticc. 


defect  energy  for  that  material  (e.g.  '.laAs  or 
.\do.';Gao.3As).  The  interiace  splitting  and  shifts  die 
out  rather  rapidly  with  distance  from  the  interiace, 
virtually  disappearing  within  three  atomic  bi-lay¬ 
ers  unless  the  level  is  p-like  with  nearby  p-iike  host 
states,  in  which  case  its  splitting  reflects  the  split¬ 
ting  of  the  host  superiattices  spectral  density. 

In  contrast  to  the  deep  levels,  which  are  rela¬ 
tively  unaffected  by  the  superiattice  order,  the  band 
edges  of  small-period  superiattices  exhibit  the  ef¬ 
fects  of  quantum  confinement,  and  so  can  lie  any¬ 
where  between  the  lower  and  upper  conduction  band 
energies  of  the  nvo  bulk  materials  which  constitute 
the  superiattice— with  the  e.xact  superiattice  band 
edge  energy  depending  on  the  ’.ayer  thicknesses,  es¬ 
pecially  the  small-gap  layer  (c.g.  GeAs).  Therefore, 
the  interplay  between  the  r.onnnement  effects  on  the 
band  edges  and  the  insensitivity  of  the  deep  levels 
to  the  superiattice  ordering  can  lead  to  shallow-deep 
transitions  which  change  the  character  of  a  dopant 
from,  say,  an  n-cype  donor  to  a  semi-insulating  deep 


trap. 

The  concept  of  shallow-deep  transitions  for  p-like 
deep  hole  traps  lying  slightly  above  the  valence  band 
ma-timum  of  ZnSe  has  been  proposed  as  a  possible 
eorplanacion  of  why  ZnTe  is  relatively  easily  doped 
p-cype,  in  contrast  to  other  II-VI  semiconductors. 
Finally  this  concept  has  been  used  to  design  GaAs/ 
ZnSe  small-period  superlanices  which  have  almost 
the  band  gap  of  ZnSe  and  almost  the  doping  prop¬ 
erties  of  GaAs.‘®' 
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Solf-Consistent  Antiferromagnetic  Ground  State  for  La2Cu04  and  CuO  Via  Energy  Band 
Theory 
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We  have  used  the  pseudofunclion  (PSF)  method  to  compute  self-consistent  spin-  polarized  energy 
bands  for  La2Cu04  and  CuO.  The  ground  state  is  found  to  be  semiconducting  and  anti-ferromagetic 
(AF)  for  both  La2Cu04  and  CuO.  For  La2Cu04  the  Cu  moment  is  0.35pB-  For  CuO,  the  moment  is  0.68 

PB  on  the  Cu  site  and  0.19  ps  on  the  0  site.  The  moment.*,  are  in  agreement  with  neutron  diffraction 
data. 


The  pseudofunction  (PSF)  method(1]  with  a 
local  density  potential[2]  has  been  used  to  compute 
spin-polarized  energy  bands  for  La2Cu04  (214) 
and  CuO  which  give  properties  in  good  agreement 
with  experiment.  The  PSF  nrtethod  uses  a  local 
orbital  basis  set  employing  the  full  potential 
including  the  core  States  and  non-spherical 
corrections  throughout  the  unit  cell.  Nine  s,  p,  d 
basis  functions  were  used  on  the  metal  atoms  with 
s.  p  sets  on  the  O  atoms.  The  basis  functions  are 
continually  changed  during  iteration  to  self- 
consistency  so  as  to  optimize  the  description  of  the 
charge  density.  The  pseudoluntions  and  the  non- 
spherical  part  of  the  potential  are  expanded  in 
plane  waves. 


In  Figure  1 ,  the  spin  polarized  bands  near  E( 
3fe  plotted  for  214.  The  band  gap  between  filled 
3hd  empty  states  varies  from  approximately  2  eV  at 
^to  0.06  eV  on  the  hexagonal  face.  The  bands  are 
ffllatlvely  flat  in  the  Z  direction  because  of  the 
'^ered  structure  and  are  not  shown.  For  the  spin- 
■  polarized  bands,  the  gap  varies  from  0.35  eV  at  X 
;.  and  M  to  0.06  along  the  M  to  N  and  N  to  S 
;®f8ctlons.  The  necessity  that  bands  be  narrow  on 
*^1,^18  hexagonal  face  thus  becomes  obvious  since 
•pin  effect  which  gives  the  semiconducting  gap 
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varies  from  0.35  eV  to  0.06  eV.  This  band  is  0.17 
eV  wide  for  the  PSF  method(3]  and  approximately 
0.5  eV  for  other  methods[4]. 

The  spin  polarized  energy  bands  for  CuO 
give  a  larger  gap  of  0.2  eV.  The  recent  crystal 
stmcture  determination  with  neutrons  was  used[3]. 
There  are  16  atoms  in  a  monoclinic  cell.  The  non- 
spherical  potential  was  expanded  in  18513  plane 
waves  in  order  to  allow  very  accurate  description  of 
the  potential. 


Figure  1)  Spin  polarized  energy  bands  near  the 
Fermi  energy  for  La2Cu04 
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In  Table  I.  the  moments  on  the  Cu  and  O 
sites  are  shown  with  comparison  to  recent  neutron 
diffraction  data.  The  agreement  is  very  close 
indicating  that  the  local  density  energy  bands  are 
giving  an  accurate  description  of  the  ground  state. 

Table  I.  Comparison  of  calculated  magnetic 
moments  with  experiment. 

LaaCuO^ 

Calculated  Experiment 

Cu  0.35pB  0.35pB 

0  0.0  small 

CuO 

Cu  0.68(xb  0.65|iB 

0  0.1  9hb  0.1 4pB 

The  importance  of  the  description  of  the 
potential  is  illustrated  in  Figure  2  where  only  the 
band  at  the  Fermi  energy  is  shown.  These 
calcualtions  are  non  spin  polarized  as  we  wish  to 
llllustrate  how  the  band  width  Is  sensitive  to 
description.  The  band  at  E(  is  very  broad  in  the  top 
panel  because  the  potential  was  limited  to  only 
4913  plane  waves.  The  band  in  the  bottom  panel 
is  very  narrow  because  15,625  plane  waves  were 
used.  The  wave  functions  and  their  plane  wave 
description  were  identical  for  both  calculations. 
This  result  apppears  to  be  a  convincing  means  of 
relating  band  widths  and  accuracy  of  description. 

Rnally.  the  CuO  bands  were  also  calculated 
with  spin  polarized  augmented  spherical 
wave(ASW)  method.  Identical  postions,  unit  cells 
and  k  point  sampling  were  used  in  the  ASW  and 
PSF  calculations.  The  ASW  method  found  the 
ground  state  to  be  metallic  with  no  magnetic 
moment  on  the  Cu  or  0  sites  contrary  to  the  results 
obtained  with  the  PSF  method.  Thus,  the  ASW  and 
PSF  methods  give  different  solutions. 

In  summary,  the  itinerant  energy  band  model 
yields  a  band  structure  which  has  a  moment  of 
approximately  the  correct  magnitude  on  the  Cu  site 
and  O  sites  for  both  214  and  CuO. 
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I  Relaxed-lattice  model  of  isolated  and  paired  isoelectronic  traps  in  GaP 

r. 

Jun  Shcn,  Shan^;  Yuan  Ren,  and  John  D.  Dow 
Dcponmrni  0/  Physics,  Unwersin'  of  Notre  Dome,  Notre  Dome,  Jndiono  *6556 
fRcceivcd  25  Ociobcr  19S9) 

A  simple  model  of  isolated  and  paired  isoelccironie  tubstiiutionaJ  traps  in  GaP  is  presented  that 
successfully  explains  the  anomalous  monoioruc-in-n  ordering  of  ihc  nth  nearest-neighbor  (NN),- 
pair  levels,  with  pair  energies  that  vary-  as  the  inverse  cube  of  the  N-N  separation.  Tne  model  is  a 
multiband  deep-level  iheorj’  that  includes  the  effects  of  lattice  relaxation  around  the  impurities.  It 
is  predicted  that  (Bi.Bi),  pairs  will  not  exhibit  the  monoionic-in-n  ordering  of  KN,  pairs,  that 
(Bi.Bil,  and  (Bi.Bi)]  will  be  resonant  with  the  valence  band,  and  that  hydrostatic  pressure  can  drive 
the  Bi  trap  into  the  valence  band. 


1.  INTRODUCTION 

Although  luminescence  associated  with  isoelectronic 
P-subsiiiuiional  defects  such  as  Bi  and  N  in  GaP  has 
been  studied  for  over  a  quarter  of  a  c:niurj‘,'“^  there  'till 
is  no  theory  which  can  explain  all  of  the  following  facts: 
(i)  isolated  produces  a  deep  holt  trap  =40  meV  abov; 
the  valence-band  maximum,*-^  (ii)  isolated  Nj,  produces 
an  electron  trap  1 1  meV  below  the  conduction-band 
minimum,*  (iii)  luminescence  associated  with  (Bi.Bi),. 
pairs  is  not  observed,  and  (iv)  discrete  luminescence  lines 
!  associated  with  (K,N)„  pairs  are  observed  with  the  energy 
of  the  nth  nearesi-ne'ghbor  pair  being  simply  related  to 
the  separation  P.,,  between  the  two  N  atoms  in  the  pair: 

£(NN%)=£(N)-/3P"-’  for  n>2  . 

Here  NN,,  refers  to  the  nth  nearesi-ncighbor  pair  (with 
both  K  impurities  on  anion  sites),  £(N)  is  the  energy  of 
isolated  K,  and  /?  is  a  constant.  Indeed,  existing  theories 
have  been  incapable  of  satisfactorily  describing  the  ob¬ 
served  ordering  of  the  NK,,-pair  levels,  and  do  no;  even 
predict  the  observed  monotonia  variation  with  the  sepa¬ 
ration  (or  the  neighbor  number  n }. 

In  this  paper  we  shall  assume  the  accepted  viewpoint 
that  the  isoelectronic  trap  ievcis  are  deep  levels*”'  due  :o 
the  central-cell  potentials  of  N  and  Bi,  and  we  shall 
demonstrate  that  (i)  a  strain-free  modal  cannot  explain 
the  observed  monoionic-in-n  behavior  of  the  KK,-piir 
lines,  namely,  if  the  impurities  and  host  atoms  are  as¬ 
sumed  to  occupy  silts  or.  the  undisioried  lattice  i.e.,  if 
the  efiifcts  o.'"  strain  due  to  the  size  mismatch  of  impurity 
and  host  are  neglected),  then  the  energies  of  the  N.N,- 
pair  lines  are  not  monotonic  functions  of  n  or  and  (ii) 
if  strain  is  incorporated  in  a  simple  model,  then  all  of  the 
above  unexplained  facts  can  be  easily  understood. 

The  successful  theories  of  the  prototypical  i.soelcciron- 
ic  trap,  isolated  substitutional  N  in  GaP,  have  all  as¬ 
sumed  that  virtually  the  entirt;  defect  potential  is  ic  cal- 
iaed  in  the  central  cell;  the  theories  differed,  however,  in 


the  number  of  host  bands  significantly  mi.\ed  into  the 
chemical  bonds  of  this  defect.  The  earlier  theories,  no;- 
ably  by  Faulkner'  and  by  Hsu  ci  cl.,”  employed  only  ore 
band  until  Swans  et  demorstraied  that  the  natural 
defect  potential  had  a  strength  so  large  (=•“  eV)  that  2 
muhtband  model  is  necessary  to  renormalize  the  defect 
potential.  .1;  is  now  generally  accepted  the;  the  minimum 
number  of  bands  required  to  produce  the  sp'  chemical 
bonding  is  eight:  The  K  trap  levc;  is  a  ‘‘deep”  levv)  origi¬ 
nating  f.'cm  the  central-cel;  potential,  as  desc^ribed  by  the 
Green’s  function  theory  ofKjalmarson  c;  0!.' 

n.  INADEQUAtTi'  Or  STRAIS’-?Rr.£  THEORIES 

Effons  to  extend  the  isc’.t'.eci-)'*  theory  to  NK',.  patrs" 
have  no;  p:'oduccd  the  observed  monoionic  ccpeitdctict 
on  r.  (or  a,,  );  set  Fig.  1.  ;>.;s  failure  might  appear, 

21  first  glance,  to  be  a  proper;)  of  the  specific  iheoreiical 
models,  closer  examination  of  the  theories  reveals  that  it 
is  a  general  propert)'  of  mcceis  which  emit  the  effects  of 
lattice  strain  and  which  i:isiead  at-tmpi  to  obtain  the 
NN,-pair  energies  from  titctronic  coupling  alone,  while 
constraining  ire  N  impuriii.;  a.-.d  the  host  atoms  to  sitta 
of  a  perfect  zinc-biend;  iztiiet.  7c  see  this,  consider  the 
Schrodinger  equation  for  an  isolcied-K  .mpurity  at  site 
A  (  A  =  0or  A  =  RJ: 


Here  is 
impurity  2; 
an  N’K,  pa 


the  deep-level  wav*  function  of  an  isolaier-.’C 
site  A,  and  is  the  defect  potential.  For 
:;,  with  one  N‘  impur.ty  .•■.i  the  origi:.  and  the 


other  at  R.,,  the  Schrodinger  ec  jaiicr.  is 


I -  k’o  -h  » lO  =  £ ;  1  N,  -.0  , 
with  approximate  bonding  and  amibonding  o:bi:als 
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FIG.  1.  Binding  enirgics  of  NN,  pairs  observed  by  Tnomas  a'ld  Hopneld  (Ref.  2),  and  calculated  with  the  present  model,  by 
raulkner  (Ref.  31,  by  B.-and  and  Jaros  (Ref.  161,  by  Li  f?  cl.  (LMX)  (Ref.  15),  and  by  Gil  tt  cl.  (Ref.  17).  The  zero  of  energy  is  the 
free*exciton  luminescence  line.  Note  the  ordering  of  levels  with  n.  The  NN|  binding  energy  for  the  present  theory  is  an  upper  bound 
estimated  by  assuming  d‘  =  d(,\  see  text. 


The  corresponding  energy  is 

£(NKj=£(N)-r(:r<do|Pol<iR.)-J- 

X(l±<(iold.R^>r' 

or 

£  ( KK  „ )  5=  £  ( N  )±  <io(  0  ( 0  )uo 

for  large  (This  approximate  result  agrres  weli  with 
more  exact  calculations.)  Here  the  overlap  integral  is 
^<iol^R, );  '*£  have  assumed  that  £,  is  moderately  large, 
and  we  have  taken  the  defect  potential  to  have  .s'rcngth 
Uq  and  to  be  localized  within  the  impurity’s  cell.  Kcnce 
the  NN-pair  energy  £{N.K*^ )  has  the  same  dependence  on 
Rf,  as  the  isolated-K  dccp-levcl  wave  function  centered  at 
Since  this  wave  function  is  'Known  to  oscil¬ 
late  '  as  a  function  of  £,  and  is  not  monoionic,  we  con¬ 
clude  mat  tne  observec  .N’K^-p.’.!:  energies,  which  are 
monotonic,  are  no:  determined  prccomir.ar.tiy  by  the 
electronic  coupling  between  the  two  impurities. 

III.  ROLEOFSTRAI.N 

Since  the  electronic  coupling  between  the  two  K  im¬ 
purities  in  the  KN’„  pairs  does  not  produce  the  observed 


monoionic-in-n  ordering  of  ;he  pair  levels,  there  must  be 
some  other  effect  which  is  larger  in  magnitude  than  the 
electronic  coupling  and  which  produces  monotonic  level 
ordering.  We  shall  show  that  elastic  strain  due  to  defor¬ 
mation  of  the  host  lattice  in  the  vicinity  of  the  defect  is 
responsible  for  the  observed  ordering  and  also  produces 
the  correct  magnitudes  for  the  energies  of  the  NN„  pairs. 

We  arc  not  the  first  to  suggest  strain  as  the  mechanism 
primarily  responsible  for  the  KN,-pair  energy  levels  in 
Gap,  Twenty  years  ago  Alien'*  proposed  that  strain,  and 
consequentially  deformation  pcteniials,  determined  the 
NK,-pair  binding  cne.-gies  in  GaP.  His  theory  fell  into 
disfavor,  however,  because  his  N-rciated  levels  were  rig¬ 
idly  attached  to  the  conduction  band  edges  and  w'ould 
have  preserved  their  energies,  with  respect  to  the  conduc¬ 
tion  minimum  in  GaASi_,P^  as  ihe  alloy  composition  x 
varied — contrary  to  the  observations.’  The  important 
demonstration  of  Allen  is  that  strain  products  an  effect  of 
the  correct  order  of  mesnitude  to  explain  the  NN„-pair 
data.  This  raises  the  qucsiior.  of  how  il.-  sirair.-frtf 
theories  were  able  to  obtain  KN’,-pair  energies  of  the 
correct  order  of  magnituc.  aloe;:  incorrectly  orn.ercd:  If 
both  the  cltctrcnic  coupling  and  strain  effects  aie  of  the 
required  order  of  magnitude  to  explain  the  data,  then  a 
correct  thto.">  must  account  for  the  cotripeiiiion  octwetn 
these  enects.  Hov'ever,  in  the  two  muliiband  theories, 
namely  ihc^e  of  1  i  ei  a!.'^  and  Jaros  and  Brand,"'  the 
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coupling  between  N  atoms  or  the  defect  potential 
strength  is  an  adjustable  parameter — and  this  parameter 
leads  to  the  large  apparent  electronic  coupling  with  in¬ 
correct  ordering.  The  theories  of  Faulkner^  and  Gil 
'  et  al.^^  are  missing  the  important  effects  of  the  valence 
band’  on  the  impurity  levels.  Thus  the  previous  theories 
.  that  attempted  to  explain  the  NN^-pair  energies  in  terms 
'  of  a  strain-free  electronic-coupling  model  considerably 
overestimated  the  effect. 


IV.  ISOLATED  IMPURITIES 

We  incorporate  the  effects  of  lattice  deformation  and 
strain  around  an  isolated  impurity  using  the  widely  ac¬ 
cepted  sp^s*  tight-binding  model  of  Hjalmarson  et  al.^ 
This  model  successfully  described  the  isolated-N  trap  in 
GaAS|_,P,  and  many  other  deep  levels  in  a  wide  variety 
of  Semiconductors.  Most  of  the  published  calculations 
based  on  this  model  did  not  incorporate  the  effects  of 
strain  or  lattice  relaxation,  however,  although  lattice  dis¬ 
tortions  are  easily  incorporated  into  the  nearcsi-neighbor 
tight-binding  model  Hamiltonian,  because  the  various 
matrix  elemenii  T  between  orbitals  centered  on  adjacent 
sites  obey  Harrison's  scaling  rule  and  so  arc  approxi.-nate- 
ly  inversely  propo.nional  to  the  square  of  the  bend  length 
d:  T  —  To(dfy/d)',  Tnerefore,  the  defect-potential  matrix 
Fof,  say  a  Kp  defect  in  GaP,  is  localized  to  the  impurity 
site  and  its  four  neighbors  and  to  the  basis  orbitals  s,  p^, 
Pf,  and  Pj  (Ref.  18)  centered  on  those  sites.  The 
impurity-site  diagonal  elements  are  the  same  as  given  by 


Hjalmarson  et  al.  for  the  strain-free  theory” 
F,=^.[u.',(N)-m,(P))  , 
and 

Here  the  lu's  are  the  atomic  energies  in  the  solid,  as  given 
by  Vogl  et 

The  wave  function  of  the  N  deep  level  is  dominantly 
hostlike  and  Ga-dangling-bond-like  in  character.  Hence, 
the  effects  of  host  d  states  are  well  simulated  by  the  sp's" 
mooei,  which  is  known  to  treat  the  host  electronic  struc¬ 
ture  well,  and  the  relativistic  effects  (which  arc  important 
for  Bi  levels)  can  be  neglected.  The  nonzero  nearest- 
neighbor  off-diagonal  elements  of  the  defect  potential  V 
arc  obtained  from  the  host-crystal  matrix  elements, 
scaled  for  the  altered  (inwardly  relaxed)  bond  length  of 
the  impurity  with  its  four  neighbors: 

^''cffd..fon.l  =  i5=ro[(tfo/t/)=-l). 

where  d  is  the  relaxed  bond  length  and  has  a  value  be¬ 
tween  df^i  =2.36  A)  and  the  sum  of  the  covalent  radii  of 
•K  and  Ga,  2.01  A.  Details  of  the  caiculational  procedure 
for  obtaining  the  N  deep-level  energy  are  given  in  Appen¬ 
dix  A.  * 

I 

V.  PHYSICS  OF  THE  RELAXATION  EFFECT 

The  qualitative  physics  governing  the  strain-relaxation 
effect  is  displayed  in  Fig.  2,  where  we  use  a  defect- 
molecule  model.  When  a  Ga  atom  (energy  Cq,)  and  a  P 
atom  (ene.'gy  tp)  are  brought  together,  a  bonding  and  an 


Defect  Molecule  Model  with  Lattice  Relaxation 
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r  -“10.  2,  Schematic  illusir.tiion  of  the  effects  of  the  centra'-ccll  defect  poteatta!  and  iaitice  relaxaticr  on  the  deep  impuntt  levels  of 
ji  Np;  (b)  Bip  tn  GaP.  See  text  for  discu.ssion. 
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antibonding  state  will  be  formed.  The  bonding  state  is 
broadened  in  the  crystal  and  forms  the  valence  band, 
while  the  antibonding  state  produces  the  conduction 
band.  When  a  N  atom  replaces  a  P  atom  (£o»> 
its  central-cell  defect  potential  produces  a  bonding  hyper- 
deep  level  in  or  below  the  valence  band  and  an  aniibond- 
ing  deep  level  in  the  band  gap.  The  depth  of  the  defect 
potential  is  related  to  the  difference  in  atomic  energies  of 
P  and  N  (Ref.  7),  and  is  =  — 7  eV.  Because  the  nearest- 
neighbor  transfer-matrix  element  Tis  almost  the  same  for 
N — Ga  and  P — Ga  bonds,  it  is  the  larger  (by  =7  eV)  en¬ 
ergy  denominator  which  causes  the  defect’s 

bonding-antibonding  splitting  [of  order  T5(eo,  — Ep)~'  in 
the  extreme  tight-binding  limit]  to  be  smaller  than  the 
host’s — leading  to  a  deep  level  in  the  gap.  Because  N  has 
a  smaller  covalent  radius  than  P.  the  Ga  collapses  inward 
and  the  transfer-matrix  element  T  increases,  owing  to 
Harrison’s  universal  rule  T  —  T^dl/d"^  (Ref.  211.  The  an¬ 
tibonding  deep  level  is  pushed  up  due  to  such  a  lattice  re¬ 
laxation  by  an  amount  of  order  ro[(do/^)^“  1  )(Eo. 
“En)  '.  When  the  impurity  atom  is  larger  than  the  P 
host  atom  it  repisc's,  as  in  the  case  of  Bi,  the  bonding 
deep  level  is  puiie.i  up  in  energy  by  appro.vimately 
T^i[^dQ/d)*  —  \  )(E(3,  — EBj)“'  as  a  result  of  relaxation. 

We  assume  that  the  relaxed  bond  length  is  a  linear 
function  of  Ar,  tne  difie/ence  between  the  covalent  radii 
of  the  impurity  and  the  host:^^ 

d  —dQ-r}.hr  . 


The  parameter  a  is  necessarily  positive  and  less  than  uni¬ 
ty,  since  the  relaxation  around  a  vacancy  (the  smallest- 
radius  ’'impurity”)  does  not  annihilate  the  vacancy.  We 
determine  the  precise  value  of  a  by  requiring  the  theory 
for  isolated  Kp  to  produce  exactly  the  observed  energy 
lc''cli_^^  namely  1 1  meV  below  the  conduction-band 
edge."’^^  Note  that  without  the  lattice-relaxation  effect 
the  N  level  would  lie  270  meV  lower  in  energy.  (See  Fig. 
3.)  It  should  be  emphasiaed,  however,  that  a  is  a  phe¬ 
nomenological  parameter  whose  precise  value  should  not 
be  overinterpreted.  We  have  performed  calculations  for  a 
variety  of  values  of  X,  using  various  models  of  the  defect 
potential,'''  and  find  the  physics  insensitive  to  the  detailed 
parametric  choice.*^ 

At  this  point  the  theory  is  completely  determined  for 
isolated  impurities,  and  we  can  compute  the  energy  level 
of  Bip — a  defect  that,  without  the  lattice-relaxation 
effect,  would  give  a  T- -symmetric  deep-ievel  resonant 
with  the  valence  band.  But  Bi  is  larger  than  P  and  causes 
an  outward  relaxation  of  its  neighboring  Ga  atoms — and 
this  relaxation  pushes  the  deep  leve’  up  into  the  gap  (see 
rig.  3),  where  it  lies  gratifyingly  close  to  the  observed  Bi 
energy. 

According  to  the  theory,  neither  .-  'p  nor  Sbp  p.'oduces 
a  deep  level  in  the  fundamcmal  band  gap  of  GaP~in 
agreement  with  experimen: 

When  hydrostatic  pressure  is  applied  to  GaP,  the  3i 
level  should  move  down  into  the  valence  band*'  (see  ~ic 


In  Ga.^s,  we  predict  that  the  N  i£\el  descends 
funcamenial  banc  gap  with  increased  pressure — 


Np  In  GiP,  A,  levil 


Blp  In  Gap,  p-llke  Tj  level 


FIG.  3.  Predicted  results  of  lattice  relaxation  v$  the  parame¬ 
ter  A  on  (a)  the  s-like  A  |  level  of  Sp  and  (b)  the  p-)ike  Tj  level 
Bi,  in  OaP.  The  conduction-band  minimun.  is  denoted  CB.M 
and  the  valence-band  n..^ximum  is  V3.M.  The  shaded  area  cor¬ 
responds  to  the  energies  of  me  host  bands. 
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•  served**  (see  Fig.  5). 

This  same  theory,  applied  to  Bi  on  an  anion  site  in 
.  GaAs  or  InP,  produces  a  level  above  the  valence-band 
1  maximum  in  InP  that  descends  into  the  valence  band 
'  with  hydrostatic  pressure  and  a  level  in  the  valence  band 
;  of  GaAs  that  remains  there — in  agreement  with  the  ob- 
f  servations.^®  Clearly,  this  simple  model  convincingly  ac- 
t-  counts  for  the  principal  experimental  facts  concerning 
I  isolated  isoelectronic  traps  in  GaP  and  other  III-V  semi- 
f  conductors. 

a 

VI.  NN,  pairs 

The  model  Hamiltonian  for  a  pair  of  impurities,  such 
as  KN„,  can  be  constructed  from  the  isolated-impurity 
Hamiltonian  as  discussed  in  Appendix  B.  One  new 
feature  arises  in  the  case  of  pairs;  one  must  (in  principle) 
determine  the  relaxed  bond  length  d'  for  each  bond  using 
elasticity  theory. 

An  important  point  is  that  the  strain  field  cicatec  by 
one  N  impurity  at  the  origin  influences  the  nearest- 
neighbor  bond  lengths  of  the  second  impurity  at  R„.  One 
N,  being  smaller  than  P,  causes  the  second  N  of  an  NN„ 
pair  to  move  toward  it  by  an  amount  approximately 
given  by  elasticity  theorj*.  The  second  N's  neighboring 
Ga  atom  that  is  closest  to  the  first  N  also  moves  toward 
the  origin,  causing  a  net  elongation  of  the  second  K — Ga 
bond  length,  given  by  the  gradient  of  the  strain  field: 

d‘  =  d{i-rA/R^)  forn>2, 

where  we  have  =2j-6a,  and  a =r ±0. 7904  a\^° 
This  dependence  of  the  strained  bond  length  on  is 
reflected  in  the  NK,. -pair  energies.  In  the  case  of  KN) 
p?irs  this  expression  for  the  bond  length,  which  is  based 


Pressure  Dependence  of  Isolated  N 
and  NN  Pair  Levels  in  GaAs 


Pressure  (koar) 

rIG.  5.  Predicted  dspenaences  on  pressure  o:':hc  Ga.As  sand 
eases  P  and  A  (heavy  solid  lines!  and  'he  isotlecironic  trap  K, 
in  vosipanson  wjih  the  cats  of  Ref.  Some  K.S',  pairs  are 
s.iso  shown,  and  have  sim.lar  pn^sure  oepenoer.ces  lo  the 
iso!aied-K  level. 


on  linear-elasticity  theory  and  the  assumption  that  the 
strain  fields  of  the  two  interacting  K  atoms  do  not  inter¬ 
fere,  breaks  down.  In  this  case,  both  N  atoms  will  relax 
toward  their  common  Ga  neighbor,  assuming  a  bond 
length  less  than  d(,,  GaP's  bond  length.  Thus,  for  «  =  1, 
we  can  only  say  d'<dQ.  To  simplify  the  calculation  of 
the  pair  energies,  we  assume  that  a  single  stretched  bond 
length  d',  that  of  the  N — Ga  bond  closest  to  the  first  N 
atom,  characterizes  all  of  the  N — Ga  nearest-neighbor 
bonds  of  the  second  N.  (In  fact  there  is  often  a  slight  an¬ 
gular  dependence  to  these  bond  lengths.)  Then  w-c  obtain 
the  energies  of  the  KN,.  pairs  using  the  Green’s-function 
method,  as  outlined  in  Appendix  B. 

The  NN„-pair  energies  are  in  excellent  agreement  with 
the  data,  as  shown  in  Figs.  1  and  6.  The  relaxation  of  the 
strain  is  responsible  for  the  and  monotonic-in-n  or¬ 
dering  of  the  KK„-pair  levels  as  obsers-ed  first  by  Thomas 
and  Hopfield.^  Note  that  the  other  theories  attempted 
unsuccessfully  to  explain  the  observed  ordering  in  terms 
of  electronic  coupling  between  the  two  N  atoms  rather 
than  in  terms  of  strain. 

It  is  noteworthy  that  the  theory  does  not  predict  that 
(Bi.Bi)„  pairs  order  monotonically  with  n  or  R„.  In  this 
case  the  isolaied-Bi  states  are  p-li):e  or  Tj  symmetric  and 
near  the  valence-band  maximum.  We  find  that  isolated 
Bi  lies  in  the ‘.gap,  as  observed,  and  that  some  (Bi,Bi), 
pairs,  for  n  >3,  lie  in  the  gap  as  well.  However,  (3i,Bi)) 
and  (Bi,Bi)j  lie  in  the  valence  band,  due  to  the  strain  asso¬ 
ciated  with  its  larger  radius  than  the  radius  of  P. 

No  (Bi,Bi)„  pairs  have  been  repontd  in  GaP,  to  our 
knowledge.  It  stands  to  reason  that  for  very  large  n  these 
pairs  will  lie  close  to  the  isolaied-Bi  level — and  so  m.’y 
not  be  resolved.  For  smaller  n  the  formation  of  enrugh 


Binding  Energies  for  NN^,  Pairs  in  GaP 
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FIG.  6.  Energies  of  KN',-pai:  luminescenrr  lints  in  GaP 
(with  respect  to  the  isolated  N‘r  ium.nesctnce)  vs  pair  separation 
divided  by  the  iaitice  consia-.:  of  Ga.*',  a.  Open  circles  .'.re 
the  pre.seni  theory  for  n  >  2:  data  are  denoted  b;.  soiic  triangle... 
Note  the  R depenrence  causer  oy  strain.  T'te  upper  bouro 
on  the  N’Nj-j.air  energy  cotaincd  by  assuttti.t:  c”  =  C{,  is 
displayed  as  a  bar. 
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pairs  10  observe  may  be  limi’.ed  by  the  Bi  solubility.  The 
interesting  feature  of  the  present  theory  is  that  it  explains 
why  (Bi.Bi),  and  (Bi.Bi),  should  not  be  observed,  even  if 
these  pairs  form  in  sufficient  concentration,  namely  they 
lie  outside  the  fundamental  band  gap.  To  a  good  approx¬ 
imation.  the  (Bi,Bi}„  levels  have  the  same  dependence  on 
hydrostatic  pressure  as  isolated  Bi,  according  to  the 
theory.  It  would  be  interesting  if  these  missing  levels  pre¬ 
dicted  by  the  theory,  and  their  pressure  dependences, 
were  measured. 
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APPENDIX  A:  ISOLATED  IMPURm' 

We  consider  an  impurity  (c.g.,  Bip)  in  GaP  at  the  ori¬ 
gin  with  nearest  neighbors  at  the  sites 

v,  =  (1,1, 1)0/4  , 

v,  =  (l,  — 1,  —  l)o/4  , 

v,  =  (- 1,1, -1)0/4  , 

and 

v,  =  (- 1,-1, 1)0/4  , 

where  a  is  the  lattice  constant.  On  the  impurity  (anion  or 
_ ^ _ ) 


"a")  site  we  have  five  orbitals  lO.o.j  ),  |0.c,;j,  ),  (O.o,;?,. }, 
!0.o,pj),and  lO.c.s*).  Since  the  defect-potential  matrix 
docs  not  involve  the  excited  state  s',  we  ignore  that  or¬ 
bital.  The  s  orbital  is  a  basis  for  the  A  j  or  s-like  irreduc¬ 
ible  representation  of  the  tetrahedral  group,  and  the  three 
p  orbitals  transform  according  to  the  T;  representation. 
Another  set  of  basis  orbitals  can  be  formed  from  the 
inward-directed  sp^  hybrids  centered  on  the  first  shell  of 
adjacent  (c.niion  or  “c")  sites  to  the  impurity.  They  are 

|l,A|,5)=(-i-|fi|  >  +  |fij>-f  |/i3)  +  |/i4>)/2  , 

|l,r},x  >=(-|/ii  )-l/i2)  +  |fi3)  +  l/i4>)/2  , 

and 

|l,r,.r>=(-|fi,>  +  |/i2>  +  |fi3>-ifi4>)/2  . 
where  the  hybrids  are 

|/i| )  =  (|v|.c,r )  -  |v,,c.x  > -  |v,,c,>' )  -  |v|.c.c )  )/2  , 

1^2 )  =  (  Kn.c,^  >  -  ivj,c,x  >  -r  |v;.c,y )  -r  Ivj.c.a  )  )/2  , 

Ihy )  =(iv,.c,r  >  -r  ivj,c.x  >  -  !v,,c,.v )  -r  |vj,c,:  >  )/2  , 

and 

1^4 ) ® {iv4,c,s )  -f  iv4,c,x  )  -f  IV4.C,;' ) -  |v4,c,: )  )/2  . 

If  the  defect  potential  matrix  V  associated  with  the  im¬ 
purity  is  confined  to  the  central  cell  and  its  coupling  is 
limited  to  first  neighbors  (the  latter  being  a  strain  effect), 
then  we  have  the  nonzero  matrix  elements  of  1' 

.:) 


lO.c.j)  ll,/<|,j>  lO.c.x  >  |l,r2.x)  !0.o,>’)  ll.r^.y)  lO.fl.:)  Il.T, 
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where  according  to  Harrison’s  rule,"  the  ini'ratomic  perturbation  matrix  elements  are 

Z),  =  [;^((do/r/)-i]  , 
and  we  have 

V,~  Tls,s)  —  ^T{sa,pc)  , 
and 


C',  =  r(sc,po)-i-r!x,x)-f  2T(x,y)  . 
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Here  and  d  arc  the  perfect  and  relaxed  bond  lencths,  respectively,  and  the  matrix  elements  are  the  same  as  in  Refs. 
20  and  7.  The  deep  impurity  levels  associated  with  defect  arc  obtained  by  solving  the  secular  equation 

det[l-G(£)n=0. 

where  G (£)=(£  is  the  host  Green's  function,  with  matrix  elements 


G(R,fc,i:R',fc’.r)=2<R,i,/kk){£-£(v.l:)]-'(v.kiR\6'./')  . 

v.k 

Here  the  sum  is  over  the  hosi  bands  v  and  ten  special  points’*  in  the  Brillouin  jonc.  Because  of  the  simple  form  of  the 
defect  matrix,  the  eigenvalue  equations  for  A  ,  and  T-,  levels  become 

;i-G(0,fl,r;0,o,j)I' )+[G(0,a,5;l,X,,s)’-G(0,a,f:0.<3.i)G{1..4,,r;l.^,.r)]£>;-2G(0,o,5:l,/I,,r)Z),=0  , 
;i-G(0.c.x:0,fl.x)RJ  +  [G(0,a,j::l,r„j:)--G(0,fl.x:0,u,x)G(l,r,.x;l,r,,A:))D;-2G(0,c,x:l,7'j,x-)D,=0  , 

threefold  degenerate). 


This  calculation  has  not  included  the  effects  of  bond- 
ength  changes  between  the  first-  and  second-nearest 
neighbors  for  two  reasons:  (i)  we  wanted  to  present  the 
.implest  possible  model  containing  the  essential  physics 
of  (K,N)  pairs,  and  (ii)  these  effects  are  very  small,  be¬ 
cause  they  affect  In  first-order  only  the  hybrid  orbitals 
-•entered  on  the  nearest-neighbor  sites  and  directed  cut- 
uard  from  the  impurity,  whereas  the  level  positions  arc 
determined  primarily  by  the  hybridization  of  the  impuri- 
y  orbitals  with  the  inioorJ-directed  hybrids. 

APPE.NDLV  B:  NK. -PAIRED  IMPURITIES 

For  paired  impurities  the  defect  potential  F,  is  a  direct 
um  matrix 


I - 

r,  =  I'(0)-f  K(R„)  , 

where  each  of  F(0)  and  F(R„)  is  an  SXg  matrix,  and 
the  bond  lengths  in  these  matrices  are  the  appropriate 
values  d'  for  the  pair  (as  opposed  to  d).  The  resulting 
secular  e.quation  is  a  16X16  matrix: 

I 

det(l-G(£)K)  =  0  , 

which  is 'solved  for  the  pair  energies  £.  The  ten  special 
points  were  extended  to  240  points  by  the  24  symme¬ 
try  operations'^  and  then  were  used  to  evaluate  the 
Green’s-function  matrix  elements. 
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tential,  we  use  the  notation  T  for  the  transfer-matrix  elements 
wheie  Vogl  c:  a!,  used  V. 
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’■The  covalent  radii  of  N.  P,  As,  Sb,  and  Bi  are  0.75,  1.06,  1.20, 
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tionality  of  impurity  bond-length  change  and  covalent  radius 
difference  has  been  demonstrated  by  J.  L.  Martins  and  A. 
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phenomenologies  differed  by  almost  a  constant  and  by  topi¬ 
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DEEP  LEVELS  IN  SUPERLATTICES 

JOHN  D.  DOW.  SHANG  YUAN  REN.  JUN  SHEN,  AND  MIN-HSIUNG  TSAI 
Depvtfticni  of  Phyiici,  liniveriity  of  Notre  Dune,  Notre  Dune,  Indinii*  4C530  U.S.A. 

ABSTRACT 

The  phyiie*  of  deep  leveli  in  temicondiieton  ii  reviewed,  with  einpliiisin  un  ll  e  Iki 
lliht  kll  jiibetituiiouM  impuriliei  produce  deep  levi-li  -  »o!iic  of  whieli  iiiiiy  not  lie  wi  hin 
the  fundiunenthl  biutd  gnp.  The  ehnrkrtrr  of  »  doputt  ehaiiGex  whei!  o;ie  of  tli'  deep  levels 
moves  into  or  out  of  the  fundtmentel  {up  in  response  to  &  perturbation  s-.irh  tun  pressure 
or  ehuige  of  host  composition.  For  exunple.  Si  on  a  G«  site  in  G»As  is  s  shuliow  donu-. 
but  become*  •  deep  trap  for  x>C.3  it.  .Af^Gai.^As.  Such  shallow.di'ep  tntu'itions  can 
be  induced  in  superlattir  ,*  by  ehansins  the  period-widths  and  ({uantuiii  rcuhi.eni  -ut. 
good  rule  of  thumb  for  doep  levels  in  tuperlattiee*  it  tha*.  the  en*rg.v  levels  with  ivspeet 
to  t'acuum  are  relatively  insensitive  (on  a  >0.1  eV  scale)  to  tuperlattic’'  pwriud-wItUiu, 
but  that  the  band  edge*  of  the  superlattiee*  ue  sensitive  to  el  uiges  ef  period.  Heuee 
the  deep  levol  positions  rtltUvt  to  (he  hand  edges  uc  sensitive  to  the  ptiiod- widths,  and 
shallow-deep  transitions  can  be  induced  by  b..nd-gap  engineering  the  superlattice  periods. 


DOPING  AND  DEEP  LEVELS  IN  BULK  SEMICONDUCTORS 

In  recent  years,  the  theory  of  doping  has  been  revised  from  the  old  cfT'  tive-mas* 
theory  |lj  which  explained  the  behavior  of  substitutional  P  i:i  Si  and  S"  <v.  an  A*  kite 
(Se:^,)  in  Ga.As  so  well.  In  the  effective-mass  model,  one  simply  aMumed  that  an  impuriiy 
such  as  Se  substituting  fb.  .As  would  have  a  singly-ionized  riable  grouiiil  stun-  Sr'^  were 
it  not  for  the  long-rangvd  Kreencd  Coulomb  potential  -e*/(r  which  binds  the  "extra” 
electron  of  Se  in  a  hydrogenic  shallow-donor  orbital  with  radius  aa:)f-r/m‘'<-*  mucii  larger 
than  a  typical  l.-i'tire  constant.  (Here  m*  is  the  effective  mass  of  the  riectiun  r  is  the 
dielectric  cunsthni  of  the  ho<:.  and  we  have,  for  kimpliciiy  of  prt-.«-ntulii' ;.  nv-  ;h«il.t-<:  tl 
anisotropy  of  m  .)  Hydrogenic  cffcciive-mnss  theory  h.is  a  significa;::  wt  :d:ii>'! in  ■■.••vcr: 
it  do'-s  ni'-.  srcdiel  ilia:  i.  |>.trticular  impurity  such  .xs  Sc  m  Ga.As  should  l>',  ,^  ^hallow  donor, 
but  rathe;  lii.ii  all  atoms  comi.-ig  from  Columns  to  the  right  of  .As  (Coiuiun-V) 

ill  the  Pv.’iodic  Tabic  arc  donors,  when  they  substitute  for  .As.  Tlic.eforc  effect ive-niass 
ihi-ory  offers  no  criicnon  for  dciermining  when  the  extra  electron  of  an  impurity  such  a.<  Sc 
in  C.i.A.'  >:.iiuld  occupy  u  deloc.iiixcd  "shallow"  l-.ydrogciuc  orbital  ru:!u  •’  tliaii  a  loc.ili.-<  'i 
"ducp”  oriiital.  bu:  rather  as.'.-.uiu-!.  tha:  a!!  impurity  Irvi-is  are  slndlnw 

Dcc.-ium;  tlic.'c  «;u.  iio  .accepted  theory  of  deep  iinpuriiy  li-vcl.->  in  n  niicoiuiuctorx  .scv-. 
cr:d  years  ago.  an  luiforiiinaic  semantic  problem  hius  ari'Cn  in  the  liicratur.-.  Sh.dlow 
enc-gy  Ic'cU  wc.'c  deiiiii-d  as  lesels  wiiam  arO.l  t-V  of  a  hiuid  edge  and  ’•iciir-  ih  ••u.'.lly 
loiuvahlc  -viiii  .sigiiirirtuil  probah.lily  at  rcMSm  temperature  Deep  levels  wi-;e  I'-vcl.  m  the 
/unce- (Oita/  kcr.d  gay  by  more  ll.i-.n  0.1  eV.  With  t!ic  adeem  of  the  'hcory  of  dirp  levels 
(2j.  tiii.-.  c.ncrgy  dcFitiitio:.  has  been  kuppla.iicd  by  one  dial  n-ii.-s  lU'irc  i---  i.bc  lo<':i!;r.ntion  of 
the  -:ale  :)iaii  ;t.s  energy  wi;h  respect  to  a  nearby  liaiul  edge  Shallow  l•■•.cl^  i-.ri.M-  from  tin- 
lung-.. meed  Coulomb  poieii'.i.a]  of  the  detect,  and  h.i-.-e  extended,  hydro  .envcfuiicti->ns 
d-  'CiiIk-.I  iiy  i-tr'  -iive-imcs  ilieory.  Di-ep  levels  originnle  from  the  ccir.i ell  piilemial  of 
:!»•  ii<  f<  c;  .As  .i  result,  all  !e-.cU  that  were  sii.i'lnv.  lij  siic  pii-Moiis  dcr.tr.- .o-i  ,iie  sh.illow 
lo.  riie  ne.'.  one.  '-.wpl  those  few  deep  levrlk  toy  ihe  new  (ii-fniition)  tii.il  .u  idrn:.!!!;.  .e 
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within  0.1  eV  of  »  bund  <d(l.  'An;^i(^^’'t^?wtrc  dcop  by  th«  old  definition  ere  »l*o  deep 
by  the  new  one,  but  now  there  m  ^viou^y  untntieipkted  deep  le%’el$  thet  mny  lie  in 
the  bendi  or  in  the  £kp  within  Oil  bV  of »  bm  .cd|e.  In  thif  pnper,  we  idw>a,vt  ndopt  the 
fullowins  n>o^m  definitiMt:  A  thidlow  donor  (eeceptor)  ir  n  defect  whore  only  enert;}' 
levels  in  the  fundementn!  bend  gnp  arc  ihallow  ef^eetiv^maa5  theory  level*.  A  deep  trap 
it  a  defect  that  hu  one  or  more  deep  leveU  in  the  gnp.  An  iioeleetronie  trap  it  a  deep 
trap  aatoeiated  with  an  itoelecironir  tubttitutional  impurity.  An  isoelectronir  rcH'imni 
tcattcring  center  it  an  itoeleeironic  imput'ty  a<hotc  deep  level*  lie  enereeticaliy  withiti  the 
host  band*  and  hence  resonantly  scatter  carriers. 

The  character  of  a  tuhsti  tut  tonal  impurity,  namely  whether  it  it  a  shallow  du'itir.  ii 
shallow  acceptor,  a  deep  trap,  an  isoeleetrunic  trap,  or  an  itoeleeironic  resonant  srn.icrins 
center  it  now  knoa'n  to  be  determined  no*  solely  by  the  position  of  an  impurity  in  the 
Periodic  Table  relative  to  the  host  atom  it  tcph.i  e*.  but  alto  by  the  strength  of  the  central. 
celt  defect  potential  and  the  energietof  the  “deep*  level*  pro,liieril  by  that  potential  relative 
to  the  band  edges  of  the  fundamental  gap,  Missing  from  the  efrcetivc.mass  moc.I  are  the 
local  bonds  that  arc  perturbed  by  tbs  d^ect  (and  their  bonding  and  anlibotiding  states), 
and  a  sjmptom  of  this  short-coming  is  the  fact  that  the  s-atomic  energies  of  Se  and  A* 
differ  by  M  rV,  indicating  that  the  central-cell  potential  which  model*  the  difference  must 
be  larged  yet.  Clearly  one  cannot  introduce  a  ar4  cV  attractive  potential  into  a  solid  atic 
have  its  ssfc  cfect  be  the  formation  of  shallow*  donor  state*  with  binding  energies  of  less 
than  10  me\';  something  must  also  happen  to  the  electronic  structure  on  tiiv  i  eV  settle. 
What  happens  is  the  formation  of  the  hyperdeep  and  deep  let-els. 

Pig.  {1}  (3.4)  illustrate*  the  case  of  K  sulwtituting  for  P  in  CaP,  but  (for  *iiii|>lir. 
ity)  ronsider*  only  a  sin(^  directed  orbital  of  each  atom.  The  Ga  and  P  atuma,  when 
brought  together  into  a  molecule,  form  boiiding  t  nd  antibonding  orbitals,  with  a  bondinr.- 
antibonding  splitting  invcr.eely  proportional  to  the  atomic  energy  dih'erenee  In 

the  solid,  these  molecular  levels  broade;:  into  band*.  Wliei,  N  substitutes  fur  P.  it*  energy 
»(<  is  as"  eV  lower  than  that  of  P;  hence  the  energy  denominator  «Ga’^N  '*  < 

than  <G)|-*p  and  tlie  impurity  bondinr-antibonding  splitting  is  corr«i>ondins!y  smaller, 
causing  the  aniil>uiiding  deep  level  lu  li.-  in  the  gap.  Tlic  .'C-lifce  iinpuriiy  h-vi-l !.-  tJic  mt 
ohiicrved  hyperdei’p  level  below  the  viilcncc  band  maxinmm.  whcrch.c  the  dec;- 1  -.cl  t:i:-.: 
lies  in  the  gap  ha.-  a  wavcfunction  tlini  is  Ga-dangling-bond-likc  (almost  inde)'.-:idcr.:  uf 
the  N  iinpiiriiy).  If  tl.c  conduction  hand  nerc  somewhat  broader  in  energy,  tlic  nevp  level 
would  be  covered  up  (as  happens  for  N  in  Ga.At). 

N.  being  iMielvctronic  to  P,  simply  produrr-  «  deep  level  in  GaP  that  iie*  in  the 
fiindiiii  '-ntid  l)and  gap.  slightly  below  the  '-..ndi:  tion  baml  edge.  Fo;  Se  in  Ga.n-  i  :ui  P 
in  Si.  ih.  ImiKiiiig.a  ihondiiig  physics  is  the  same,  bni  llier<-  is  an  extiii  cieciro;.  (tin- 
“diiiiur"  eleriruii)  that  will  occupy  ihi.«  deep  stale  if  it  is  within  the  gap.  oi  alter:. aiivcly. 
if  tile  deej)  level  lies  ab"--:  the  conduction  band  edge,  t}:*-  extra  electron  will  tali  to  th<- 
riiiidiiciiuii  band  eiigc.  produce  Sef  or  P~.  luid  be  bouinl  to  the  ionized  impurity  i..  ti.e 
hydrogeiiic  effective. mass  donor  state. 

For  Se  ill  G.a.A.s  (or  for  P  ii,  SI  ili-  ey;,erdee}>  or  iauidn-.e  levels  hav.  «.like  and  p-iil.e 
syiiiiiieiry.  arc  folly  occupied  by  electro.,-.  an<:,  i:‘-  vtil  below  the  valence  b.nu!  mnxitr.'iiii 
loo!  are  lo-  ihserted  and  norinaliy  miii!!ei<.:ii::c  Tn-  n:.;il>utiiii;:g  or  I'-vels  of  Se 
ill  i..i.As  are  also  s-iike  luid  p-lisc.  awi  he  e.-'-ec  tie  (.unc.i  lioa  in-.nd  n.imiimin.  Thes.- 
ilii-p  !i-v..|s  ill  parneolar  arc  missinc  (lom  ihe  erfec:!-.,  -tn.,  >  modi  1.  ami  ciii-i’.i-.'  •.■I’i/ 
rrfjiirt  tn  tin  noiLurtien  band  rail-  d>  iintmn  tin  dnfini)  charnctrr  o!  l!i'  nipuTiii/. 

Ill  i:.-effr.eii\e-ii!;u..>i model. !  ioiio:.  -•:lisliti::!ii;.ii]  ilop.i.-its oxyg* soh'iil .  si  leiinim. 
.iml  ii'lli.riiioi  -iiiiuK!  .-d!  he  shallon  doim:-  in  G.i.n-.  GaP.  Ga.A-j.jjl'.^  A'e:  ii'.;  ci':i 
h.is  bi-eii  nil-.  ; Mil  In  he  a  dee)>  iraji  III  G.  is  iliougii!  lo  la-  dee;,  m  Ga  a-'  w,l! 
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Fig.  {1).  Schematic  cncrgj-  levels  of  Ga,  P,  molecular  GaP,  and  solid  GaP,  togctlicr 
with  the  encrgj*  levels  of  a  N  anion-substitutional  impurity  “defect  molecule"  in  GaP,  after 
Refs.  [3]  and  (4). 
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Fig.  (2).  Energx'  level  diagratii  illustrating  the  •iifT-rcnce  bciwccii  shallow  iiujuirilies 
such  lu*  S  and  Sc  in  Ga.^s.  which  have  their  d'-ep  Icvrl.s  (''f  .Aj  and  To  symmetry  for 
substitutional  defects  m  zincblendc  hosts)  outside  the  funthunenlal  band  pap  and  deep 
impurities  such  as  0  in  GaAs  or  GaP,  which  have  at  Ic.xs;  one  tieep  level  in  the  gaj). 


fS).  In  contrut  to  oxygen,  3  tuid  Se  we  ihnllow  donors.  The  rewon  for  the  diffrrent 
doping  chweeter  in  GnAs  of  S  wid  Sc  on  the  one  hand  and  oxygen  on  the  other  is  that 
llie  sdike  (Aj'tymmelrie)  antibonding  deep  lc\*e)  of  oxj'gcn  lies  i«l*w  the  conduction  bund 
minimum,  but  the  corresponding  deep  levels  of  S  and  Se  lie  sieve  the  conduction  baud 
edge.  (See  Fig.  (3}.}  As  a  result,  the  ground  state  of  the  neutral  oxygen  impurity  U 
GaAs  has  the  deep  level  in  the  gap  oeeupied  by  one  electron.  This  level  lies  about  0.4  c\' 
|5]  below  the  conduction  band  edge,  and  so  is  not  thermally  ionizable  -  and  does  not 
contribute  to  the  semieonduetivity  of  GaAs.  MoreoN'er,  the  neutral  oxygen  impurity  can 
trap  another  electron  of  opposite  spin,  and  so  rather  than  donating  an  electron  and  mailing 
GaAs  tt’type,  it  can  remove  one  and  make  the  host  semidnsulating. 

In  contrast  to  oxygeti,  S  and  Se  are  shallow  donor  impurities  that  make  GaAs  i|. 
tyiie,  because  their  Ajsymmetrie  antibonding  deep  levels  lie  sieve  the  conduction  band 
minimum.  As  a  result,  their  e.rtra  electrons  we  autoionised,  fall  to  the  conduction  band 
edge,  and  we  then  bound  (at  wro  temperature)  in  the  hydrogenic  shallow  donor  levels 
pr^ueod  by  tlie  ionist-d  S'**  or  Sed*.  These  shallow  levels  lie  witldn  0.1  eV  of  the  conduction 
band  edge  and  so  can  be  thermally  ionised  rather  ewily  at  room  temperatures.  Since  the 
A}  deep  level  lies  well  above  the  Fermi  'tnergy  and  the  Coulombic  potential  of  S'**  or 
Se*^  it  screened  by  the  donor  electron,  each  of  these  impurities  is  incapable  of  binding  a 
second  electron  (6).  Thus  the  shallow  donor  impurities  S  and  Se  are  thermally  ionised  and 
contribute  to  the  mtype  semieonduetivity  of  GaAs.  In  contrast,  the  deep  impurity  oxygen 
traps  electrons  and  detracts  from  condtietiWty. 

A  central  notion  of  the  modem  theory  of  doping  is  that  ikt  ckersetcr  of  s  iopent 
u  dstermiaee'  kfi  the  rtUlivc  yesitisns  •/  its  deep  levels  end  (he  bond  edges  (and  .’'ermi 
energy). 

Fig.  (3)  shows  i.ow  the  doping  chwacter  of  an  impurity  such  as  Sc,^^  in  Ga.As  could  be 
changed  from  a  shallow  donor  to  a  deep  trap  to  a  falw  valence  shallow  uccepto: ..  provided 
one  could  perturb  the  band  edges  of  the  host  Ga.As  relative  to  titv  deep  level  s<  that  hist 
the  conduction  band  edge  and  then  the  valence  band  maximum  pa.«sed  through  the 
symmetric  decjj  level.  In  the  ca«cs  considered  in  Fig.  (3).-  the  bonding  s-like  and  p-like 
level.'!  we  fully  occupied  by  electrons,  lie  well  below  the  '-alencc  baud  maximum,  and  an- 
electrically  inactive  and  uninteresting.  Tlie  ami-boix^sig  p-likc  levels  lie  well  alxivc  the 
conduction  band  minimum,  and  are  unoccupied  and  also  uniiiiercsting.  Thus,  foi  Se,\5  it: 
Ga.As,  the  interesting  deep  level  is  the  s-likc  or  .A] -symmetric  deep  level  and  it  lies  above 
the  conduction  band  edge.  Dut  wr  could  imagine  applying  a  very  larg"  pressure,  winch 
would  cause  the  Ga.As  cunducii(>:i  hand  cdp.  to  jwiss  above  the  Sc  dir;'  level  -  changing 
the  doping  ch.irncirr  of  this  impurity  from  a  shalinw  donor  to  a  drep-irnp.  and  causing 
the  coiidiiclivity  to  (!ro]>  pri'C!)>ii<iusly.  One  could  uIms  imagine  ajiplying  it  uniaxia!  stress 
.so  strong  that  the  Sc  ilirp  level  would  p.ass  i:!io  the  valence  hand.  (The  m.'igniiudo  of  the 
stress  would  actually  exceed  the  breaking  i>oiiii  of  Ga.As.)  In  this  case  Sc.^j,  would  bcc»:nc 
a  sliallow  acceptor  l)ecause  the  liole  in  tin-  deep  level  would  bubble  up  to  the  valence  h-md 
iniiximinn  or  Fcrini  energy,  leaving  the  siahie  ion  of  Sc  to  l>e  Sc’’.  .A  liolc  could  tl!c:i  |••oit 
the  negative  ion  in  a  hydrogenic  cffective-tna-'S  state.  AVc  cal!  this  case  tlie  case  of  false 
valence  It  does  not  nee  :r  often,  liccaii.si  it  goiicrally  r-tpiircs  a  small  b.'ind-gau  host, 
whose  v.icjincy  luis  no  dirp  level  in  tlie  g.ap.  NeierlhelcfS  when  li.*'  co:i(hiioiis  for  fal.'c 
v.-ilciicc  are  met.  an  impurity  such  as  Sej^,.  from  C«>Imnn-Vl  can  appear  to  ..rt  as  if  it  "tie 
from  Cohnnn-lV  (from  the  viewpoint  of  soin*  >iie  familiar  only  with  cffcciiM'.iii.i.ss  lin-ory) 
In  gem  ral  ihe  false  laleiice  will  differ  from  liie  true  valence  by  two  nr  six.  (lepe  id.ng  on 
whether  llic  cics-p  ie\c!  'li.il  hs.s  crossed  lie-  fuiniaiiiciital  band  gap  is  an  s-siaic  nr  a  _.  /tii:e. 


Shallow  Daip 
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Fig.  {3}.  Energ}*  level  diagram  illustrating  how  the  relative  positions  of  a  deep  level 
and  the  band  edges  could  affect  the  character  of  a  dopant.  Band  edges  arc  long  solid  lines; 
shallow  levels  are  duhed,  and  the  deep  level  is  a  short  solid  line.  The  relevant  deep  level  is 
assumed  to  be  s-like  and  capable  of  holding  two  electrons  of  op|>osite  spin  This  impurity 
is  assumed  to  come  from  one  Column  to  the  right  in  the  Periodic  Table  of  the  atom  it 
replaces  (such  as  Sej^g  in  GaAs).  If  the  relevant  deep  '•  vel  is  in  the  gap  fb).  it  is  occupied 
by  an  electron  (solid  circle)  and  a  hole  (open  triangle),  and  the  neutral  impurity  is  a  deep 
trap  for  an  electron  or  a  hole.  If  the  level  is  in  the  conduction  band  fa),  the  electron  is 
autoionized.  falls  to  the  conduction  band  edge,  and  is  trapped  in  the  shallow  donor  hwel 
(dashed).  If  the  level  is  in  the  valence  band  (e),  the  hole  will  bubble  tip  into  the  shallow 
acceptor  level. 


F*fi.  (•«).  Tin-  baud  r.ipcs  al  the  X-poiii.  .and  the  r.poini  of  the  B.-illoinii  yone.  the 
N  di-cp  level,  the  O  d(a;i  I.-m  I.  and  the  S  .and  Sc  sliallow  levels  in  G.a.Asj.jjP,  '.crMi' alloy 
coiir.K, 5111011  .after  fief-  '.  .aid  [?' 


altering  the  character  of  a  dopant 

IntmtUng  phytic*  ocain  when  »  b*nd  «<)pe  {Mutoi  thropih  »  dcej>  1ml.  bcctuic 
lh»!  cliitnicter  dopnnt  citiut|e*.  This  effect  wtt  ffrti  cxpliciU,v  recognited  by  Wolford 
et  nl.  (5),  who  ttudied  wiiop-tiii  N1  0,  S,  *nd  Se  in  GnAtj.xPx*  impurity-hoti 
combinotiont  were  etpedelly  im  mtii*  r,  beceute  0.  S,  wtdf  e  «»e  «l  ColumnA'I  impurii  !<•» 
tnd  were  nelvoly  cximted  to  Imw  ila*  tnme  doping  ehurnctcr  m  iltnllow  donor*,  but  N  i» 
k  ColomnA'  itocleetronie  defect  and  thould  not  haw  produced  any  level  according  to  tlie 
then<ronventional  thinking.  ln*iead,  only  S  uid  Se  pr^ured  rhnllow  donor*.  N  prorhiccd 
an  “iaoelertronic  trap"  in  GaP  that  actually  lie*  flo»er  to  the  conductiun  band  edge  than 
the  thallow  donor  »tatct.  We  mm  know  that  thl*  itoelectrott'C  trap  ttate  it  in  fact  a  deep 
level  that  by  accident  lie*  only  /lightly  below  the  conduction  band  edge,  but  whe-  it  wu* 
firtt  ditcovt^in  Gap.  the  itoelcrironic  trap  concept  wa*  novel  and  no  one  eoiiMtlercd  the 
pottibility  that  the  trap  level  wa*  in  fact  a  deep  level  (with  ‘■deep”  defined  in  a  loealiiatiun 
«ente)  although  it  wa*  energetically  *hallower  than  the  (hallow  leveh'  Oxygen  in  GaP 
produced  a  level  energetically  deep  in  the  gap,  ml  eV  below  the  conductiun  band  edge. 

The  key  to  the  difference  between  S  and  Se  on  the  one  hand  and  oxygen  on  the 
other  i*  provided  by  the  data  for  N  in  the  alloy  GaAt^.xPx  l^l*  ^ 

0  energy*  level*  behave  (itvilarly,  bdng  roughly  linear  function*  of  alloy  compotit'on  x. 
while  the  S  and  Sv  level*  are  ratlier  obvioudy  “attached"  to  the  band  edge*  and  *0  vary 
quadratieally  tFig.  {4}).  Tin*  fact  indicated  that  the  phyeic*  reeponeible  for  the  N  and 
0  level*  ii  the  aame,  and,  riiice  N  i*  Ifoeleetronic  to  kt  and  P  and  hence  ha*  no  lone- 
ranged  defect  potential,  that  the  long-ranged  or  Coulombic  part*  of  the  defect  potentiui 
arc  not  retpontible  for  the  0  Ivvcl.  Thi*  wa*  an  important  fact  becauie  it  thowed  that 
the  0  level  wa*  not  a  (hallow*  level  that  *«mch<m*  appeared  energetically  lower  in  the  gap 
-  lower  due  to  *oni”  extra  phytic*  that  it  unintportam  in  mo»t  aliallow  donor  problem*. 
The  data  demon*: 't.tvd  that  the  0  deep  level  originated  inttcad  from  the  ccntrul-ccll 
potential  of  the  deu'ct.  Thi*  meant  that  liim*  were  two  teparate  limit*  to  the  immitity 
problen:;  ['.)  ordinary  (hallow- donor,  onc-band,  eff«*eiive-m***  theory  in  which  the  long* 

ranged  Coulomb  potential  of  the  impurity  created  hydrogenie  (tates  loraliied  in  k  •(pure 
near  a  hand  cxtrcnvim  and  delr>caii*rd  in  7*-sphce.  and  (2)  deep-level  theory  wbicli  relied 
upon  tin*  loealixcd  renirul-ccll  )>otciiti;tl.  crcateii  lionding  atid  aniibundiug  statvs.  u::t!  Iiad 
sp^-chtiracter  wavefuiictions  involving  at  least  '-ight  energy  bands  -  tlie  w.tvcfunction*  o: 

(urh  levels  arc  delocalixed  in  k  -spuee. 

A  direct  rontcquencc  of  tliis  way  of  thinking  about  impuriiic*  w.ts  the  notion  ilitti 
every  *-  and  p-lx>nded  impurity  produces  four  localized  bonding  and  four  uitiiboiidtng 
deep  levels  ju-sorinied  with  tlie  impurity^  jxTturbet!  Imi.ds.  IjiiifiiPj  the  am  It" ‘tiding 
levels  of  an  ehTironcgaiive  defect  arc  nciir  the  band  gap.  wliili*  tlie  building  le.••;*  are 
in  iir. below  the  valenre' btnid,  are  termer!  liyperdirp  levc!>  |2).  and  nurmidl;.  an  t.i’itlier 
electrically  active  nor  interesting.  Tlierefore.  ;!.••>•>•  must  be  liolii  an  s-like  and  a  llifx-fold 
degenerate  p-like  ticep  level  of  S  (or  Se)  in  G;.*'*  that  lies  aioei  •.lie  conuuriion  btu.ri  edg' 
and  varies  ajiproxiinatcly  linearly  with  x  in  G.i.Asj.^Pji.  -  as  the  X  und  0  (s-like)  I'-  !s  do. 

The  dtlti  for  N  in  Ga.Asj.jjPj;  indi  ate  il.tit  deep  icvci.  nt.mely  an  an:!ii(>:i.img  ('•• 
like,  in  this  ca.se)  )c'el  can  be  cither  i:.  tlie  gtiti  or  in  the  band  (for  x<0.2).  Imtied  the 
level  passix  rather  readily  into  tlie  rominctioii  iiaiid  near  .\a:0.2.  witliout  a  Itirge  Ftmo 
resunaiire  effect,  indieating  tiini  its  wavefunrtion  lia.s  ratlier  difT'-rcnt  ciiiiraeter  from  ilic 
.sliii'low  levels  or  llie  Gi.Asj.j.Pjj  coiidurtin:;  band  m'nie.mm.  tnecall  tiial  CaP  liav  ,n 
iii'Iirert  band  gap  til  the  .\  point  of  the  Orihunn.  /o:,.*.  in  liie  (001;  rlireetioii,.  wlide  Ga.-ts 


DEEP  LEVELS  IN  SUPERLATTICES 


3» 


Sliitllovi'derp  <ra».<iiii>n,t 

Since  the  relntive  poiition*  of  the  *u)i!itilution»)  impurit}*'*  deep  IcveU  imd  the  bitnd 
edges  of  the  iiost  determine  the  ehitraetcr  of  n  dopant,  the  physics  of  deep  levels  in  ho* 
perlnttiees  is  particularly  interesting  IT-Mj.  The  energ)'  of  a  deep  level  with  respect  to 
vacuum  is  relatively  insensitive  (on  a  scale  of  >0.1  eV)  to  the  structure  nf  the  superinttice, 
because  the  spatial  extent  of  the  c1ee|>  level  wavefunetion  is  only  of  order  sclO  A  and  only 
infrequently  comes  itito  contact  with  the  interfaces  of  all  but  the  smalle.'*  jteriod  superlat- 
tices.  Therefore  the  wsei./uncn'ott  sad  c»crjtc.«  «/  s  deep  •inpwntp  in  a  ruperlaitir.i'  do  uol 
change  nvrh  v>htn  the  pen’o.'  a/  Me  ."upcrfstliee  ehonget.  In  contrast  to  the  deep  levels 
themselves,  ihe  faperlatlke  land  tdgtf  art  readily  periurhed  by  changes  of  the  superlsittice 
periods,  exhibiting  the  well-known  effects  of  quantum  confinement:  for  rxnr.iplc.  redu.ing 
the  width  of  the  Gn.^s  layers  in  a  GaAs/Af^Gaj.^As  superlatticc  causes  n  tpinniutn-wcll 
effect,  which  rcsult>  tn  the  conduction  band  edge  of  the  superlatticc  moving  up  in  ei'erp.v 
while  the  valence  bwtd  inaximvun  dwe<nds  to  lower  energ)’.  (Sec  Fig.  (5).)  Hence  it  is 
possible  to  “band-gap  engineer"  the  p"rk»d  sixes  of  a  superlatticc  so  that,  fur  exaniplr ,  the 
small-period  superlattiee's  cunductio:i  uand  maximunt  lies  above  the  deep  level  of  a  dopant 
-  although  the  corresi>onding  large-period  superlattiee’s  conduction  band  edge  lies  belew 
the  deep  level.  This  means  that  the  doping  character  of  the  intpurity  changes  with  period 
site  |7|. 

To  be  specific,  consider  GaAs/AfojGno  3AS  supcrlatticcs  grown  iti  th-.  (001) 

direction  and  doped  with  Ga-titc  Si  near  the  center  of  the  Ga.As  quantum  wells  (Fij:  (G)  1. 
Here  Nj  is  the  number  of  Ga/is  lasTrs  in  the  superlatticc  period,  and  No  is  the  numhi-:  of 
Afg.yGaQ  3  layers.  For  the  ISxlS  large-period  superlatticc,  the  Si  dcei>  level  lies  above  the 
conduction  band  edge,  resonriut  with  the  conduction  band,  so  that  the  “extra"  Si  ch’ctron 
is  autoionixed.  making  Si"*"  -  and  the  Coulon.hic  potential  of  Si"^  binds  the  electron  lai  low 
temperature)  in  a  shallow  donor  level.  However,  for  the  small-period  2>:34  .supe:l.->itire, 
the  conduction  band  minimum  lies  at  higher  cncrg.v  that  the  Si  do?p  le'il  (i,crc:i'..'g  to 
theory  |7.S.15])  and  so  the  Si  lieep  level  lie.-,  in  the  gap.  where  it  can  trap  an  ndditioiui' 
electron  of  opposite  spin  to  its  own  electron,  becoming  St'.  In  the  thick-period  superlatticc, 
^'Ga  center  of  a  Ga.A#  well  is  a  shallow  donor,  making  the  .iiantiim  well  n-type, 

but  in  the  thin-period  2x3-1  siniciurc.  Sig;,  i.s  a  deep  irai>  wliich  tends  to  make  the  Cii.A-- 
weii  scmi-insiii.'iting.  This  t<n:i-  -.ion  from  shallow  donor  to  deep  trap  behavior  of  t!.r 
Si  substitutional  impurity,  as  a  funriion  of  Ga.As  layer  thickness  (Fig.  (7)),  is  one  of 
many  example  of  shallo»-to-dccp  iriui.-.ilions  o.  siibsiitiitional  impurities  in  scniicoi.diu  i.i: 
supcrlatticcs  Tin-  fact  that  the  mos:  common  donor  in  the  most  widely  studied  siipcrlaitm- 
exhibits  this  bebi."ior  should  be  disturbing  to  pco[>lc  wl.o  have  assumed  thai  impurities 
do  not  change  tbeir  doping  c!-:>:::ctcr  in  supi  rlattices 

Indeed,  one  ran  show  tbeorctirallj  that  in  tvn-dtmtnfional  oemtfondnrtor  ho.^tf  vvt  -y 
subslilutiona)  inipiirili  iia.s  at  lea-st  one  deep  !e\el  m  tlie  fundanientai  'oand  g.ip.  Tli.it  is. 
in  two.dinien.sions  all  in.purities  are  deep  trap®  with  levels  in  the  gap  Tins  beln-'-iur  is 
to  be  contrasted  with  the  situation  for  tlm-e  dimensions  most  iinjnn  iiies  do  not  utoducc 
deep  IcM-ls  in  tie-  gap.  Thus  mth  incrcaonj  ijiiantam  confincvi'nt.  many  fhatlow  di'iiorr  nr 
acceptor t  ran  hr  rxvrrlcd  tn  brenmr  drep  trap* nd  tn  rra.'C  prni  tding  free  cnrr-cr,«.  In  otle-r 
words,  thin  quantum  well  structures,  will  base  a  greater  tendency  to  exhibit  >■  .tii-insulai'iig 
liehavior  r.iilier  than  doping  behavior. 
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Fi(.  {S}.  Predicted  caerides  (in  eV)  of  the  AuperlAttiee  ronduction  bund  minimum 
and  valence  band  maximum  'vith  rcrpect  to  the  valence  band  maximum  of  bull:  GhAa  for 
a  GaAs/Af](Ga].j(At  [001]  auperlattice  veraut  reduced  layer  thickneaaea  and  Nj  fur 
variovu  KjxNj  jOOl]  GaAt/Af^Gaj.^At  lupcrlaUkct,  with  xaO.7  and  Nj+N2  fixed  to 
be  20,  after  Ref.  [7].  Note  the  broken  tcale  on  the  ordinate.  The  potiiioni  of  the  band 
extrema  of  bulk  GaAa  at  F,  L,  and  X  are  tlmu’n  on  the  ri);ht  of  the  Ccure.  a:  N]a20. 
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Fig.  jC}.  Iliustrntiug  llic  qunrilum-wcll  effee:  on  the  ha:;.i  c.Hp  £p.p(SL)  of  a  ,N’j>  .Vi 
Ga.As/.AfQ jG.io.s.Ai  superlailicc.  after  Ref.  J7|;  (o'  Nj=N2=lS;  aii.l  (1))  N)=2.  No'-Ol 
The  hand  odgc.s  of  tlie  supcrlatnro  are  denoted  by  rli.Mitod  hiu"..  For  this  ivlhy  com¬ 
position  the  su)>crlatlicc  gap  is  indirect  for  c.t.sc  (b|,  with  tlie  coiidiictioii  blind  edge  at 

k  =(2!r/itL)( i/2. 1/2.0).  Note  the  broken  energy  .sr.ile.  Tin-  zero  of  energ.v  is  the  iideiirc 
li.nid  m.ixiiiiiini  of  Ga.As. 


Sij,  Dt<P'i’iotle«  trsniltian 


in  N,XI0  GeAt/Al^jGijg  jA»  »uperlatlic« 


Fig.  {"}.  Illustrating  the  dcep*to-*h«ll<iw  transition  at  a  function  of  C^.^t  layer 
thieknets  Nj  in  a  GaAt/AfxGaj.xAi  NjxlO  tupcrlatiice  (SL)  with  xaO.7  for  *.  Si  impurity 
on  a  Column-Ill  titc  in  the  center  of  a  GaAs  layer  of  the  tiipcrluttice  host,  after  Ruf.  |«). 
The  conduction  band  edge  (CBM)  and  valence  band  maximum  {\'DM)  arc  indienicd  by 
light  solid  lines.  The  Si  deep  level  is  denoted  by  a  heavy  line,  which  is  solid  wlicn  the 
level  is  in  the  gap  but  dashed  when  the  level  i.<!  resonant  with  the  conduction  band.  The 
deep  level  in  the  hand  gap  for  Nj<6  is  covered  up  by  the  conduction  band  ll^  a  result  of 
changes  in  the  host  for  Tiic  impurity’s  deep  level  lies  in  the  gap  fur  Nj<C  and 

is  occupied  by  the  extra  Si  electron;  the  Si,  in  this  ease,  is  thus  ti  “deep  impurity."  For 
N)>C,  the  deep  level  lies  above  the  conduction  band  edge  as  a  resonance.  The  duugnter 
electron  ftom  the  Si  impurity  which  wh.«  destined  for  this  deep  level  is  autoionixed,  sj  tils 
out  of  the  deep  resonance  level,  and  falls  to  the  conduction  hand  edge  (light  soil,,  i.tic) 
where  it  is  subsetjucnlly  bound  (at  low  icinperaturc)  in  a  shtillow  level  associated  with  the 
long-ranged  Coulomb  potential  of  the  donor  (in<liraicd  by  the  short  tliushcd  line}.  •  is 
important  to  rcaliitc  that  both  the  deep  level  and  the  shallow  levels  coexi.-'t  and  arc  uis:  met 
levels  with  qualitatively  dilTercnt  wavcfunctioiiE  The  issue  of  whether  at.  impurity  is  “deep 
01  ‘■shallow“  IS  (letermincd  by  whether  or  not  a  deep  level  n'sociaii-d  with  >he  impurity  lies 
in  the  hand  gaji.  The  rompiiled  dee]i-shallow  iransiiioii  occurs  for  NjSiC  !au*is  \\  iiia*  the 
qiialiiati'o  physics  is  coiii)>>ct(dy  reliable,  it  is  poisibli-  that  l!ic  iraiisilioii  layer  iliicaiiess 
may  differ  somewhat  from  N]=0  m  real  supcrlaiticcs.  The  predicted  final.,:. leiital  iiiunl 
g.ip  of  llic  supcrlallice  is  iiidmct  for  I<N)<S.  .All  energies  arc  with  rcspcci  to  tiie  valence 
band  maxiiiimn  of  G.i.As. 


Dependence  of  deep  level*  «n  lietiivn 


Altliough  nhnlloW'deep  trimtitiont  th»t  :v:eur  when  a  quAntum-conBiicd  bund  cd^e 
pnssof  Ihrough  a  deep  level  ure  the  most  druniAtk  superlettiee  efTeet.  the  inlerfueet  between 
levers  in  superluttices  liuve  tlio  cfect  of  shifting  sad  splitting  deep  levels  on  a  seulc  of  order 
0.1  eV.  Fig.  {S}  shows  the  predicted  |8]  dependmee  on  site  of  an  As  vACuney's  deep  levels 
in  A  10x10  GaAs/AfgjGao  3AS  superluttiee  grown  in  the  |111]  direction.  Note  tlini  the 
valence  band  edges  of  GsAs,  Af^j  ^Gso  ^As,  imd  the  sut>erlAttice  are  aligned  uecordinr,  to 
the  measured  bsitd  offset  |1CJ,  and  that  once  this  is  done,  the  deep  levels  are  relati'-vly 
independent  of  site,  except  for  a  variation  on  the  scale  of  0.1  eV.  Tlie  bulk  p-like  wlcnce 
band  maximum  is  split  into  an  upper  x<likc  and  lower  edike  band  edge  due  to  the  reduced 
symmetry  of  the  superlattice.  This  splitting  is  reflected  in  different  energies  for  x>like  and 
edike  deep  levels,  even  if  the  parent  impurity  it  distant  from  an  interface.  The  T2  or  p*lik( 
deep  levels  split  the  most  at  an  interface  into  one  or  two  levels  corresponding  to  hybrid 
orbitals  directed  primarily  into  GaAs  (which  orbitals  have  energies  near  the  bulk  GaAs 
vacancy  energy]  and  hybrids  directed  primarily  into  AfQjGao^gAs  with  corresponding 
energies  !7,8].  In  addition  to  this  effect,  the  band  offset  between  the  two  materials  it 
reflected  in  l»th  the  p*like  deep  level  and  the  t'like  level. 

The  sdike  deep  levels  of  cation-site  Si  in  a  SxlO  GaAt/Afo,7Gao.3At  superlattice  are 
given  in  Fig.  {9},  and  show  (i)  how  undramatic  the  effects  of  the  superlattice  arC  on  the 
level  itself  and  (ii)  how  an  impurity  can  produce  a  shallow  donor  and  n-type  doping  in 
bulk  GaAs  but  yield  a  deep  trap  and  semi-insulating  behavior  in  the  superlattice.  Si  (or  a 
defect  involving  Si)  is  thought  to  be  the  DX  center  in  Af^Gai.^As  |18). 

These  shifts  and  splittings  of  the  deep  levels  as  a  function  position  in  the  superlattice 
will  manifest  themselves  experimentally  as  inhomogeneously  broadened  deep  level  energies. 

Form»li*m 

Tlie  formalism  employed  here  it  rather  simple,  and  has  been  discussed  in  Refs.  |7} 
and  |8]  in  detail.  Basically,  one  solves  the  secular  equation 

del  (  1  -  GV  )  =  0, 

where  GstE-H)'I  is  the  Green's  function  operator  associated  with  the  host  Hamiltonian 
H  and  K  is  the  deep  level  energ)'.  The  defect  potential  operator  is  V,  and  we  have  obtained 
solutions  to  this  equation  using  'iic  lipht-binding  model  of  Vogl  e:  al.  (17j.  The  bi-ic 
approarn  is  the  same  as  that  usc<'  oy  iljuimarson  et  al.  |2|  tor  substiiutiona!  impurities  in 
bulk  scnu.'onduciors. 

SUMMARY 


Deep  impurity  levels  in  semiconductor  superlattice.s  will,  on  an  absolute  cnefg.v  scale, 
be  shifted  only  slightly  and  split  (or  inhomogeneously  breadened)  from  their  bulk  semicon¬ 
ductor  energies  This  relative  insensitivity  of  the  levels  to  superlattice  geometry,  combin 'd 
with  the  supcrlattic,>'.«  property  of  being  sensitive  to  quantum  coi.Sniment  effects  on  its 
band  edg'-s  should  lead  to  slialiow-deep  transitions  and  novel  dupm  .  piuperties  of  imiiu- 
riiics  ii.  liiin  quantum  wells.  The  inhomogeneous  brondeniiig  of  deep  let  els  will  be  a  less 
dramatic  effect  on  the  cli-ctromc  structure,  but  will  nonetheless  have  important  conse¬ 
quences  in  iire.-es  such  us  cxciton  tran.';>ort  and  re.sonant  energy  transfer  of  excitons  Imuiid 
to  deep  iinpurili.'s. 


Fig.  {8}.  Predicted  energj*  level*  of  $at  Ai'V»e»ncy  in  *  (GnA»)jo(Afo.7G*0,3A»)io 
(111)  *up*rUttiee,  m  *  function  of  jj,  tW  position  of  the  vocMicy  (even  v»lue*  bf  $  co^* 
respond  to  As-sites),  *fter  Ref.  (8).  Note  the  splitting  of  the  T2  level*  st  snd  new  the 
interfaces  {fi  «  0, 20,  and  40),  and  that  the  T2-derived  vacancy  levels  lie  at  higher  energy 
in  an  AdQ  ^Gao  ^As  layer  than  in  a  GaAs  layer.  The  aj  and  e  ordering  change*  at  succes¬ 
sive  interface*.  The  aero  of  energ>‘  is  the  viUenre  band  maximum  of  bulk  GaAs,  and  the 
corresponding  s-alence  band  (VBM)  and  conduction  band  (CBM)  edges  and  deep  level*  in 
bulk  GaAs  and  bulk  Af0,7Qao.3As  are  given  to  the  left  and  right  of  the  eeniral  figure, 
respectively.  The  top  of  the  central  figure  is  the  conduction  ba,  d  edge  of  the  supcrhstiicc, 
and  the  bottom  corresponds  to  the  split  valence  band  in  the  superlattice  the  valence 
band  maximum  of  the  suj  erlattiee  being  of  c  symmetry  (p^)  and  tlic  splii-ofi  aj  (p„) 
band  maximum  lying  O.OST  eV  lo»vcr  in  energy  Th.-  Aj  level  in  the  Af^Gaj.y.As  layer 
of  the  supcriattice  is  lower  than  the  corresponding  level  in  the  GaAs  layer  because  of  tin- 
band  offset  of  0.334  eV.  The  electron  (hole)  occupancies  of  the  deep  levels  in  bulk  GaAs 
and  bulk  Afg  7Gag,3As  are  denoted  by  solid  circles  (open  triangles). 


Si  A  Ittdi 


Fig.  (9).  Predicted  Aj-drrived  dcqi  levels  of  Si  in  Ga.A.s,  lii  a  3xlli 
Ga.As/Afp  7GaQ  3AS  superlatticc,  afn-r  J.  D.  Dow,  S.  Y.  Ren.  and  .1.  Shn,.  Rrf.  |1(>' 
as  a  fui'.rtion  of  the  position  of  the  Si  in  the  supcrh-ltice.  ari<l  in  th"  bulk  iiiatcria).. 
Interfaces  correspond  to  sitc.s  /3=  0.  C,  and  2C.  Note  th.it  in  hulk  G  aAs  Si  is  ii 
Lmu  that  it  is  a  deep  trap  iti  this  suprrKMiirc.  arcordin^  lo  thr  ilicory. 
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Abstract  '... 

A  molecule  whole  electronic  itnieturw  are  oooriy  undemood.  in  lorae  ones,  can  be  thoufht  of  a 
another  perhapi  wcll-undcnitood  "hoit  chemical  molecule”  with  jubitiiutional  “impurities."  From  uhii 
point  of  view,  the  Green's  function  technique  in  solid  state  ph^aies  ean  be  extended  to  i3>’supie  the 
chanyy  due  to  substitutional  impunties  in  large  chemical  or  biological  moitcuis  of  the  froober  eiecuont: 
orbiials—the  highest  occupied  molecular  orbitals  I  homo)  and  the  Ictvest  unoccupinJ  moieedar  orbials 
I  bUMO ).  These  orbitals  are  the  ones  most  active  in  chemical  reactions,  and  i«  concepnialiy  very  aniiix' 
to  the  vaienee  band  and  conduction  band  states  of  Kmieonduetori  Substitutional  impunties  wiU  often 
introduce  new  “deep"  electronic  states  into  the  pp  between  the  HOMO  and  bUMO  stales  of  tbe  oripsal 
molecule  and  these  states  will  have  the  followmi  properties,  (r)  The  likelihood  that  an  impjriiy  uiu 
produce  a  state  in  the  pp  depends  on  the  site  of  the  impuniy.  HI)  For  impurities  on  a  parjedar  site, 
the  wav-efunetion  of  the  pp  state  will  be  relatively  independent  of  the  impurity.  I  Hi)  Only  a  smaU  sarass 
of  the  deep  level  wavei'unction  ( lyrically  tOSi)  will  lie  within  the  impurity’s  ed*  ^  bSe  la.’ps:  par. 
(about  40- 50%)  will  be  on  the  impurity’s  nearest*neighbor  atomb  (/>•)  The  derivanve  of  ike  sssrsy  of 
the  impurity  state  with  respect  to, me  ceniral<ell  impuniy  potential  rf£/^l'is  approxitsauiy  ecaal  to 
the  probability  of  the  pp>sute  electron  bei.nt  found  on  the  impurity’s  site. 


The  frontier  orbitals  of  chemical  or  biological  molecules,  namely  the  highest 
occupied  moiccuiar  orbitals  (  homo)  and  the  lowest  unoccupied  moiecu^  orbiah 
( LUMO).  arc  the  most  active  electronic  orbitals  in  chemical  reactions  and  biological 


processes  ( I ).  A  theory'  of  how  these  orbitals  change  when  the  chemical  compostion 
of  the  macromoiecule  is  altered  would  be  particularly  useful.  For  reasonably  small 
molerjics.  with  fewer  than  100  atoms,  the  stai  uard  quanium-chemisiry  ^culadoas 
can  provide  the  needed  information  about  electronic  struaure.  But  for  very*  l^e 
macromolccufes.  such  as  DNA  or  proteins,  such  calculations  are  currently  tm- 
practical.  and  so  one  must  find  an  alternative  scheme  for  predicting  how  the  froauer 
orbitals  of  large  molecules,  their  energy  levels,  and  their  electronic  prorc.'*-ss  C4.ang: 
with  chemical  composition.  Such  a  scheme  is  suggested  by  the  fact  that  inaay  Urg- 
chemical  or  biological  molecules  are  ciccirically  scmicondu'.tive  [2],  '''•’‘It  homo 
and  LUMO  states  that  .arc  conceptually  similar  sr  the  valence  bands  and  conduction 
bands  in  semiconductors.  Therefore  replacement  of  one  atom  in  a  macremolscuit 
is  analogous  to  iniroductirn  of  a  suhsiitutionai  impurity  into  a  semiconcuctor.  and 
so  the  qualitative  concepts  of  the  theory  of  deep  impurity  levels  in  scmiconou«ic:s 
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FRONTIER  ELECTRONIC  ORBITALS 


Diazabicycloocianc  can  be  considered  as  2.2.2-bicyclooci3nc  with  ‘Two  substTu- 
tional  impurity  nitrogen  atoms’*  at  carbon  sites  I  and  4  (with  two  H  atoms,  each 
of  them  bonded  to  each. of  the  two  C  atoms,  removed). 

The  Green's  function  technique,  which  so  successfully  described  the  wavefinc- 
tions  1 6)  and  the  electronic  structure  of  deep  impurity  leVels  in  semiconductors 
can  be  readily  e-ttended  to  investigate  the  change  of  the  electronic  structure  of  the 
frontier  orbitals  due  to  one  or  more  substitutional  impurity  atoms  in  a  large  chemical 
or  biological  molecule.  That  is  to  say,  for  very  large  chemical  or  biological  molecules, 
although  it  may  be  very  difficult  to  calculate  the  complete  electronid  structures  of 
the  frontier  orbitals  themselves,  the  changes  in  electronic  structure  due  to  an  "im¬ 
purity"  can  be  predicted  rather  easily.  Moreover,  in  some  molecules  so  large  that 
calculations  of  even  the  electronic  structure  of  a  single  "perfect  host"  is  impractical, 
qualitative  predictions  arc  nevenhclcss  still  feasible,  based  on  the  general  structure 
of  the  theory  of  deep  levels.  For  example,  many  e,\trcmc!y  Large  biological  molecules 
have  several  P  atoms.  The  theory  of  deep  levels  (without  any  calculation)  would 
predict  that  molecules  with  one  P  atom  replaced  by  an  .^s  atom  should  have  similar 
electronic  structures  to  the  ‘  host.”  but  that  N  substitution  might  yield  a  deep  level 
in  the  homo-lumo  gap. 

A  substitutional  impurity  produces  a  defect  potential  (9),  the  short-ranged  pan 
of  which  is  normally  associated  with  the  “deep"  impurity  levels  in  the  gap.  Thi 
states  associated  with  s-  and  p-orbitals  on  each  site  will  be  especially  penurhed  b> 
this  potential  (but  the  number  of  states  near  the  hcnK'-lumo  gap  will  not  be 
altered )  and  one  or  more  of  the  states  'ying  in  the  bands  of  the  "host"  may  move 
into  the  gap  as  a  result  of  the  penurbation.  The  gap  states  or  "deep  levels"  interest 
us  in  this  paoer..  In  par.icular,  we  are  especially  interested  in  determining  whether 
a  panicular  "impuritv"  produces  a  “deep  level"  in  the  HOMO-LUMO  gap.  and  if  it 
does,  what  would  be  the  propenies  of  such  a  le\cl--£  problem  similar  to  the  deep 
impurity  problem  in  semiconductors. 

The  approach  wc  use  to  investigate  the  wavefunciions  of  deep  impurity  levels  is 
a  semiempirical  Green’s  function  approach  (3-6)  with  a  nearest-neighbor  tight- 
binding  Hamiltonian  describing  the  electronic  band  sinictures  of  the  "hos’"  ens- 
lalline  semiconductor  or  macromolcculc.  In  the  molecular  case  the  major  di.Tcrence 
is  that  chemical  or  biological  molecules  generally  do  not  have  periodic  symmetry 
(10.11).  Here,  for  illustrative  purposes,  wc  treat  a  ralhct  small  molecule.  2.2.2- 
bicyclooctanc.  as  our  "host."  It  has  S  carbon  atoms  and  14  hydrogen  atoms.  We 
use  a  nearcsi-ncighbor  iieht-binding  Hamiltonian  to  describe  the  electronic  sirar.ure 
of  the  "host  molecule."  2.2.2-bicycloociane.  Wc  include  one  s  orbital  and  three  p 
orbitals  for  each  carbon  atom  and  one  s  orbital  for  each  hydrogen  atom.  The  resulting 
light-binding  Hamiltonian  is  an  8  X  4  -r  14  =  46-dimcnsional  mairi-x.  To  obtain 
the  tight-bindinr  parameters,  we  hi  the  ah  iniro  scf-lcaO-MO  calculation  oi  Lchn 
and  Wipff  (!2).  obiaining  the  results  of  Figu.'e  2. 

This  Hamiltonian  has  46  energy  eigenvalues,  each  of  which  can  contain  tuo 
electrons  of  opposite  spin.  The  -16  valence  electrons  of  this  "host  molccuic  v^.l 
occupy  the  ?3  lower  energy  levels,  and  the  23  higner  energy  levels  arc  unoccupied 
at  zero  lempemturc.  We  nnd  a  -t.-'c  eV  energy  gap  between  the  highest  occusita 
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[3-61  should  have  “"^'opues  in  jnacrornolecular  physics.  In  panicular.  it 
should  be  possible  "  hich  ‘-imnuritics"  or  chemical  changes  arc  likcK  to 

introduce  a  deep  level  into  the  homo-ll'mo  gap. 

In  crs-stallin®  semiconductors  such  as  Si  or  Ga.As.  an  impurity  produces  "deep" 
levels  «sociated  with  the  pans  of  the  “ccntral-ccll”  impurity  potential  at  and  near 
the  impurity  s'*®  four  such  levels  should  appear  in  the  energetic 

vicinity  of  the  fundamental  band  gap.  if  the  impurity  is  s-  and  p*bondcd:  one  s* 
like  level  and  three  p-like  levels  associated  with  the  penurbed  impurity  bonds.  The 
theory  of  deep  impurities  is  concerned  with  predicting  whether  one  or  more  of  the 
deep  levels  lies  within  the  fundamental  band  gap.  in  which  case  the  level  often 
alters  the  electron  dynamics,  by  trapping  carriers,  by  serving  as  a  center  for  non- 
radiative  recombination  of  electrons  and  holes,  or  by  providing  an  initiation  site 
fora  chemical  reaction.  Similar  phenomena  should  occur  in  macromolecules. 

Two  cases  amenable  to  treatment  by  the  theory  of  deep  impurities  in  macro- 
molecules  are; 

Several  macromolecules  with  the  same  crystal  structure  are  so  large 

■  that  it  is  practical  to  compute  the  complete  electronic  structure  of  only 
,  one  prototype,  in  which  case  the  other  molecules  arc  considered  as  the 

■  prototjpe  "host”  with  "impurities  " 

(i7)  The  macromoleculcs  are  sc  large  that  it  is  impractical  to  compute 
the  electronic  structure  of  even  a  prototype, 
in  this  latter  case  one  can  often  make  a  crude  approximation  to  the  Green’s  ‘unction 
of  the  "host"  and  then  estimate  how  the  frontier  orbitals  of  the  host  are  afTccied 
by  “impurities."  For  example,  if  take  the  rather  small  2.2.2-bicycloociane  (7) 
molecule  as  the  "host  molecule"  (see  Fig.  1 ).  then  quinuclidine  can  be  considered 
as  2.2.2-bicyclooctane  with  "a  substitutional  impurity  nitrogen  atom"  at  carbon 
site  1  (with  the  H  atom,  which  is  bonded  to  the  C  atom,  also  taken  out  [8)). 


H  H  H  H 


H  H  H  H 

Fipire  1.  Tn;  gtomcincal  s'-.-uriurr  of  ihc  2.'.'-bir\cloor:3nc  molecule  ;;';er  Rcfcrerce 
1 ' inCicaiir.g  me  sue  bocimf 
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Fifure  2.  Comparison  of  the  valence  electronic  siruciurc  or2.2.2-hicv'clooctane  ohiameO 
with  ah  iniiiii  $CF  theor)'  |l2)—and  the  ti{ihl>bindin|t  parametrieaiion  of  that  theor)- 
(•--).  The  LUMO  state  is  the  highest  energy  At,  state  and  the  iiumo  state  is  the  hignesi 

stale  in  this  figure. 


-•nsrgy  level — the  25rd  eigenvalue — and  the  lowest  unoccupied  energy  level — the 
24’.h  eigenvalue.  This  gap  is  analogous  to  the  band  gap  of  semiconductors,  and  the 
ovter  23  energy  levels  are  very  similar  to  the  valence  bands  of  semiconductors. 
v.".iie  the  higher  23  unoccupied  energy  levels  correspond  to  the  conduction  bands. 

If  a  substitutional  impurity  atom  is  introduced  into  this  host  molecule,  the  dif- 
irence  in  Hamiltonians  will  be  the  “impurity  potential"  of  the  system  [4-6).  i-iere 
'  *  consider  only  tbc  cfTcci  of  the  short-ranged  potential,  in  analogy  with  the  case 
if  crystalline  semiconductors,  for  which  the  short-ranged  potential  rather  than  the 
ong-ranged  Coulomb  potential  is  dominant  in  producing  states  in  the  cap  more 
han  0.1  eV  from  a  band  edge.  We  model  the  defect  potential  in  a  tight-binding 
tasis  of  s  and  p  orbitals  centered  on  each  atomic  site.  Following  the  theory  of  deep 
mpurity  levels  in  semiconductors  (4-6).  we  take  the  defect  potential  to  be  shon- 
•anged  in  such  a  basis  (assuming  cither  that  the  basis  states  on  cifTcreni  sites  do 
'01  overlap,  or  that  the  basis  states  are  symmctricaliy  onhogonalizcd ):  The  cfTcctivc 
-lOn-ranged  potential  may  be  expressed  as  an  “on-site"  potential  on  the  "impurity 

(faj  ^\jby  =  I'j  if  /  =  y  =  0.  and  c  -  h  =  .t. 

=  if  i  =  j  =  0.  and  a  =  h  -  .\  or  y  or  :. 

0  othertK-ise 


(I) 
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I'KON’TinR  cLECtROSiC  orbitals  ,7' 

Here !  and  y  denote  aiomic  sites  in  the  molecule,  and  «  and  ^  arc  symbols  of  orbitals 
s.  p,.  p,  .  and  p..  With  this  approximation,  the  ene  gy  level  associated  xviih  the 
impuriiv  in  the  host  molecule  will  be  determined  by  the  fclli  w-mg  secular  equation 
H.5) 

dciil  -  OT|  «  0  (2) 

Here  O’  is  the  Green’s  function  of  the  host  mol-culc: 

<«/|0'|y/0  «  I  0'</l«></i|;/>>/(/r  -  £J  (3) 

ft 

and  we  arc  interested  in  the  levels  E  in  the  gap.  Here  n  represents  one  of  the 
solutions  of  the  tight-binding  Hamiltonian  of  the  perfect  host  molecule.  For  many 
eases,  especially  when  the  impurity  has  an  attractive  shon-ranged  potential,  r, 
plays  a  more  imponani  role,  and  the  cfTcct  of  I';,  can  be  neglected.  In  this  case,  the 
energy  level  of  the  gap  state  is  determined  by. 

<0.vlO’|0.v>i;  -  1  (4) 

The  3.ssociatcd  wavefunction  is  [6] 

\^)-lC,Jia)  (5) 


where  we  have 


Co,  •  -<Os|C|0i><Oj|C|Oi>/(</<Oi|C|O.f>/</£) 


C,«- Co.</o|0’104>/<0i|0’|0i>.  .  (7) 

We  take  2.il.2-bicyclooctanc  as  an  illustrative  e.\ample  of  a  "host"  although  it  is 

not  a  very  large  molecule.  It  has  eight  car'oon  atoms,  but  only  two  of  the  carbon 

sites  arc  incquivaicnt.  namely  those  at  site  I  and  site  2  I  Sec  Fig.  I  )•  1"  Rsttre  3  we 

show  the  energx  levels  of  the  gap  state  for  site  I  and  site  2  as 

impurity  potential  I',  (sec  Eq.  (4 )).  •  • 

For  an  impurity  at  sue  I.  an  impurity  potential  of  about  -0.12  eV  t*®”  *0 

produce  an  impurity  state  in  ibv  homo-lumo  gap.  while  for  an  itnpunty  at  sits 

2.  a  much  stronger  impurity  potential.  -2.0  eV  is  needed  to  produce  «n  im^nix 

state  in  the  gap.  The  folloxving  atoms,  substituting  for  C  at  either  site,  should 

a  deep  level  in  the  gap:  the  rare  gases,  the  halogens,  and  N.  0.  S.  an<J  ®  ^ 

At  is  the  only  s-  and  p-honded  impunty  with  an  impunty  potential  berweeo  .  - 

and  -2.0  eV  when  it  substiiuics  for  C  (S  j.  and  so  .At  produces  a  dtsP 

il'.c  iiOMO-Li.'MO  gap  when  u  occupies  site  1.;  but  not  w-hen  it  occupies  sdc  -  (  J. 

In  Figure  4  we  show  that  the  on-siic  cocifieient  C'h,  for  the  j** 

1  depends  xerx  little  on  I as  determined  bx  £q.  ( 6 1.  This  is  an  examp  e^c. 

level  "pinninc"  (.x).  It  means  that  the  deep  level  in  the  M()MC'LL'M®f^P'**r*®'^‘^’‘ 

,  '  ,  .  ■  ..„,vnt  of  the  irr.pun:-.: 

with  a  wavcrirr.ction  and  an  er.ercx  mat  ..:e  almost  in.:epenoc“‘ 

...  .  .  ,  Thishapp-ens  tx-causc 

tlilieren:  impu.'inex  proouce  uim.Kt  i.ne  same  xvaxetunctions.  ' 
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Vs  (eV) 

muitS.  Thctncrtylcvaloftheiapiuteforiheiiie  K— ttindiiie2(->-)ufui>nient 
of  th«  impurity  ihon*nn|«  potcntiiil  i;,  determined  by-  Eq.  (2).  We  ehaee  the  lowest 
cneriy  lt%tl  of  the  lumo  states  of  the  host  molKule  as  the  enersy  aero. 
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the  deep  level  in  the  gap  associated  with  a  sironfily  attractive  defect  is  an  antibondins 
Kate  and  has  a  predominantly  host-like  character. 

The  charge  density  associated  with  the  deep  level  has  typically  40-50%  of  its 
charge  on  the  atomic  sites  of  the  ncarcst-ncighbors  to  the  impurity,  and  only  about 


0.0 


■20 
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Vc  (9V) 


Fiiure  4.  The  on-s:i:  u-3%;runr.ion  cj'.-fiijicr.i  C^,  pf  ’.ht  ysp  I'.sic  for  the  silt  1  3S  3 
function  of  ;.nt  impuaix  jhor.-rartft  poieniisl  I determined  by  Eq  (6) 
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Vs(eV) 

Fi|urt  5.  The  M-tlie  prabability  (>  -  •)  Z.  ICVI*.  where  i  •  s,  p„  Rr> 
corTetpentfinsMarett'neiBhberprebahiliiyt— loTthep^tMielbrihetiu  I  aefiinrJon 
of  the  impurity  thon-rente  potential  I  determined  h>'  Eqi  ( O  end  ( 7 ), 


10%  on  the  impurity's  site  (see  Fig,  S).  The  10%  on  the  impurity  site  for  this 
molecule  is  similar  to  the  corresponding  charge  densities  for  deep  levels  in  elementary 
cr)‘$talt:ne  Si.  Hertce.  the  basic  ph>‘sics  of  deep  levels  in  semiconductors  carries  over 
to  macromolecules:  the  chaige  density  in  the  impurity  cell  is  rather  smaller  than 
an  uninformed  person  might  have  guessed,  and  is  related  to  the  derisitive. 

</£/</IXCu,)*  .  .  (8) 

by  the  Hdimann-Fcynman  theorem  (131.; 

In  summary,  ihe  calculation  for  the  rather  modcrate>sized  2.2JJ*bicycloocunc 
molecule  reveals  that  deep  impurity  levels  can  be  produced  in  the  HOMO-LU,sto 
gap  with  similar  properties  to  those  found  for  deep  impurity  levdl  in  crystaJIin? 
semiconductors.  Since  this  physics  holds  for  both  the  rather  small  2.2.2*bic>xlooc. 
tane  molecule  and  for  vciy  large  semiconductor  crysuls.  it  appears  to  be  'xp’  general, 
and  hence  can  be  cxpioiicd  for  elucidating  the  electronic  structures  of  wdt  classes 
of  macromolecules  that  correspond  to  a  "host”  plus  one  or  more  subsututioaal 
“impuniies." 
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An  invene  formula  is  given  for  a  new,  generalized  Mttbius  transform  namely 

Fix.  y.  ••. )  ■  Xr.  I S  S.  I  n)pim)G{  n"x,  m*y.  i.  In  the  case  of  one  dimension  and  a«  - 1 ,  this  reduces  to  the  Chen-Mbbius 
transform  which  hu  been  applied  to  the  inverse  blacfc4)ody  radiation  problem.  For  cta-e  I  this  becomes  the  Hardy-Wright- 
Miibius  transform. 


Recently  Chen  [1]  modified  Mdbius'  number* 
theory  transform  [2]  and  showed  how  it  could  be 
used  to  solve  a  class  of  inverse  problems  in  physics, 
including  the  problem  of  determining  the  tempera* 
ture  distribution  of  a  black  body  from  iu  radiation 
spectrum:  given  a  quantity  F(x)  that  is  related  to  a 
function  G(x)  by 

G(x)~  i  F{x/n) .  (la) 

ttm  I 

the  inverse  function  F  is  given  by 

I  (lb) 

ttml 

where /i(n)  istheM6biusfunaion:j<(R)al  ifnsl. 
;t(/t)  at  ( *•  1 )'  if  n  is  the  product  of  r  distinct  prime 
factors,  and /r(/t)s0  otherwise  [2].  For  example,  in 
the  phx’sics  of  black  bodies.  G(p)  may  be  related  to 
the  radiation  spectrum  and  F(T)  may  be  related  to 
the  areal  thermal  distribution  of  the  body  (1  ].  As 
shown  by  Chen,  this  theorem  is  particularly  useful 
for  inverting  data,  for  example,  obtaining  the  un¬ 
known  areal  distribution  of  tcmintratures  on  a  black 
body  from  its  radiation  spectrum. 

Here  (i)  we  combine  Chen's  result  with  another 

'  On  leave  from  the  University  of  Science  and  Technology  of 
Oina.  Hefei  230036.  China. 

’  Present  address. 


inverse  Mdbius  transform  from  number  theory  [2]*' 
to  obtain  a  more  general  inverse  formula,  and  (ii) 
we  extend  this  generalized  formula  to  two  and  higher 
dimensions.  The  new  and  more  general  inverse 
Mdbius  transform  formula  for  one  dimension  is: 

If  we  have 

(/(.t)«  I  F{n^x)  (2a) 

Ha  I 

then  it  is  the  case  that 

/■(.v)»  I /i(/t)C(n“v) .  (2b) 

I 

Hardy  and  Wright  proved  this  theorem  for  a»  1 
("theorem  270”  (2)  *')  and  Chen’s  Mdbius  theo* 
rem  corresponds  to  ass  - 1.  We  show  that  the  theo¬ 
rem  is  valid  for  any  real  a,  and  hence  there  are  entire 
classes  of  Mdbius  transforms,  each  corresponding  to 
a  different  choice  of  o.  Our  proof  uses  "theorem  263” 
[2]  for  the  Mdbius  function, 

£/i(/i)»l,  forA'ssl, 

m\k 

sbO.  otherwise.  (3) 

Here  the  sum  is  over  all  the  different  divisors  n  of 
the  integer  k.  including  1  and  k.  The  proof  of  the  gen- 

•'  In  fact  due  to  cq.  (31 )  of  ref.  ( 1 1.  cq.  (32)  in  ref.  ( 1 1  is  a 
direct  result  of  theorem  270  of  ref.  (2). 
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eralized  formula  is  obtained  by  evtajuating  ” 

£  n{n)G(n"x)==  X  £ 

ttaal 

=  I  m“A-)I/i(«)=F(.v).  (4) 

km\  n\k 

In  this  derivation,  we  performed  the  sum  over  A:=mn 
and  n,  instead  of  over  m  and. n;  this  reordering  will 
converge  if 

I  I  |f(»n“«“v) I  =  S  d(k) \Fik<‘x) I  (5) 

convfjrges,  where  d(k)  is  the  number.pf  divisors  of 
k  [2],  Note  that  for  a=0  this  convergence  critcfion 
is  normally  not  met,  and  hence  asO  is  to  be  avoided. 
Furthermore  the  converse  of  eq.  (2)  is  true  if 

I  l\G(m‘'n‘'x)\^ld(k)\G(k«x)\  (6) 

Nail  mail  7 

converges. 

Generalization  of  forn,Ja  (2)  to  higher  dimen¬ 
sions  is  straightforward;  for  example,  in  two  dimen¬ 
sions  for  F(x,  y)  of  the  form 

<7(a-..v)=  f  lF{mOx,nfy),  (7a) 

mal Hml 

we  have 

nx,  .V) »  X  f  y(k)y{l)  Gik^x,  l>y) .  (7b ) 

km  \  im\ 

k  sufficient  condition  for  this  inverse  formula  to  be 
valid  is  the  convergence  of 

Hll\nm’>k^x,n'‘l*y)\ 

k  m  t  N 

-II</(<VO)|f(/".v,/r)|.  (8) 

The  new  inverse  Mdbius  formula  (2)  can  be  applied 
to  problems  in  mathematics  and  physics,  and  a  con- 
venii'nt  choice  of  the  indices  a,  jJ, ...  can  be  selected 
to  best  suit  each  problem. 

The  onc-dimensional  version  leads  to  the  new 
mathematical  identity 


r^:=-  £  fi(n)n-'’  for p>  L-,.  (9) 

v.tPr  <1-1 

wliw  C(p)  =  -?.ii«“'’  |s  the.Riema'nn  zeu  func- 
'tion.  It  also  iMdsio^the  identities 

I  \fiOi)Hnx/n)f^(Kx/n)-l)  (10) 

and  ;  ' 

-b2  £  («V-4»rjt-)-'«),  (M) 

where  we  have  INn < nx<  (2S+  2 )«,  Jo  is  the  cylin¬ 
drical  Bessel  function.  f^{x)-co\\i{x),  and 
/_(-Y)=cot(.Y).  and  we  have  used  the  identities  [31 

/i()LY)=:-^,+^  £  {x-±n-)-', 

ItX  It  n.l 

and 

£  Mrtx)^-\  +  T.+2  £  (.Y*-4m-K^)-''* 

n—  I  -  »t.  I 

for  2Nn<x<  (2/V+2)Jt. 

There  are  undoubtedly  many  more  u.cs  for  this 
inverse  formalism  in  both  mathematics  and  theo¬ 
retical  physics  that  will  be  uncovered  as  it  is  used 
more. 

We  are  grateful  to  the  U.S.  Office  of  N'aval  Re¬ 
search  and  the  .^ir  Force  Office  of  Scientific  Re¬ 
search  (Contract  Nos.  N00014-89-J-I136  and 
AFOSR-89-0063)  for  their  generous  support. 
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We  compare  the  forces  generated  by  various -ans)e>depehdeii,t  po(entiali(>and  'by.a6  initio  band- 
structure  calculations  with  a  lirniied  number,  ori  points.  Iti  ceJls.with  32  and  54  atoms  we  nnd  sub¬ 
stantial  errors  for  all  of  the  angle-dependent-force  models,  as  well  as  for  the  band-structure  forces 
with  very  few  k  points. 


Recentiy  there  have  been  a  nuntber  of  computer  simu¬ 
lations  pf  amorphous  solid  Si  (o-Si)  and  of  liquid  Si  (/-Si) 
using  tWo  distinct  methods: '(i),  empirically  , fitted  anglb- 
dependent  forces  (ADF)  (Refs.  1-8),  and  (ii)  first- 
principles  band-structure-derived  forces /BSF).’*'°  Finite 
Si  clusters  and  reconstructed  Si  surfaces  have  also  been 
studied.  Similar  studies  have  been  performed  Or  are  be¬ 
ing  performed  for  other  materials."  Both  of  these 
methods,  .^DF  and  3SF,  have  inaccuracies  that  have  not 
been  systematically  evaluated  for  liquid  and  amorphous 
geometries.  The  ADF  models  are  fit  to  a  fairly  small  da¬ 
tabase  so  that  one  might  question  how  well  they  predict 
the  forces  in  an  amorphous  solid  of  liquid.  This  is  espe¬ 
cially  true  since  relatively  small  errors  in  total  energies 
can  lead  to  rather  large  errors  in  forces.  The  BSF 
methods,  which  use  first-principles  molecular-dynamics 
simulations  within  a  supercell  geometry  are  not  fit.  How¬ 
ever,  becau.se  of  the  computational  comple.xity  in  actually 
implementing  them,  they  have  a  different  set  of 
difficulties.  One  car,  ask  how  important  cell  size  is 
and/or  what  the  dependence  is  on  the  band  curvature  or 
number  of  k  points  in  the  Brillouin  zone.  In  this  com¬ 
munication,  we  Briefly  investigate  these  questions.  For 
the  .ADF  models  we  have  limited  our  investigations  to 
the  potentials  of  Stillinger  and  Weber  (SW).'  the  two  po¬ 
tentials  of  Biswas  and  Hamann  [BH-1  (Ref.  S)  and  BK-2 
(Ref.  6)1,  and  the  model  of  Carlsson,  Fedders,  and  Myles 
(CFM).’  Concerning  the  BSF  method,  the  cell  size  and 
k-point  sampling  weic  investigated  using  the  ab  initio 
total-energy-molecular-dynamics  scheme  of  Sankey  and 
Niklewski.  '• 

For  the  BSF  model,  we  find  that  the  use  of  only  the  F 
(k^O)  point  in  .samples  of  from  32  to  54  atoms  leads  to 
typical  errors  of  0.3-0.4  eV/A.  These  are  small  in  com¬ 
parison  with  typical  force  scales  in  Si  that  are  of  order 
several  eV/A.  However,  errors  of  this  magnitude  caused 
one  sample  to  rem:un  amorphou.s  upon  annealing,  when  a 
more  accurate  calculation  led  to  crystallization.  Further, 
we  have  generally  found  fewer  defects  in  .samples  an¬ 
nealed  with  a  more  accurate  procedure.  This  suggests 
that  the  energy  barriers  lor  necessary  forces'  for  crystalli¬ 
zation  and  defect  healing  are  much  lower  than  the  natti- 
ra!  scale  in  Si.  Increasing  the  number  of  points  to  eight 
can  still  lead  to  non-negligible  error.  Ti:  •  .ADF  models 
all  yield  errors  several  times  larger  than  the  F^pomt  BSF 


method.  . 

We /first'  examine  the  accuracy  of.  the  BSF  method. 
Ab  initio  molecular-dynamics  simulations  usually  in¬ 
volve  the  use  of  “supercells,"  because  it  is  computational¬ 
ly  impossible  to  simulate  a  cluster  large  enough  to, ignore 
surface  states.  In  addition,  ‘supc.'cells  are  convenient 
since  it  is  quite  straightforward  to  use  the  apparatus  of 
crystalline  solids  to  simulate  approximately  an  infinite 
system.  One  chooses  a  particular  Bravais  lattice  with  a 
large  unit  cell,  and  periodically  repeats  the  cell  through 
■  all  space.  The  aim  of  this  cor.:>truction  is  to  eliminate 
spurious  surface  effects,  while  providing  a  cell  large 
enough  to  include  a  representative  sample  of  local  mi- 
crostructurcs.  The  properties  of  this  infinite  system  are 
then  easily  handled  with  k-spuce  methods.  A  crucial  as¬ 
pect  of  this  procedure  is  the  choice  of  geometry  of  the 
cell. 

If  the  supercell  is  large  enough  so  that  intetcell  interac¬ 
tions  are  negligible,  we  are  still  left  with  the  tequirement 
of  handling  the  “supercell  bands"  properly.  Because  the 
supercell  lattice  constant  is  much  larger  than  typical  in¬ 
teratomic  spacings,  the  supercell  bands  are  much  nar¬ 
rower  than  for  a  system  with  s  smaller  unit  cell.  In  this 
paper  we  examine  the  effects  of  curvature  and  dispersion 
of  the  supercell  bands  on  the  atomic  forces,  in  serne  com¬ 
monly  used  supercell  geometries.  The  systems  we  treat 
arc  disordered  and  have  a  finite  electronic  density  of 
states  at  the  Fermi  level.  The  Fermi  level’s  crossing  the 
supercell  bands  is  e.xpcotcd  to  enhance  the  effects  of  the 
band  dispersion  on  the  total  energy.  If  ail  bands  were 
completely  filled  or  empty,  then  only  the  center  of  gravity 
of  a  band,  and  not  its  width,  would  affect  the  total  ener¬ 
gy.  This  is  not  true  for  puriially  filled  bands.  Typical 
.ADF  method.-!  yield  15CJ-  or  m.orv  defects,  which  yields  u 
.sizable  density  of  states  at  the  Fermi  level. 

It  has  become  a  common  practice  in  ab  initio 
mulecular-d\ mimics  .simulations  10  use  one  k  point  in  the 
Brillouin  zone  to  calculate  energies  and  forces.  In  efl'eci 
this  assumes  that  the  supercell  bands  are  flat.  For  com¬ 
putational  reasons.  >he  F  point  tk»OJ  is  usually  selected 
for  this  purpose.  To  evaluate  this  procedure,  vve  have 
calcuiated  the  supe.'cell  bands  for  a  54-atom  fee  structure 
that  is  highly  disordered.  There  is  signifieiint  dispersion 
in  the  hands  for  the  5-i-.iiom  ceil,  wiili  widths  of  order  0. 1  j 
eV.  \Vc  performed  .i  sin-.ilar  eaiculalion  for  the  54-iitvMr  | 
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FIG.  1.  Radial  distribution  Tunciiun  of  samples..-!  and  £,  in 
reduced  units  of  </=  2.35  A. 

ever,  starting  with  the  same  near  equilibrium  sample, 
when  it  was  annealed  further  with  four  k  points^  it  be¬ 
came  a  crystal.  Thus,  the  structures  that  emerge  are  very 
sensitive  to  the  interatomic  forces.  Although  (as  men¬ 
tioned  above)  the  annealing  sequence  that  was  used  is.cer- 
tainly  far  removed  from  the  actual  physical  growth  pro¬ 
cess  for  fl-Si,  we  see  no  reason  why  the  physical  growth 
process  should  be  less  sensitive  to  the  forces.  Sample  L 
was  obtai.ied  in  the  same  fashion  as  sample  A,  but  the 
procedure  was  stopped  when  the  kinetic  energy  was  still 
several  thousand  degrees  per  atom.  Its  radial  distribution 
function  is  shown  in  Fig.  1,  and  is  basically  liquid-like. 
However,  the  radial  distribution  function  is  quite 
different  from  the  e.tperimental  one.  ”  This  is  not  a  seri¬ 
ous  problem,  since  these  samples  were  created  in  order  to 
test  the  potentials  and  not  to  simulate  real  /-Si.  Thus  our 
two  samples  include  a  sample  that  is  close  to  being  a-Si 
and  a  sample  that  is  still  quite  liquid-like.  Our  results  on 
the  forces  are  summarized  in  Table  I  (for  the  32-atom 
sample  A)  and  Table  II  (for  the  54-atom  sample  L).  We 
have  defined  A/*  as  a  root-mean-square  deviation 


where  F(/)  is  the  computed  force  on  atom  i,  F,(/)  is  the 
e.xact  (16-point  BSF  method)  force  on  atom  /,  and  N  is  the 
number  of  atoms  (32  or  54).  The  quantities  6,F  in  the 
tables  are  the  errors  in  the  forces  for  the  3  or  5  worst 
atoms.  The  quantity  AF'  is  the  same  as  AF  c.xcept  that 
the  worst  10%  of  the  atoms  were  taken  out  of  the  sum 
and  thus  .V  is  29  or  49.  The  rms  force  on  the  .atoms  in 


Since  sa.mple  A  had  been  e.xtcnsisciy  annealed,  the  e.x- 
ac;  force.s  are  sniall  and  thus  the  r-;!ative  error  in  the 
fc:,  *s  is  rather  large  compared.  lOAii'i  rms  total  force.  I: 
is.mtcresting  to  note  Uve  atom’s  bn  which  tfi'e  errors  were 
the  .'greatest'  are  highjy  -correlated  not  only,  among  the 
Various  ADF  methods  but  also  between  the  ADF 
methods  and  the  P-point  BSF  method.  For  this  reason 
"afefe  numbers  of  the  worst  atom  arc  included  in 
parentheses  in  Tables-I  and  II.  Atoms  no.  1'2.  24,  and  29 
allhave  one  bond  angle  of  less  tluan  90*.  atom  no.  31  has 
bond  angles  between  and  127*,  an^  atom  no.  9  has  a 
dangling  bond.  Evidently  almost  allbf  the  methods  have 
more  trouble  with  badly  strained  bonds  than  with  dan- 
"glihg  borids.  Further,  from*  the  tables,  one  ban  see  that  a 
substantial  part  of  the  rms  error  comes  from  10%  of  the 
atoms. 

For  the  54-atorri  more-liqiiid-likc  sample  L,  we  note 
that  all  of  the  ADF  methods' predict  forces  that  arc  in  er¬ 
ror  by  about  the  magnitudb'bf  the  exact  forces,  while  the 
BSF  method  is  much  betteK  even  with  one  k  point.  Since 
we  have  found  no  correlations  between  the  methods  on 
which  atoms  were  treated  worst,  we  have  not  included 
the  atom  numbers  in  this  case.  We  have  qlso  noted  some 
dynamical  differences  between  the  methods.  For  exam¬ 
ple.  if  one  plots  rHt)  versus  t  with  the  e.xact  BSF  method, 
one  obtains  a  reasonably  straight  line  indicating  the  ex¬ 
pected  diffusive  behavior.  However,  with  only  one  k 
point,  r’  is  supralinear,  suggesting  a  partially  ballistic  be¬ 
havior. 

Some  further  general  trends  are  evident  as  one  goes 
from  the  least-accurate  methods  (the  ADF  models)  to  P- 
point  BSF  calculations  and  finally  to  exact  BSF  calcula¬ 
tions.  One  trend  that  is  evident  among  all  of  the  methods 
is  that  more  accurate  methods  tend  to  favor  crystals  and 
methods  with  larger  errors  tend  to  favor  disorder.  As 
noted  earlier,  we  found  samples  that  had  converged  to  an 
amorphous  sample  using  the  P-point  BSF  method  but 
that  became  crystalline  when  annealed  further  with  more 
k  points.  Further,  it  has  become  clear  from  other  investi¬ 
gations'*"^'  that  the  ADF  methods  invariably  give  very 
large  (or  order  15-20%)  concentrations  of  defects. 
•Another  trend  is  evident  on  a  shorter  time  scale.  This  is 
the  observation  that  distortions  from  perfect  crystal  or 
from  tetrahedral  coordination  are  relaxed  faster  and 
more  completely  with  the  better  methods. 

In  order  to  understand  better  the  origins  of  the  errors 
caused  by  using  a  finite  number  of  k  points,  wc  have  per¬ 
formed  some  additional  calculations.  It  is  well  known 


TABLE  1.  Comparison  of  the  errors  in  the  32-atom  cluster,  sample  .•!.  The  root-mcan-square  error 
AFis  deiined  in  the  text  ai'.J  6,Fare  the  errors  in  the  forces  for  the  worst  three  atoms.  The  numbers  in 
parenihe.ses  refer  to  the  one  atom  number  in  the  cluster  an:!  all  units  are  cV/.A. _ 


Method 

AF 

AF' 

f>,F 

b:h' 

5,F 

1)S  tk=0l 

0.34 

0.23 

0.74.12! 

0.62'21t 

0.541291 

DIM 

1..‘4 

1.41 

2.SO.-1- 

2.2"  24> 

•>  •'SllOl 

mi-2 

0.74 

0.57 

2.0hil2’ 

l.4')-3ll 

l.2Si2')i 

S\V 

I.O'f 

O.hn 

.Vi 2  l.'ti 

2.4Sii:! 

CFM 

1  01 

O..XI 

2.43il2- 

2.i)"'31. 

1.S4!24‘ 

'-TABIi^Ejnl.^EWot-)  in  the  fo 
denned, jMhejl|xt;a;hd!6^^ 

f«>Tof-.tns;5.~atjiTO>clu.«cri.!.arrif»ie  L.  fhe;root-nK-a:.-square  errur  A  A  v. 
,fhe  erfbrsMriaheAorcesfdftih'r.w6rsi5five,a;rins.  All.unitVare.e'V.A. 

hfethqd;  -.: 

2^''.  'll'". 

'  's.F  lA  ,  biF  h.F 

IJS'»i45k>poimsi 

sb)15 

0il4 

WM 

0.2?  'i 

0.20  ‘  0.20  0.20' 

DS(-A-=t'0i;  .  *■'  ' 

■0.43 

.  0..39 

0.74;^  0.72  0.68 

BHU 

•  ,15;'54' 

.  748 

65:65 

61.41 

■?l..t8  .  30.63  17.11 

2.45. 

2,19. 

440 

4.40 

4.39?  ,  ^.OU  3.85 

SW 

4101. 

3.52. 

8.45 

7.95 

0.91  6.40  6..M 

GFM  1 

6;98. 

r4;54.-,,. 

25.69 

23.99 

12,96  10.97  9.09 

that  in  order  to  obtain  accurate  Brillo.uin-zone  ,(BZ)  in¬ 
tegrations,  one  needs  more  k  points  for  rnaterials  with, .a 
finite  density  of  states  at  the  Fermi  energy  than  for  ma¬ 
terials  with  a  gap  at  the  Fermj  energy.  To  establish  the 
importance  of  thisi  effect,  we  performed  calculatipns  for . 
the  energies  and  forpes  on  crystals  that  were  allowed  to 
disorder  slightly  by  small;  randotn  motions.  In  case  I  we 
allowed  rms  excursions  of  roughly  0.1  A.  This  simple  re¬ 
tained  a  gap  of  about  1  eV.  In  case  II  we  allowed  rms  .ex¬ 
cursions  of  roughly  0.2  A,  and  the  gap  barely  closed. 
Both  samples  contained  54  atoms.  In  both  cases  the  er¬ 
ror  obtained  by  using  only  one  k  points  was  0.26  eV  per 
atom.  However,  the  rms  error  in  the  force  in  case  I  was 
0.13  eV/A,  while  the  error  was  0.53  eV/A  in  case  II. 
This  suggests  that  the  majority  of  the  error  is  in  the 
bands  crossing  the  Fermi  surface.  Further,  we  note  that 
the  energy  differences  between  different  crystal  structures 
arc  often  less  than  0.26  eV  per  atom. 

One  can  understand  this  result  in  the  following  way. 
We  view  the  use  of  different  numbers  of  k  points  in  the 
BZ  as  corresponding  to  different  interpolations  for  £(k) 
versus  k  between  the  chosen  k  points.  The  approximate 
£(k)  then  oscillates  about  the  exact  £(k)  over  a  distance 
^k  in  k  space  and  with  an  amplitude  proportional  to  Ak. 
For  bands  crossing  the  Fermi  surface  the  integral  is  cut 
short,  leading  to  errors  of  order  (Ak)*  (a  length  times  a 
height).  The  first  factor  of  Ak  comes  from  the  number  of 
wave  functions  that  arc  affected,  and  the  second  one 
comes  from  the  distance  of  these  wave  functions  from  th^ 
Fermi  level.  On  the  other  hand,  the  error  in  the  charge 
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density,  which  determines  the  electric  field  and  thus  the 
forces,  is  of  ordeV  Ak,  Thus  the  total  energy  is  in  error  b\\ 
terms  of  order  (Ak)N  while  the  forces  arc  in  error  by 
terms  that  are  linear  in  Ak. 

In  conclusion,  we  have  shown  that  the  ADF  methods 
studied  yield  forces  for  badly  disordered  samples  far  from 
equilibrium  that  are  in  error  by  about  the  mfignitude  of 
the  forces.  In  terms  of  absolute  error,  these  methods  arc. 
much  better  for  samples  that  are  nearer  to  equilibrium 
but  are  still  only  qualitatively  correct.  Whether  they  can 
yield  reliable  defect  or  surface  structures  is  open  to  ques¬ 
tion.  The  r-point  BSF  method  for  small  clusters  is  worst 
when  far  from  equilibrium.  For  such  clusters  the  forces 
arc  in  error  by  about  0.4  eV/.A  and  in  error  by  about 
0.1 -0.2  eV/A’  when  using  four  k  points.  The  absolute 
error  m  the  forces  when  using  only  one  k  point  is  rather 
small  for  nearly  converged  samples,  but  this  small  error 
can  make  the  difference  between  att  amorphous  and  crys¬ 
talline  sample. 
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function,  but  rather  perform  the  integrals  numerically 
from  a  wave-function  tabulation. 

The  total  energy  in  Eq.  (1)  can  be  conveniently  rewrit¬ 
ten  as  a  sum  over  eigenvalues  of  the  single-particle  Ham¬ 
iltonian, 


E,Jn)==Eos{n  -i-[U,i-V„{n,n)]+bUJn) ,  |3a) 

where  E^{n)  is  the  "band-structure”  energy  and  is  given 
by  a  sum  of  one-electron  eigenvalues  over  occupied  states 
i, 


£bs<«)=2  2e,  .  (3b) 

I 

The  quantities  e,-  are  the  single-particle  Hamiltonian  ei¬ 
genvalues  that  satisfy  the  Schrodinger  equation  in  matrix 
form 


2  ^2;v^(v,a')”e;  "2  S“;;‘fl,(v,a') ,  .  (4a) 

v,a’  v,a' 

where  the  single-partic!?  Hamiltonian  matrix  elements 
are 

and  the  overlap  matrix  is  ' 

S“:e‘»<i|:^0(r-r„)|^P^0(r~v)>  .  (4c) 

The  standard  single-particle  LDA  Hamiltonian  operator 
is 

^(«)“^+2if'ton(r-r„)+K„,(r-r,)] 

<1 

.  (5) 

where /ij,j(fl ) =d(n  £,j( fl  )j/</n  is  the  exchange-correlation 
potential.  The  other  tei  ms  in  Eq.  (3a)  arc  the  "short- 
ranged"  repulsive  potential 


(6) 


and  an  exchange-correlation  correction, 

5^xe“/«(r)[e„(n)“/ixc(«))«(^''  •  (7) 


where 


tt  (r)  2^«lom(E  • 

n 


(9) 


The  neutral-atom  density  for  Si  is  taken  to  be  a  spherical¬ 
ly  symmetric  s'p^  configuration.  To  first  order  in  5«,  the 
energy  functional  of  Eq.  (3a)  simplifies  to 


where  Eb's’  is  the  “band-structure"  energy  determined 
from  the  single-particle  Hamiltonian  of  Eq,  (5)  with  n  re¬ 
placed  by  n,  viz.,  h{n%  This  functional  neglects  terms 
of  order  5/i  \ 

Our  approach  is  to  use  this  approximate  total-energy 
functional  rather  than  the  fully  self-consistent  total- 
energy  functional.  The  advantages  of  this  method  are  (1) 
the  electronic  eigenvalue  equation  only  needs  to  be  solved 
once  for  each  atomic  configuration  instead  of  ~  10  times 
as  in  a  fully  sclf-consistent  calculation  and  (2)  four-center 
Coulomb  integrals  do  not  need  to  be  evaluated  since  they 
appear  only  in  the  bn  *  terms,  which  are  neglected.  This 
non-self-oonsistcnt  method  has  been  tested  by  Harris,^' 
Polatoglou  ei  a/.,^'  and  Rend  et  al.’^  on  metals,  semicon¬ 
ductors,  and  an  ionic  compound  (NaCl)  and  it  is  found  to 
give  surprisingly  good  agreement  with  self-consistent  cal¬ 
culations. 

Out  fourth  major  approximation  involves  the  evalua¬ 
tion  of  the  matrix  elements  of  the  l^arious  terms  that 
make  up  ^(/t®).  The  kinetic  energy,  overlap,  and  nonlo¬ 
cal  part  of  the  pseudopotential  are  easily  evaluated  exact¬ 
ly  (numerically)  for  each  geometry.  A  table  of  results  on 
a  fine  grid  of  separations  is  constructed  so  that  the  matrix 
elements  for  any  separation  between  the  atoms  is  accu¬ 
rately  interpolated.  The  three  center  matrix  dements  of 
the  neutral-atom  potential  given  by 


are  accurately  approximated  by  an  r-dependent  expan¬ 
sion  in  multipole  moments.  The  nonlinear  exchange- 
correlation  matrix  elements  are  calculated  within  the 
average  density  approximation  described  in  Ref.  23.  All 
matrix  elements  are  calculated  in  real  space. 

The  forces  are  determined  by  differentiating  the  total 
energy  given  by  Eq.  (10),  F,.*-0£ji,'/3r,..  The 
most  difficult  term  comes  from  the  band  structure.  Its 
derivative  is 


The  Schrodinger-like  Eq.  (4a)  follows  from  a  variation¬ 
al  principle  for  the  total  energy  of  Eq.  (1).  The  single- 
particle  Hamiltonian,  Eq.  (5),  defines  a  problem  requiring 
a  self-consistent  solution.  To  avoid  the  difficulties  associ¬ 
ated  with  iterating  to  self-consistency,  we  adopt  an  ap¬ 
proximation  suggested  by  Harris.^’  TTius  our  third  major 
approximation  is  to  consider  a  sum  of  neutral-atom 
spherical  charge  densities  as  a  zero-order  approximation 
to  the  self-consistent  density  of  the  cluster,  and  to  keep 
only  first-order  changes  from  this  density  in  the  energy 
functional.  Thus  we  write 

«(r)=n°(r)+Srt{r) ,  (8) 


-3£‘bs 


or, 


=  -2 
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(iKC) 
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H.v 


a.  a 


-1 

M.V 

a.u 


y/iv 


,3(/j®);;"' 

_ £an' _ tiL. 


3r, 


'/(V 


3r, 


(12) 


where  p®,  and  are  the  density  and  energy-density 
matrices,  respectively. 
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similar  to  that  recently  proposed  by  Foulkes  et  al?* 

In  this  paper,  we  use  this  method  to  determine  the 
equilibrium  structures  and  vibrational  spectra  of  small  Si 
clusters  (Si„  with  «=2-7).  The  technique  of  molecular 
dynamics  is  employed,  in  which  the  many-body  classical 
equations  of  motion,  F=m  dh/dt^,  is  solved  in  time, 
and  the  subsequent  motion  of  the  atoms  is  determined 
from  some  set  of  initial  conditions.  Equilibrium  struc¬ 
tures  were  found  by  simulated  annealing  and  dynamical 
quenching  (in  which  the  system  is  periodically  quenched 
by  remo-'ing  kinetic  energy  from  its  atoms).  In  finding 
the  vibrational  spectra,  the  velocity  autocorrelation  func¬ 
tions  arc  Fourier  transformed.  The  theory  and  tech¬ 
niques  described  in  this  paper  are  immediately  transfer¬ 
able  to  condensed-matte*  systems,  which  will  be  dealt 
with  in  future  publications.  , 

II.  ELECTRONIC-STRUCTURE  THEORY 

In  this  section,  we  briefly  describe  the  theory  used  to 
determine  the  electronic  structure  and  to  obtain  the  total 
energy  of  these  systems  and  forces  between  atoms.  We 
will  describe  four  major  approximations  that  simplify  the 
electronic  structure  tremendously,  so  that  medium  scale 
( -- 100  atoms  or  more)  molecular-dynamics  simulations 
can  easily  be  performed.  Only  a  brief  description  of  the 
theory  will  be  given,  as  a  more  complete  description  can 
be  found  in  Ref.  23. 

The  theoretical  foundation  used  is  density-function 
theory.  Within  this  rigorous  ground-state  theory,  we 
make  our  first  major  approximations.  The  approxima¬ 
tions  are  the  use  of  the  Hohenbsrg-Kohn-Sham  local- 
density  approximation  (LDA)  and  the  pseudopotential 
approximation.  The  LDA  replaces  the  exchange- 
correlation  energy  functional  by  a  local  function  of  the 
density.  We  use  the  local  exchange-correlation  function¬ 
al  of  Ceperley  and  Alder  as  parametrized  by  Perdew  and 
Zunger."  The  pseudopotential  approximation  replaces 
the  core  electrons  by  an  effective  potential  that  acts  on 
the  valence  electrons.  Accurate  nonlocal  (angular- 
momentum-dependent)  pseudopotentials  of  the  norm- 
conserving  Hamann-Schliiter-Ciiiang  type  are  used.^* 

In  the  LDA  and  nonlocal  pscudopotential  approxima¬ 
tion,  the  total-energy  functional  is  given  by 


^.o.  =  r,.  +  / 


,  r  n  *■ )  .  I 

+  -—J  - - +e„i 

2  r-r'  ' 


,(«) 


n(r)£/V  , 


(1) 


where  the  individual  terms  arc,  respectively,  the  kinetic 
energy,  the  ionic  local  pscudopotential  of  the  ion  at 
the  nonlocal  pseudopotential,  the  electron-electron  Har- 
tree  repulsion,  and  the  exchange-correlation  energy, 
which  is  a  functional  of  the  electron  density  «. 

The  electronic  energy  eigenstates  t/'/  aff  the  expanded 
in  a  tight-binding-like  linear  combination  of  pseudo- 
atomic-crbitals  (PAO's): 


!;&/(/•)>=  :Eof(,s.a)l^i:^°(T-rJ)  .  •  (2) 

/i.a 

The  PAO’s  are  self-consistently  determined  eigenfunc¬ 
tions  of  the  valence  electrons  of  the  free  atom  in  the  non¬ 
local  pseudopotential  approximation  and  are  nodeless. 
The  pseudoatom  contains  only  the  valence  electrons,  so 
that  the  orbital  types  for  Si  are  s,  Px,Py,  and  p,. 

Our  second  major  approximation  is  motivated  by  the 
need  to  reduce  the  range  of  interaction  between  atomic 
orbitals,  and  hence,  greatly  reduce  the  number  of  neigh¬ 
bors  each  atom  (or  pair  of  atoms)  interacts  with.  This 
reduction  leads  to  a  corresponding  reduction  in  computer 
time  required  for  the  calculations.  We  do  this  by  slightly 
exciting  the  PAO's  by  imposing  the  boundary  condition 
that  they  vanish  at  and  outside  a  predetermined  radius  r^. 
(As  the  atom  approaches  the  ground  ctate.)  The 

value  of  Cf  that  we  use  here  and  in  previous  work  is 
S.Ofl^.  The  motivation  for  this  value  is  that  it  rigorously 
yields  a  third-neighbor  model  for  matrix  elements  of  the 
single-particle  Hamiltonian  in  crystalline  Si.  We  have 
found  that  our  results  are  not  critically  dependent  on  the 
precise  value  of  r^  as  long  as  r^  is  not  too  small.  The  ki¬ 
netic  energy  associated  with  the  confinement  of  the  elec¬ 
tron  in  the  atom  begins  a  sharp  increase  at  values  of 
slightly  less  than  Sog.  A  plot  of  the  Si  s-orbital  wave 
function  in  the  pseudopotential  approximation  for  vari¬ 
ous  values  of  is  shown  in  Fig.  1 .  The  curve  for  r,  “  iog 
is  virtually  identical  to  the  ground-stat®' wave  function. 
Notice  that  for  r.  -Sog,  the  wave  function  in  the  bond¬ 
ing  region  is  well  represented,  but  the  long-ranged  tail  is 
eliminated.  In  the  multicenter  integrals  needed  in  this 
work,  we  do  not  fit  the  wave  function  to  any  analytic 


FIG.  1.  The  s  pseudo-atomic-orbital  of  Si  using  various 
valuc>  of  r. .  The  wave  function  for  r,  “Rofl  i-.  wry  close  to  the 
pscudouiom  ground  state.  The  bondi'ig  region  '.defined  to  be 
half  the  nearest-neighbor  distance  in  bulk  Si)  is  shown  by  the 
vertical  arrow. 
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function,  but  rather  perform  the  integrals  numerically 
from  a  wave-function  tabuiation. 

The  total  energy  in  Eq.  (D  ean  be  conveniently  rewritf^ 
ter.  as  a  sura  over  eigenvalues  of  the  single-particle  Ham¬ 
iltonian, 

Et^{n)-E^{nt[U,i-U„{n,n)]+bU^e{nl,  {3a) 

where  Ebs^'*)  ^  '‘band-structure"  energy  and  is  given- 
by  a  sum  of  one-electron  eigenvalues  over  occupied  states 
i. 


ET^{n)^2'2,ti  •  (3b) 

/ 

The  quantities  S;  are  the  single-particle  Hamiltonian  ei¬ 
genvalues  that  satisfy  the  Schrodinger  equation  in  matrix 
form 


2  2  S^^‘a,{v,a') ,  .  (4a) 

>',0'  ViO* 

where  the  single-particle  Hamiltonian  matrix  elements 
are 

*;:5'-<^Hr--r,)|Mn)1^^0(r-r„.)>  ,  (4b) 

and  the  overlap  matrix  is  ' 

•  (40 

The  standard  single-particle  LDA  Hamiltonian  operator 
is 

(5) 

where )]/<//!  is  the  exchange-correlation 
potential.  The  other  terms  in  Eq.  i3a)  are  the  "short- 
ranged"  repulsive  potential 


Z,'  •'  Ir-r'l 


«(r') 


(6) 


and  an  exchange-correlation  correction,  bU^, 

5(^w“/«(r)[eic(«)~Mxe(«l]<(^'' •  (7) 


where  ''  < 

«®(r)“2«.<«n(r-f«l.*'  '  (9) 

n 

The  neutrahatom  denOty  for  Si  is  taken  to  be  a  spherical¬ 
ly  syihmetric  x:^p.  ^configuration  5/i,.the 

energy  functional  of  1^;  (3a)  simplifid  to 

whe^c “i^4nd,-structu^^^^  ®^®rsy  <l®l«t:tnlned 
from  the  sihgie-phhlcIe^Hamiltonian  witH  fl:re- 

placed  by  n^^  vi2;,;fi(h®).-:'11»is-fuhctional  neglect^ 
oforderM*;  '  ■  .r 

Ourapproach  isvtp  use^this  approximatentotal-energy 
functional  rather  than  the  fully' teif-ednsUtent  total- 
energy  functional.  The  advantages  of  this  method  are  (1) 
the  electronic  eigehVaiuP|quatlph'Phly  needs  to  be  solved 
once  for  each  atomic  configuration  instead  of  ~  10  times 
as  in  a  fiilly  self-consistent  calculation  and  (2)  four-center 
Coulomb  integrals  do  not  heed  to  be  evaluated  since  they 
appear  only  in  the  5n^  terms,  which  are  neglected.  This 
non-self-oonsistent  method  has  been  tested  by  .Harris,^' 
Polatoglou  ei  al.}*  and  Read  et  al.^  on  metals,  semicon¬ 
ductors,  and  an  ionic  compound  (NaCl)  and  it  is  found  to 
give  surprisingly  good  agreement  with  self-consistent  cal¬ 
culations. 

Our  fourth  major  approximation  involves  the  evalua¬ 
tion  of  the  matrix  elements  of  the  l^arious  terms  that 
make  up  h(n®).  The  kinetic  energy,  overlap,  and  nonlo¬ 
cal  part  of  the  pseudopotential  are  easily  evaluated  exact¬ 
ly  (numerically)  for  each  geometry.  A  table  of  results  on 
a  fine  grid  of  separations  is  constructed  so  that  the  matrix 
elements  for  any  separation  between  the  atoms  is  accu¬ 
rately  interpolated.  The  three  center  matrix  elements  of 
the  neutral-atom  potential  given  by 

)'NA{r-r„)-Kte„(r-r,)+?7~~^dV  (11) 

are  accurately  appro.ximated  by  an  r-dependent  expan¬ 
sion  in  multipole  moments.  The  nonlinear  exchange- 
correlation  matrix  elements  are  calculated  within  the 
average  density  approximation  described  in  Ref.  23.  All 
matrix  elements  are  calculated  in  real  space. 

The  forces  are  determined  by  differentiating  the  total 
energy  £}^’  given  by  Eq.  (10),  Fy*"-8£jj,'/9ry.  The 
most  difficult  term  comes  from  the  band  structure.  Its 
derivative  is 


The  Schrodinger-like  Eq.  {4a)  follows  from  a  variation¬ 
al  principle  for  the  total  energy  of  Eq.  (1).  The  single¬ 
particle  Hamiltonian,  Eq.  (S),  defines  a  problem  requiring 
a  self-consistent  solution.  To  avoid  the  difficulties  associ¬ 
ated  with  iterating  to  self-consistency,  we  adopt  an  ap¬ 
proximation  suggested  by  Harris.*^  ITjus  our  third  major 
approximation  is  to  consider  a  sum  of  neutral-atom 
spherical  charge  densities  as  a  zero-order  approximation 
to  the  self-consistent  density  of  the  cluster,  and  to  keep 
only  first-order  changes  from  this  density  in  the  energy 
functional.  Thus  we  write 

«(r)=«°(r)+5rt(r) ,  (8) 


“2  2 
I 

(ncc) 


3r, 


•r  /i.v 

a.a* 


-“2 


*4,V 

a»<i 


d(h°C' 

dr, 


fart' 

ftv 


(12) 


where  p®?'  and  are  the  density  and  energy-density 
matrices,  respectively. 
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/ 

(occ) 

■^“v  =  2  £,a*(^,a)o,{v,a') , 

i 

iocc) 

where  o,(^,a)  are  the  expansion  coefficients  of  the  wave 
function  in  Eq.  (2). 

After  the  forces  arc  evaluated,  the  equations  of  motion 
can  be  solved  and  the  positions  and  velocities  of  the 
atoms  updated.  The  energy  functional  has  been  tested  in 
bulk  Si  and  found  to  give  bulk  static  properties  that  arc 
in  good  agreement  with  experimcnt.^^  It  is  the  purpose 
of  this  paper  to  extend  these  calculations  to  the  static  and 
dynamic  properties  of  small  Si  clusters  using  an  ab  initio 
molecular-dynamics  technique. 

III.  MOLECULAR  DYNAMICS 

The  equilibrium  structures  and  vibrational  spectra  of 
Si  clusters  have  been  determined  using  the  technique  of 
molecular  dynamics.  The  force.s  are  calculated  quantum 
mechanically  from  the  total  energy  [Eq.  (10)]  and  its 
derivatives  [c.g.,  Eq.  (12)].  The  nuclear  coordinates  arc 
moved  in  time  according  to  the  classical  equations  of 
motion,  F„  d-T/di^,  These  equations  are  integrated 
using  the  Gear  predictor-corrector  algorithm,^  with  a 
time  step  of  ~  1.5  fs. 

The  ground-state  equilibrium  structures  are  found  by 
one  of  two  techniques;  dynamical  quenching  or  simulat¬ 
ed  annealing.  In  either  of  these  two  techniques,  the  clus¬ 
ter  is  started  in  some  trial  structure,  generally  with  very 
little  symmetry. 

In  the  technique  of  dynamical  quenching,  the  atoms  in 
the  cluster  are  allowed  to  respond  to  internal  forces  and 
are  accelerated.  A  kinetic  temperature  Ti;  is  defined  as 
the  average  classical  kinetic  energy  of  the  atoms, 

i/2kT^=[]/N)2iymiV^,  where  /  =  1,2 . n  is  the 

atom  indc.x.  As  the  atoms  accelerate,  the  kinetic  temper¬ 
ature  increases  until  a  maximum  is  reached.  The  system 


1  2  3  4  5  6  7  8 

dusts,  size  (n)  buir  — ► 


FIG.  2.  The  energy  dillerencc  per  atom  between  Si  tluMcrs  of 
size  =2-7  and  bulk  S  The  solid  dots  correspond  to  emiilibn- 
um  siructure.s  e.splaiiii  in  the  text,  and  oilitr  nu’i.oi.ible 
configurations  are  also  shown  for  ii  =4-7. 


duyLef  siz«  ii 

FIG.  3.  The  n-dependent  fragmentation  energy  (A£f„j)„  as  a 
function  of  n.  The  dip  at  n=S  indicates  that  fragmentation  is 
more  likely  to  occur  for  this  cluster. 


is  then  quenched  by  setting  all  the  velocities  to  zero,  thus 
removing  all  the  kinetic  energy.  The  atoms  again  are  al¬ 
lowed  to  accelerate  and  the  quenching  process  is  repeat¬ 
ed.  After  several  dynamical  quenching  cycles,  a 
minimum  energy  configuration  is  obtained.  The 
minimum  may  be  only  a  local  minimum, that  the  pro¬ 
cedure  needs  to  be  repeated  several  times  with  different 
initial  configurations. 

In  the  simulated  annealing  technique,  the  energy  is  re¬ 
moved  more  gradually  so  that  a  gl  'bal  mi:.'  '  .um  is  more 
likely  to  be  obtained.  The  technique  iias  been  described 
by  Kirkpatrick  ei  al}^  and  uses  the  Monte  Carlo  algo¬ 
rithm  of  Metropolis  et  Monte-  Carlo  steps  are  taken 
where  the  coordinates  of  each  particle  are  changed  as 
Tj-'r, +6r/.  Here  6r,=/-o6R,  wlide  5R  is  a  triplet  of 
random  numbers  between  -  1  and  1  for  each  atom  /, 
and  To  is  the  maximum  step  size.  .\n  energy  ..'iffercnce 
6£  =  — F,.-6r,  is  computed  and  the  Is  taken  if  8£  <0, 
while  if  6£>0  the  step  is  taken  with  probability 
e  where  is  an  annealing  temp-erature.  This 

ensures  that  after  many  moves,  the  en.semble  lends  to  the 
Boltzmann  distribution.  The  lemperature  is  gradual¬ 
ly  reduced  so  that  the  system  settles  into  the  ground 
state.  We  have  found  that  ihe  dynamical  quenching  finds 
the  ground  state  much  more  .'cadily  than  simulated  an¬ 
nealing  for  these  small  clusle:.^.  However,  simulated  an¬ 
nealing  is  more  likely  to  find  the  true  ground  state  than 
dynamical  quenching.  In  eii.ier  case,  however,  a  variety 
of  initial  configurations  musi  be  tried  before  one  can  be¬ 
lieve  that  the  true  ground  state  has  been  found. 

The  vibrational  spectrum  i^  found  directly  from  the 
molecular-dynamics  simulatioi;  without  annealing  or 
quenching.  We  start  in  a  groui'd-siat,.  configuration  and 
have  constrained  the  motion  s-i  that  the  cc  iter  of  mass  of 
the  system  of  atoms  is  stationary  and  there  is  no  angular 
momentum  about  the  cer.iei  of  mass.  To  accompli.sh 
this,  we  first  give  the  atoms  random  velocities  commensu- 
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rate  with  a  given  initial  kinetic  temperature.  In  general, 
however,  the  angular  momentum  (L)  is  nonzero.  To 
keep  the  system  from  rotating,  we  adjust  the  velocities  by 
subtracting  the  quantity  v,  ={l~'*L)Xr,  from  the  veloci¬ 
ty  vector  Vj  (i  =  1,2, . . . , « ),  where  I  is  the  inertia  tensor. 
This  procedure  removes  the  angular  momentum  while 
still  allowing  for  randomness  in  the  velocities.  Finally 
the  velocities  are  uniformly  shifted  so  that  there  is  no 
center-of-mass  motion. 

The  frequency  spectrum  is  determined  by  Fourier 
transforming  the  velocity  autocorrelation  function, 

*  <v,(0)-v,{0))  • 

Here  n  is  the  atom  index  and  the  brackets  {  •  •  •  )  indi- 


d(A) 


FIG.  4.  (a)  The  one-electron  LDA  energy  eigenvalues  vs  sep¬ 
aration  distance  d  for  Sij.  The  doubly  degenerate  (neglecting 
spin)  lir,  level  crosses  the  singly  degenerate  la,  level  at  about 
d=2. 12  A.  (b)  Total  energy  per  atom  vs  separation  distance  d 
for  Sij.  A  secondary  minimum  occurs  because  of  the  level 
crossing  shown  in  (a). 


d  =  2.189  A 


FIG.  5.  The  ground-state  configuration  of  Sis. 


9  (deg) 


FIG.  6.  (a)  The  one-electron  LDA  energy  eigenvalues  versus 
opening  angle  6  fer  Sij,  where  the  bond  lengths  d  are  kept  con¬ 
stant  at  their  equilibrium  value  (2.189  A).  Note  the  level  degen¬ 
eracies  which  occur  at  9=60*,  which  is  the  highly  symmetric 
case  of  an  equilateral  triangle.  The  Fermi  level  separating  the 
occupied  from  the  unoccupied  levels  is  s.heniatic.illy  shown,  (b) 
Total  energy  per  atom  vs  opening  angle  0  for  Sij  ill-istrating  the 
bond  bending  forces  of  the  Si]  molecule. 
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FIG.  7.  Spectral  density  g{a)  of  Si3  in  its  equilibrium 
configuration  at  a  low  (lower  figure)  and  a  high  (upper  figure)  ki¬ 
netic  temperature.  The  kinetic  temperatures  are  50  and  500  K, 
respectively.  The  modes  soften  at  higher  excitation  levels,  and 
there  is  added  structure  due  to  anharmonicity. 


w  (can"^) 


FIG.  8.  (a)  A  snapshot  of  the  planar  motion  of  an  Si3  mole¬ 
cule  at  high  excitation  corresponding  to  a  kinetic  temperature 
of~25(X)K.  The  chaotic  motion  very  much  anharmonic.  (b) 
The  spectral  density  function  gfw)  for  the  motion  in  (a).  The 
three  peak  structure  is  lost  and  the  spectrum  becomes  continu¬ 
ous. 


FIG.  9.  The  ground-state  configuration  of  Si4. 


cate  that  an  ensemble  average  over  all  atoms  is  taken. 
This  ensemble  average  is  defined  by 

-  )fK ) , 

tn 

where  tj  is  the  time  at  the  ;th  stei)  and  j  =0, 1,2 . M, 

with  M  being  the  total  number  of  time  steps  in  the  simu- 


Q  (deg) 


Fill.  10.  (a)  The  one-elcctron  LDA  energy  eigenvalues  vs  in¬ 
terior  angle  6  for  Sij.  In  the  ground  .state,  9=62'.  The  side 
length  was  kept  constant  at  its  ground-state  value.  Degenera¬ 
cies  and  Jahn-Teller  instabilities  occur  at  90'.  (b)  Total  energy 
per  atom  vs  interior  angle  Q  for  Sij. 
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FIG.  11.  Spectral  density  ^((u)  of  Si4  in  its  equilibrium 
configuration  for  high  and  low  kinetic  temperatures.  The  peaks 
in  the  spectrum  are  labeled  according  to  the  symmetry  of  the 
displacements. 


lation.  The  spectral  density  glru)  is  the  Fourier  cosine 
transform  of  the  velocity  autocorrelation  function  gU), 

g(co)=(l/r)fjg(tmi}coslat)di  , 

where  T  is  the  total  time  of  the  simulation.  iy(i}  is  the 
Blackman  window  function^^  used  to  reduce  oscillations 
due  to  finite  time  sampling. 

IV.  RESULTS  FOR  Si  CLUSTERS 

In  this  section  we  individually  discuss  the  results  of 
clusters  from  Si,  with  «  =2-7.  We  begin  by  summariz¬ 
ing  the  energetics  of  the  various  geometries  found 
by  quenching  and  annealing.  In  Fig.  2  we  show  the 
energy  difference  per  atom,  =(£„,, /atoml^iu,,,, 

/atomlfcuii;,  between  bulk  (diamond)  Si  and  the  Si 
clusters  in  their  ground  state  and  metastable  equilibrium 
configurations. 

The  solid  circles  indicate  the  ground-state 
configurations  (to  be  described  individually),  while  the 
open  circles,  squares,  and  triangles  represent  local- 
minimum  (metastable)  structures.  The  energy  per  atom 
of  the  cluster  approaches  the  bulk  binding  energy  mono- 
tonically,  but  there  is  a  marked  decrease  in  slope  in  going 


TABLE  I.  Comparison  of  the  harmonic  frequencies  (in 
cm”')  of  the  Su  molecule  in  this  work  with  quantum-chemistry 
SCF  results. 


Mode 

symmetry 

This  work 

6-31  G* 

(Ref.  43) 

6-31  G 

116 

130 

137 

Bl, 

243 

223 

157 

d, 

346 

380 

328 

464 

453 

371 

B,. 

495 

543 

472 

A. 

546 

503 

425 

FIG.  12.  Ground-state  configuration  of  Si5.  This  structure  is 
a  trigonal  bipyramid,  consisting  of  an  equilateral  triangle  with 
one  atom  above  the  plane  of  the  triangle  and  one  atom  below. 


from  n~4  to  n=5.  This  can  be  seen  more  clearly  by 
considering  the  fragmentation  energy  as  defined  by  Raga- 
vachari  and  Logovinsky^^  according  to  the  reaction 
Si^-t-Si^-i+Si.  The  energy  cf  this  reaction  is 


where  (A£f„j),=(n  — l)A£,_j+nA£,.  The  last  term 
in  £q.  (13)  is  a  constant  (independent  of  n).  This  term  in¬ 
volves  the  energy  of  an  isolated  atom,  which  is  not  well 
described  in  our  method  because  of  spin  effects  and  the 
use  of  compact  orbitals.”  The  trends  with  n  of  (£(nf}n 
are  entirely  contained  in  the  first  term  (A£f„f ),  that  is 
plotted  in  Fig.  3.  The  lesser  stability  of  Sij  compared  to 
neighboring  clusters  is  clearly  evident.  The  results  are  in 
very  close  agreement  with  the  fragmentation  trends 
found  by  Raghavachari.” 


Ci?  (cm  •' ) 

FIG.  13.  Spectral  density  gtoi)  of  Si5  in  its  equilibrium 
configuration  at  low  kinetic  temperature  (T*® 50  K). 
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FIG.  14.  Ground-state  configuration  of  Si^.  This  structure  is 
a  square  bipyramid,  consisting  of  a  square  with  one  atom  above 
the  plane  of  the  square  and  one  atom  below. 


A.  Sii 

Dynamic  quenching  of  Sij  gives  the  equilibrium  bond 
distance  at  d=2.27  A,^in  good  agreement  with  the  exper¬ 
imental  value  of  2.24  A.”  The  bulk  bond  length  of  crys¬ 
talline  .Si  is  2.35  A,  while  the  present  theory  gives  2.38  A. 
Thus  we  find  especially  good  agreement  between  theory 
and  experiment  for  the  large  reduction  in  nearest- 
neighbor  distance  in  going  from  bulk  Si  to  Sij  (0.11  A  ex¬ 
periment  and  0.1 1  A  theory). 

The  single  harmonic-vibrational  stretch  mode  is  deter¬ 
mined  to  be  o)=531  cm"',  which  compares  well  to  the 
experimental  result  of  511  cm"'.^’  In  this  case  at  least, 
our  results  are  comparable  to  the  considerably  more  so¬ 
phisticated  quantum  chemistry  self-consistent-field  (SCF) 
results  of  Ref.  36  who  find  568  cm"'. 

In  Fig.  4(a)  are  shown  the  eigenvalues  of  the  one- 
electron  LDA  Hamiltonian  for  Sij  as  a  function  of  Si-Si 
distance.  The  Im^  electronic  level  is  fourfold  degenerate 
(including  spin)  and  contains  two  electrons  in  equilibrium 
(d=2.27  A).  The  2a ^  level  crosses  the  l7r„  level  as  the 
distance  is  decreased.  This  is  in  excellent  agreement  with 


FIG.  15.  Ground-siaie  configuration  of  Sij.  This  structure  is 
a  pentagonal  bipyramid,  consisting  jf  a  pentagon  with  one  atom 
above  the  plane  of  the  pentagon  and  one  atom  below. 


the  spin-polarized  results  of  Northrop  ei  al.,^*  who  find 
that  the  l7r„  level  crosses  the  doubly  degenerate  2ff,  level 
at  about  4.00^  (2.12  A),  resulting  in  a  new  ground-state 
configuration  at  smaller  separations.  We  also  find  the 
lrr„  level  crossing  the  2a,  level  at  d=  2. 12  A.  This  re¬ 
sults  in  a  potential-energy  surface  with  two  minima  as 
shown  in  Fig.  4(b).  A  metastable  configuration  exists  at 
<f=2.03  A. 

B.  Sij 

The  Sij  ground-state  configuration  we  find  by  dynami¬ 
cal  quenching  is  an  isosceles  triangle  with  the  two  equal 
sides  of  length  d=2. 189  A  and  an  opening  angle  of 
0=78.8*  (see  Fig.  5).  These  results  are  in  good  agreement 
with  the  more  rigorous  quantum-chemistry  SCF  calcula¬ 
tions  of  Diercksen  et  al.^^  who  find  d=2. 196  A  and 
0=80.6*.  Grev  et  al.*^  who  find  d=2. 160  A  and 
0=78.1*,  and  Raghavachari  et  al?*  who  find  d—lAT 
and  0=77.8*. 

The  energy  eigenvalues  as  a  function  of  the  angle  0 
(with  the  distance  d  kept  constant  at  its  value  at  equilibri¬ 
um)  are  shown  in  Fig.  6(a).  The  Fermi  level  lies  between 
the  sixth  and  seventh  levels.  These  two  levels  become  de¬ 
generate  at  0=60*  (an  equilateral  triangle),  which  make 
this  geometry  Jahn-Teller  unstable.  The  total  energy  as  a 
function  of  0  (again  with  d  fixed)  is  shown  in  Fig.  6(b). 
Two  equivalent  minimum  are  seen  corresponding  to 
0=78.8*  and  0=281.2*.  The  system  is  ijlearly  unstable  at 
60*,  and  a  relatively  large  barrier  separates  the  two 
equivalent  minimum  as  the  molecule  passes  through  180*, 
where  the  molecule  forms  a  linear  chain  as  occurs  for  Cj. 
Curiously,  we  find  a  very  shallow  local  minimum  at  the 
linear  chain.  This  is  to  be  contrasted  to  the  SCF  results 
of  Raghavachari,'^  which  find  the  linear  chain  to  be  un¬ 
stable.  The  imaginary  frequency  found  in  that  work  (82/ 
cm"')  is  very  small,  indicating  a  relatively  flat  potential 
surface  similar  to  ours. 

The  spectral  density  function  g(a))  in  the  harmonic 
tr^-  =  50  K)  and  the  near-harmonic  (Tf;^i00  K)  regions 
for  Sij  vibrations  arc  show.n  in  Fig.  7.  Three  peaks  arc 
apparent  corresponding  to  the  3w  -6  normal  modes  of  vi¬ 
bration.  The  width  of  the  peaks  is  due  to  the  finite  time 
of  the  simulation.  The  position  of  the  peaks  indicates  the 
frequency  of  the  normal  mode,  while  the  relative  heights 
reflect  the  (randomly  chosen)  relative  amplitudes  in  each 
of  the  modes.  The  spectrum  at  low  temperature  is  pure, 
indicating  very  little  anharmonicity,  while  at  higher  tem¬ 
peratures,  small  contributions  to  anharmonic  affects  (sum 
and  difference  frequencies)  are  evident.  Also  notice  the 
overall  shift  toward  lower  frequencies.  The  “harmonic" 
frequencies  at  50  K  are  202,  474,  ami  596  cm"'  for  the 
A,,  Bj,  and  A^  symmetric  modes,  respectively.  These 
can  be  compared  with  those  obtained  by  Grev  ci  of 
157,  570,  and  574  and  Raghavacliari^'  of  206,  560,  and 
582.  We  note  that  both  of  our  A  |  modes  are  in  good 
agreement  with  these  quantum-chemistry  SCF  calcula¬ 
tions,  while  the  Bi  mode  is  only  in  fair  agreement,  with 
our  mode  lying  ~  I67<j  lower  than  their.v  This  is  not  un- 
reason.nble  as  typically  calculated  frequencies  of  mole¬ 
cules  are  accurate  to  — 10- 15 
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The  molecular-dynamics  technique  makes  no  assump¬ 
tions  about.harmonic  potentials  and  does  not  construct  a 
dynamical  matrix.  Very  harmonic  systems  are, treated  in 
exactly  the  same  way  as  an  anharmonic  system.  We 
show  in  Fig.  S  the  results  of  a  simulation  of  Sij  at  a  very 
high  temperature  of  ~2500  K.  The  simulation  was  of 
.1600  time  steps  and  there  is  no  indication  that  the  mole¬ 
cule  is  about  to  disassociate,  'i  he  motion  is  two  dimen¬ 
sional  and  a  snapshot  of  the  motion  is  shown  in  Fig.  8(a). 
Within  the  chaotic  motion,  one  sees  a  nearly  sixfold  pat¬ 
tern  c.iming  from  a  molecule  containing  only  three 
atoms.  This  occurs  because  the  atom  that  acts  as  the 
apex  of  the  ground-state  triangular  molecule  is  being 
changed  from  one  atom  to  'he  next.  Also  there  is  consid- 
erabie  classical  tunnelingli’^e  behavior  of  the  apex  atom 
through  the  center  line  connecting  the  other  two  atoms. 
The  spectral  density  function  g(a>)  (Fig.  8(b))  has  lost  the 
three-peak  structure  of  the  harmonic  system,  and  takes 
on  a  continuous  spectrum.  The  frequencies  depend  on 
the  amplitude  of  each  cf  the  “modes,"  and  the  ampli¬ 
tudes  are  continuously  changing. 

C.  Si4 

The  equilibrium  configuration  determined  for  Si;  is  a 
rhombus  with  side  length  2.32  A  and  a  minor  diagonal 
length  of  2.39  A  'tjc  Fig.  9'.  This  compares  very  well 
with  the  results  of  Tomanek  and  Schliiter,'*  who  find^a 
rhombus  of  side  length  2.3  A  with  a  diagonal  of  2.4  A, 
and  the  results  of  Raghavachari  and  Logovinsky,^  who 
also  find  a  rhombus  of  side  2.30  A  and  diagonal  2.40  A. 
Other  possible  geometries  include  the  square  and 
tetrahedron.  The  square  forms  a  metastable  minimum- 
energy  configuration  with  a  side  of  length  2.32  A,  and  the 
tetrahedron  is  metastable  with  minimum  energy  of  side 
length  2.46  A.  The  energies  of  these  two  structures  are 
shown  in  ‘.-ig.  2  and  are  significantly  higher  in  energy 
(=s0.61  eV.-'atom  and  =*0.53  eV/atom,  respectively)  than 
the  rhombus. 

The  one-electron  LD.A  energy  eigenvalues  are  plotted 
as  a  function  of  the  angle  6  in  Fig.  10(a).  As  in  the  case 
of  Si},  there  is  again  a  crossing  of  levels  at  the 
configuration  of  highest  symmetry  (in  this  case,  0=90*). 
This  causes  a  Jahn-Teller  unstable  system  that  shows  up 
as  a  cusp  in  the  total  energy  per  atom  also  as  a  function 
of  the  angle  0  in  Fig.  10(b). 

The  vibrational  spectrum  of  Si;  is  shown  in  Fig.  1 1  for 
a  high  and  low  kinetic  temperature.  There  arc  3»  — 6=6 
normal  modes  for  Si;  and  the  high-temperature  spectrum 
shows  the  mode-softening  characteristic  of  anharmonici- 
ty  at  higher-excitation  levels.  The  normal  coordinates 
that  give  rise  to  the  various  peaks  in  the  spectrum  are  la¬ 
beled  according  lo  symmetry.  The  lowest  frequency 
mode  is  £}„,  which  is  the  “bond-bending"  mode  in  which 
alternate  atoms  move  on:  of  and  into  the  plane  of  the  pa¬ 
per  (see  Fig.  9'.  The  other  bond-bending  mode  is  flju  in 
which  the  atoms  move  in  the  plane  of  the  paper.  The  fre¬ 
quencies  found  from  an  examination  of  Fig.  1 1  are  com¬ 
pared  with  those  of  Raghavichari  and  R.ohlfing'''  in  Table 
I.  Overall  agreement  with  these  far  more  sophisticated 
calculations  is  excellent. 


D.  Si, 

We  find  the  equilibrium  structure  of  Si,  based  on  an¬ 
nealing  and  quenching  to  be  a  trigonal  bipy^ramid,  con¬ 
sisting  of  an  equilateral  triangle  of  side  3.02  A  capped  on 
top  and  bottom,  with  the  distance  between  the  corners  of 
the  triangle  and  a  cap  being  2.30  A  (see  Fig.  12).  This 
structure  agrees  with  tho.se  found  by  Tomanek  and 
Schluter,*'*  who  find  a  trigonal  bipyramid  with  the 
relevant  lengths  of  3.1  and  2.4  A,  and  Raghavachari  and 
Logovinsky,^'*  who  find  the  same  structure  with  distances 
of  2.34  and  3.26  A.  Metastable  structures  for  Si,  include 
the  pentagon  and  pyra.mid.  The  energetics  of  these  struc¬ 
tures  are  shown  in  Fig.  2,  where  the  pentagon  is  =*1.00 
eV/atom  higher  and  the  pyramid  is  ==0. 13  eV/atom 
higher. 

The  vibrational  spectral  density  function  of  Si,  is 
shown  in  Fig.  13  at  low  temperature  ( T = 50  K).  The  five 
peaks  in  the  figure  come  about  because  three  frequencies 
are  doubly  degenerate,  while  an  additional  frequency  is 
“accidentally”  degenerate,  and  unresolvable. 

E.  Si; 

The  equilibrium  structure  for  Sij  is  that  of  a  bipyram¬ 
id,  consisting  of  a  square  with  side  length  2.71  A  capped 
on  top  and  bottom,  witli  the  distance  between  a  vertex  of 
the  square  and  a  cap  being  2.36  A  '.see  Fig.  14).  This 
structure  agrees  with  that  found  oy  Tomanek  and 
Schliiter,'^  who  also  find  such  a  "distorted  octahedron" 
with  the  relevant  distances  of  2.6  and  2.3  A.  Raghava- 
chari,’^  however,  finds  the  ground-state  structure  to  be  an 
edge-capped  trigonal  bipyramid. 

The  energies  of  higher-energy  metastable  states  for 
n=6  are  shown  in  Fig.  2.  Three  structures  investigated 
arc  structures  higher  in  energy  than  the  equilibrium 
structure:  (1)  a  hexagon  ( =*  1.22  cV/atorn  higher  in  ener¬ 
gy),  (2)  the  singly  capped  pentagon  (=*0.50  eV/atom 
higher),  and  (3)  a  sixfold  buckled  ring  which  is  a  bulk 
fragment  ( =*0.22  eV/atom  higher). 

F.  Si, 

The  ground-state  structure  found  for  Si}  is  a  bicapped 
pentagon,  consisting  of  a  pentagon  of  dde  length  2.46  A 
capped  on  top  and  bottom,  with  the  distance  between  a 
vertex  of  the  pentagon  and  a  cap  being  2.45  A  (see  Fig. 
15).  This  is  in  agreement  with  Ballone  el  al.,**  who  also 
find  a  pentagonal  bipyramid  as  the  equilibrium  structure. 
In  addition,  when  just  dynamical  quenching  is  performed, 
the  system  invariably  settles  into  a  capped  and  strongly 
reconstructed  version  of  the  sixfold  buckled  ring.  This 
metastable  structure  is  relatively  close  in  energy  ( =*0. 16 
eV/atom  higher)  to  this  equilibrium  structure. 

V.  CO.N'CLUSIONS 

We  have  used  an  approximate  electronic-structure 
method  to  determine  the  electronic  and  vibrational  spec¬ 
trum  and  the  equilibrium  structures  of  small  Si  clusters. 
The  method  used  is  a  simplified  ab  initio  tight-binding 
model  based  on  local-density  thv>iry,  and  has  no  adjust- 
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able  parameters.  We  have  compared  a  number  of  our  re¬ 
sults  with  those  of  others  and  find  substantial  agreement 
with  more  rigorous  calculations.  Based  on  these  results 
for  small  molecules,  we  are  optimistic  that  the 
molecular-dynamics  technique  and  energy  functional 
used  here  will  be  applicable  to  bulk  crystalline  and  amor¬ 
phous  covalent  semiconductors  and  to  semiconductor 
surfaces  as  well.  The  method  is  computationally  fast  and 
efficient,  and  should  permit  us  to  study  covalent  systems 
of  previously  incalculable  complexity. 
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An  inverse  formula  is  given  for  a  new,  generalized  Mdbius  transform  C(.r.;\...)»Sr»iS5.i-F’{«*X/»r|f....).  namely 
F(.v.  .f. ... )  ■  Ir.  I S  5. 1  ...pin  )p(m  )G(  m‘'.v.  m>y. ...  1.  In  the  case  of  one  dimension  and  aa  - 1 .  this  reduces  to  the  Chen-M8biu$ 
transform  which  has  been  applied  to  the  inverse  black*body  radiation  problem.  For  a®  + 1  this  becomes  the  Hardy-Wright- 
Mubius  transform. 


Recently  Chen  ( 1  ]  modified  Mdbius’  number- 
theory  transform  [2]  and  showed  how  it  could  be 
used  to  solve  a  class  of  inverse  problems  in  physics, 
including  the  problem  of  determining  the  tempera¬ 
ture  distribution  of  a  black  body  from  its  radiation 
spectrum:  given  a  quantity  F{x)  that  is  related  to  a 
function  G{x)  by 

(7(a-)=  I  f(.Y/rt).  (la) 

fim  1 

the  inverse  function  F  is  given  by 

Pix)=  I  M»)  G{x/n) ,  (lb) 

ttw  1 

where /t(n)  is  the  Mdbius  function: /tCft)*!  1  ifn=l. 
;i(n)=!(- 1 )'  if  n  is  the  product  of  r  distinct  prime 
factors,  and  jt  ( « )  =  0  otherwise  ( 2 ) .  For  example,  in 
the  physics  of  black  bodies.  G{i')  may  be  related  to 
the  radiation  spectrum  and  F{  T)  may  be  related  to 
the  areal  thermal  distribution  of  the  body  (1  ].  As 
shown  by  Chen,  this  theorem  is  particularly  useful 
for  inverting  data,  for  example,  obtaining  the  un¬ 
known  areal  distribution  of  tcmiviratures  on  a  black 
body  from  its  radiation  spectrum. 

Here  (i)  we  combine  Chen's  result  with  another 
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inverse  Mdbius  transform  from  number  theory  [  2  ] *' 
to  obtain  a  more  general  inverse  formula,  and  (ii) 
we  extend  this  generalized  formula  to  two  and  higher 
dimensions.  The  new  and  more  general  inverse 
Mdbius  irr'.nsform  formula  for  one  dimension  is: 

If  we  have 

G(,x)^  I  ^("“v)  (2a) 

4ai  I 

then  it  is  the  case  that 

Fl\)vr.  X  n{n)G{n«.K) .  (2b) 

nm  1 

Hardy  and  Wright  proved  thts  theorem  for  aal 
("theorem  270”  [2]  ”)  and  Chen's  Mdbius  theo¬ 
rem  corresponds  to  a »  —  1.  We  show  that  the  theo¬ 
rem  is  valid  for  any  real  a,  and  hence  there  are  entire 
classes  of  Mdbius  transforms,  each  corresponding  to 
a  different  choice  of  a.  Our  proof  uses  “theorem  263” 
[2]  for  the  Mdbius  function, 

J]jt(rt)a«l,  forAratl, 

=0.  otherwise.  (3) 

Here  the  sum  is  over  ail  the  different  divisors  n  of 
the  integer/:,  including  1  and  A’.  The  proof  of  the  gen- 

•'  In  fact,  due  10  cq.  (3l)ofrcf.  (I].cq.  (32)i'>ref.  [1]  i$a 
direct  result  of  theorem  270  of  ref.  (2). 
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eralized  formula  is  obtained  by  evaluating 
f  ;i(/j)G(«"v)=  f  f  /i(«)f(m“n“v) 

nml  nmltuml 

=  f  m-“A-)I/t(n)=f(.v).  (4) 

km\  n\k 

In  this  derivation,  we  performed  the  sum  over  k= inn 
and  n,  instead  of  over  m  and  /i;  this  reordering  will 
converge  if 

I  I  \F(m«n«x)\=.ld(k)\F{lc‘‘x)\  (5) 

m  St  i  /I  M  1  k 

converges,  where  d{k)  is  the  number  of  divisors  of 
A-  [2].  Note  that  for  asO  this  convergence  criterion 
is  normally  not  met,  and  hence  o!=:0  is  to  be  avoided. 
Furthermore  the  converse  of  eq.  (2)  is  true  if 

I  l\G(m<>n^x)\=rld{k)\G{k‘‘x)\  (6) 

n  «  t  >N  ss  1  k 

converges. 

Generalization  of  formula  (2)  to  higher  dimen¬ 
sions  is  straightfonvard;  for  example,  in  two  dimen¬ 
sions  for  F{x,  y)  of  the  form 

,1')  =  I  f  ,  (7a) 

tnm  I  fl«  I 

we  have 

nx,  ,v)  =  I  I  H{k).u(l)  Gik-x,  1“}’) .  ( 7b ) 

ill  /•ii 

■A  sufficient  condition  for  this  inverse  formula  to  be 
valid  is  the  convergence  of 

k  in  i  n 

^ZldV)dU)\FU'’x,/y)\.  (8) 

•  J 

The  new  inverse  Mdbius  formula  (2)  can  be  applied 
to  problems  in  mathematics  and  physics,  and  a  con- 
venii*ni  choice  of  the  indices  ot,  ...  can  be  selected 
to  best  suit  each  problem. 

The  one-dimensional  version  leads  to  the  new 
mathematical  identity 


■—=  f.  y{n)n-''  forp>l.  (9) 

■.(P)  n-l 

where  C(p)  =  I,f.i is  the  Riemann  zeta  func¬ 
tion.  It  also  leads  to  the  identities 

I  i/i{«)((tt-v/«)/=(it-v/rt)-l]  (10) 

X  M  i  ntm  i 

and 

Mx)=  f  +  i 

nmi  \  - 

-1-2  i  {n-x--4m-n-)-''-),  (11) 

where  we  have  lNn<nx<  (2.V+2)jt.  Jq  is  the  cylin¬ 
drical  Bessel  function.  /,.(.v)=coth(.').  and 
/_(.v)=scot(.v).  and  we  have  used  the  identities  [3) 

/-(>«•■)= -^.-1-^  I  (.v-±rt-)-'. 

7t  /iM  I 

and 

I  yo(n.v)*-:^-l--.+2  I  (.v--4;n=it')-''2 
/t«i  -  A  mv) 

for  2A'71<a:<  (2;V+2)7c. 

There  are  undoubtedly  many  more  u.cs  for  this 
inverse  formalism  in  both  mathematics  and  theo¬ 
retical  physics  that  will  be  uncovered  as  it  is  used 
more. 
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